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1 12.  SPONSORING  MILITARY  ACTIVITY 


U.  S.  Army  Materiel  Command 


The  operation  of  a  f lueric  a  -<rtex  rate  sensor  is  discussed  for 
flows  in  the  fully  developed  range,  the  ranqe  in  which  the  entrance 
conditions  at  the  periphery  of  the  chamber  do  not  appreciably  in¬ 
fluence  the  radial  velocity  profile  throughout  most  of  the  chamber. 
Operation  in  this  range  is  far  different  from  that  in  the  small 
boundary  layer  range  usually  considered  in  previous  investigations. 

An  intuitive  examination  shows  that  the  frequency  response  and  rate 
threshold  of  the  sensor  are  inferior  in  the  smal3 -boundary -layer 
range,  chosen  originally  for  its  desirable  gain  features. 

The  Navier-Stokes  equations  as  they  apdy  to  flow  in  the  chamber 
for  the  fully  developed  range  are  simplified  through  very  sliqht 
approximations,  and  a  method  of  solving  the  resulting  equations  is 
outlined.  A  good  approximation  of  this  lengthy  method  is  used  to 
obtain  closed  form  expressions  that  predict  the  gain  and  freauencv 
response  of  the  sensor.  These  predictions  show  close  oualitative 
agreement  with  data  obtained  from  a  model  of  the  present  design. 

The  only  important  source  of  error  originates  from  the  viscous  losses 
in  the  sink  region,  which  are  not  fully  analyzed.  However,  it  is 
shown  that  this  error  does  correlate  with  the  go’/erning  parameters, 
so  that  the  performance  can  be  predicted  to  within  about  10  percent 
once  a  single  experiment  to  measure  the  losses  in  the  sink  is 
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conducted  on  a  given  pickof f-drain  design.  This  is  verified  with 
several  experiments. 

The  results  predict  that  the  frequency  response  of  the  sensor 
can  be  increased  by  about  an  order  of  magnitude  over  that  of  the 
present  designs  without  reducing  sensitivity.  Although  the  flow 
field  in  the  fully  developed  range  could  be  less  noisy  than  in 
the  boundary -layer  range,  the  noise  level  of  the  present  test 
model  is  not  reduced  in  the  fully  developed  range  because  of  its 
design . 

Final  discussion  covers  design  changes  that  should  signifi¬ 
cantly  reduce  the  noise  level  and  make  the  sensor  more  rugged. 

It  also  mentions  planned  future  analytical  work  aimed  at  de¬ 
veloping  a  sensor  that  compensates  for  environmental  changes. 
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ABSTRACT 


The  operation  of  a  flueric  vortex  rate  sensor  is 
discussed  for  flows  in  the  fully  developed  range,  the  range 
in  which  the  entrance  conditions  at  the  periphery  of  the 
chamber  do  not  appreciably  influence  the  radial  velocity  pro¬ 
file  throughout  most  of  the  chamber.  Operation  in  this  range 
is  far  different  from  that  in  the  small-boundary-layer  range 
usually  considered  in  previous  investigations.  An  intuitive 
examination  shows  that  the  frequency  response  and  rate  thresh 
old  of  the  sensor  are  inferior  in  the  small-boundary-layer 
range,  chosen  originally  to  obtain  maximum  gain. 

•  The  Navier-Stokes  equations  as  they  apply  to  flow 
in  the  chamber  for  the  fully  developed  range  are  simplified 
through  very  slight  approximations,  and  a  method  of  solving 
the  resulting  equations  is  outlined.  A  good  approximation 
of  this  lencthy  method  is  used  to  obtain  closed  form  expres¬ 
sions  that  predict  the  gain  and  frequency  response  of  the 
sensor.  These  predictions  show  close  quantitative  agreement 
with  data  obtained  from  a  model  of  the  present  design.  The 
only  important  source  of  error  originates  from  the  viscous 
losses  in  the  sink  region,  which  are  not  fully  analyzed.  How 
ever,  it  is  shown  that  this  error  does  correlate  with  the 
governing  parameters,  so  that  the  performance  can  be  pre¬ 
dicted  to  within  about  10  percent  once  a  single  experiment 
to  measure  the  losses  in  the  sink  is  conducted  on  a  given 
pickoff-drain  design.  This  is  verified  with  several  experi¬ 
ments  . 

The  results  predict  that  the  frequency  response  of 
the  sensor  can  be  increased  by  about  an  order  of  magnitude 
over  that  of  the  present  designs  without  reducing  sensitivity 
Although  the  flow  field  in  the  fully  developed  range  could  be 
less  noisy  than  in  the  small  boundary  layer  range,  the  noise 
level  of  the  present  test  model  is  not  reduced  in  the  fully 
developed  range  because  of  its  design. 

Final  discussion  covers  design  changes  that  should 
significantly  reduce  the  noise  levei  and  make  the  sensor  more 
rugged.  It  also  mentions  planned  future  analytical  work 
aimed  at  developing  a  sensor  that  compensates  for  environ¬ 
mental  changes. 
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PREFACE 


This  paper  presents  some  of  the  rate  sensor  work 
done  at  HDL  since  1968 .  At  that  time  the  basic  solution 
for  flow  in  the  chamber  of  a  particular  rate  sensor  concept, 
which  was  first  tested  in  1965,  was  completed. 

Although  this  paper,  as  HDL's  rate  sensor  program, 
is  mostly  theoretical,  it  may  have  wide  practical  impact. 
New  concepts  directly  suggested  by  the  analysis  are  (1)  a 
rate  sensor  with  effectively  no  time  delay  and  a  corner 
frequency  in  excess  of  1  kHz,  (2)  a  rate  sensor  with  much 
less  signal  drift  than  the  present  model,  (3)  a>  sensor 
whose  rate  sensitivity  is  immune  to  environmental  factors, 
and  (4)  a  single  sensor  that  measures  both  angular  rate, 
and  angular  acceleration. 

These  would  seem  to  have  both  commercial  and  mil¬ 
itary  applications;  however,  corroborating  experimental 
data  is  largely  absent. 

This  early  disclosure  is  being  made  should  some 
other  laboratory  be  in  a  better  and  more  urgent  position 
to  take  advantage  of  these  findings  and  pursue  the  neces¬ 
sary  exploratory  and  developmental  efforts. 
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INTRODUCTION 


Background 


There  has  been  much  interest  in  the  flueric  device 
known  as  the  vortex  rate  sensor  since  its  invention  about  a 
decade  ago.  In  this  device,  fluid  is  forced  radially  into  a 
short  cylindrical  chamber  through  the  porous  outer  wall  (the 
coupler)  and  exhausts  through  a  drain  in  one  or  both  end 
walls  at  the  chamber  axis  (fig.  1) .  This  sensor  produces  a 
signal  that  is  a  function  of  the  rate  at  which  it  turns  in 
inertial  space  about  its  longitudinal  axis.  In  most  designs 
the  signal  is  obtained  in  the  drain  from  an  air  foil  pickoff 
whose  output  is  a  function  of  the  streamline  angle  generated 
by  the  rotation. 

There  are  three  quantities  of  principal  interest  in  a 
rate  sensor  output:  (1)  sensitivity,  or  signal  gain;  (2)  rate 
threshold,  i.e.,  minimum  angular  rate  that  can  be  resolved 
above  noise;  and  (3)  the  frequency  response.  Signal  gain  has 
been  of  primary  interest  in  the  previous  literature  (ref  1-5)*; 
noise  (ref  6)  and  sensor  dynamics  (ref  7)  have  undergone  less 
extensive  investigation.  None  of  these  studies  has  examined 
thoroughly  the  analytical  relationships  among  these  quantities. 

For  many  applications  in  vehicle  guidance  or  stabiliza¬ 
tion,  the  present  vortex  rate  sensor  has  two  shortcomings,  a 
high  noise  level  and  an  unsatisfactory  frequency  response. 

An  effort  is  made  here  to  determine  whether  these  limitations 
on  noise  level  and  frequency  response  are  inescapable. 


-DRAIN 


CHAMBER 


o  o' 
< 

rwi 

%.* 

.o  * 

4 

*  . o  ( 

o-°‘ 

| 

0C*O 

•  °C 

I  0  * 

20' 

j 

3.  0 

•o*; 

•o- 

La 

-COUPLER 


PRESSURE 

PICKOFF 


Figure  1.  Vortex  rate  sensor. 


See  pp.  72-73  for  references. 
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Nearly  all  previous  analytical  studies  of  the  rate 
sensor  have  been  concerned  with  its  operation  when  viscous 
effects  can  be  treated  from  boundary  layer  concepts.  These 
presentations  imply  that  viscous  dissipation  in  the  chamber 
should  be  small  since  it  reduces  the  vortex  strength  and 
thus  reduces  signal  gain.  However,  the  effect  of  rela¬ 
tive  boundary  layer  size  and  viscous  dissipation  on  the  fre¬ 
quency  response  and  the  noise  level  should  be  considered  before 
such  a  restriction  is  placed  on  rate  sensor  design.  The 
following  examination  of  the  rate  sensor  operation  when  the 
boundary  layers  are  small  compared  with  the  chamber  height  in¬ 
dicates  that  both  the  sensor  dynamics  and  signal-to-noise  (S/N) 
ratio  should  be  relatively  poor. 


1.2  Dynamic  Characteristics  when  the  Boundary  Layers  are 

Small 

The  sensor  dynamics  can  be  characterized  by  a  simple 
time  delay  between  a  rate  input  and  the  sensor  output  if  the 
flow  in  the  chamber  is  essentially  inviscid  (that  is,  if  the 
boundary  layers  are  small  compared  with  the  chamber  height) , 
if  the  radial  velocity  profile  at  the  coupler  is  uniform  and 
axisymmetric ,  and  if  the  transport  times  along  the  various 
streamlines  in  the  vicinity  of  the  drain  entrance  are  small 
compared  with  those  in  the  outer  portion  of  the  chamber.  Since 
the  flow  can  be  considered  incompressible  except  possibly  near 
the  drain,  the  delay  time  T  is  given  approximately  by  a  quo¬ 
tient  of  the  chamber  volume  and  Qq,  the  volumetric  flow  rate 
in  the  outer  portion  of  the  chamber: 


t  =  (i) 

Qo 

where  Rq  =  chamber  radius,  and  h  =  chamber  height. 

If  no  dissipation  occurs  in  the  chamber,  the  velocity 
components  in  the  drain  are  independent  of  h  if  the  flow  rate 
is  fixed.  This  implies  that  the  signal  gain  of  a  pickoff  in 
the  drain  is  independent  of  h. 

If  h  is  made  small,  T  decreases  and  gain  does  not 
change  if  the  boundary  layers  are  small.  However,  for  fixed 
Q0 ,  v  and  RQ,  the  boundary  layer  size  will  become  significant 
if  h  is  less  than  some  as  yet  unspecified  value.  In  that 
case,  (1)  will  no  longer  hold,  and  the  gain  will  depend  on  h. 


This  simple  approach  indicates  that  the  height  should 
be  made  so  small  that  the  relative  size  of  the  boundary  layers 
becomes  significant,  because  this  results  in  greater  frequency 


response  at  no  sacrifice  in  gain  if  Q  is  fixed.  A  further 
decrease  in  height  should  cause  a  further  improvement  in  fre¬ 
quency  response  ibut  result  in  a  sacrifice  of  gain)  since 
the  flow  would  become  strongly  coupled  to  the  top  and  bottom 
plates  and  would  therefore  receive  angular  rate  information 
that  is  an  update  of  the  information  it  received  at  the  coupler. 


1.3 


Resolution  of  tile  Present  Sensor 


1.3.1  Streamline  Angle  Resolution 


The  problems  concerning  the  S/N  ratio  of  present  sen¬ 
sors  can  best  be  seen  by  using  an  inviscid  analysis  to  obtain 
an  approximate  value  for  the  angle  between  a  fluid  streamline 
in  the  chamber  and  a  radial  line.  Empirical  values  can  then 
be  used  in  this  relationship  to  indicate  the  order  of  magnitude 
of  the  disturbance  that  limits  the  resolution  of  this  angle 
in  the  present  sensor. 


I 


When  the  sensor  is  turning  at  a  constant  rate  S20,the 
angle  a  that  a  streamline  makes  with  a  radial  line  is  given  by 


a~  orctan 


|Vr| 


(?) 


where  Vg  =  tangential  velocity  at  a  radial  distance  r  from 
the  longitudinal  axis  of  the  chamber,  and  Vr  =  radial  velocity 
at  r.  If  Vg  is  small,  then 


a  =  ^_ 
“  Ivrl 


(3) 


If  the  sensor  is  turning  at  a  constant  rate 
and  if  the  flow  is  incompressible  and  inviscid,  then 


vs 


,  I20R02 


(4) 


and 


r  r  2Th  r 


where  Uq  =  average  radial  velocity  at  RQ. 


(5) 


3 

5 


Sr- 


Therefore, 


|a|Jft^o=2|flc|T  (6) 


A  value  for  the  minimum  change  in  angle  of  the  vortex 
streamlines  in  the  chamber  that  can  be  resolved  by  the  pressure 
pickoff  can  be  obtained  by  substituting  in  (6)  the  empirical 
values  claimed  for  the  rate  threshold  and  response  time  of 
typical  vortex  rate  sensors.  Values  obtained  in  this  manner 
are  in  the  range  from  0.01  to  0.001  deg. 


1.3.2  Fluid  Phenomena  Which  Cause  Output  Noise 

There  are  three  factors  that  inhibit  resolution  of 
the  streamline  angle  in  present  rate  sensors:  high  frequency 
noise,  drift  (low  frequency  noise) ,  and  null  shift  (a  dc 
shift  of  the  entire  curve  of  output  signal  versus  angxe  rate 
during  operation,  or  a  lack  of  repeatability  of  the  curve  when 
the  power  supply  is  shut  off  and  then  turned  on  again)  .  Al¬ 
though  the  exact  origin  of  these  disturbances  is  not  fully 
understood,  an  intuitive  examination  can  shed  some  light  on 
the  nature  of  the  problem. 

Most  of  the  high  frequency  noise  probably  results 
from  turbulence  in  the  drain  region.  If  the  gain  of  the  de¬ 
vice  is  to  be  high,  the  mass  flow  rate  must  be  large  and  the 
drain  radius  small.  Since  the  Reynolds  number  in  the  drain  is 
directly  proportional  to  the  mass  flow  rate  and  inversely  pro¬ 
portional  to  the  drain  radius,  it  is  invariably  true  that  in 
a  high-gain  rate  sensor  of  the  present  design  the  Reynolds 
number  in  the  drain  is  higher  than  the  critical  value  for  the 
onset  of  turbulence. 

A  portion  of  the  high  frequency  noise  is  probably 
due  to  turbulence  generated  when  the  probe  is  inserted  into 
the  flow.  Even  if  the  undisturbed  flow  were  laminar,  separa¬ 
tion  would  occur  at  some  point  downstream  of  the  sensing  ports 
on  the  surface  of  a  typical  pickoff  at  a  typical  operational 
condition.  The  propagation  of  this  noise  in  the  upstream 
direction  would  then  lead  to  noise  at  the  sensing  ports. 

The  signal  drift  probably  is  due  in  part  to  a  change 
in  the  distribution  of  the  fluid  at  the  outer  portion  of  the 
chamber.  A  change  of  environmental  conditions,  power  supply 
level,  or  downstream  load  could  cause  a  change  in  the  degree 
of  distortion  of  the  field  caused  by  the  manifolding  system 
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or  the  coupler,  a  ready  source  of  anomalies.  Because  the 
pickoff  in  the  present  design  extends  across  the  origin  in 
the  drain,  the  output  of  the  sensor  is  very  sensitive  to  the 
degree  of  nonaxisymmetry  in  the  field. 


The  signal  drift  probably  is  also  related  in  another 
fashion  to  variations  in  power  supply  level.  The  pickoff  ir. 
the  typical  sensor  is  of  the  push-pull  variety,  one  that  has 
sensing  ports  on  either  side  of  the  stagnation  point.  A 
change  in  angle  of  attack  thus  generates  pressure  changes  of 
opposite  signs  at  the  two  ports,  and  a  differential  signal  is 
then  a  measure  of  the  angle,  or  a  function  of  the  angular  rate. 


It  is  not  possible  to  orient  the  pickoff  so  that  the 
differential  output  is  perfectly  zero  when  the  sensor  is  sta¬ 
tionary.  As  a  result  the  differential  signal  at  zero  angular 
rate  is  not  zero  and  in  particular . is  dependent  upon  the 
flowrate.  Thus,  variations  in  flowrate  inhibit  the  measure¬ 
ment  of  small  changes  of  angular  rate  about  zero  rate. 


The  null  shift  probably  results  from  flexing  or 
shifting  in  the  pickoff  drain  assembly.  Present  rate  sensors 
typically  have  a  very  small  pickoff  mechanically  mounted  in¬ 
side  a  small  drain,  an  assembly  which  can  easily  be  susceptible 
to  mechanical  or  thermal  disturbances.  Since  the  pickoff  is 
situated  at  a  position  where  the  magnitude  of  the  tangential 
velocity  varies  sharply  with  radial  location,  even  a  very  small 
movement  can  result  in  a  significant  output  shift. 


It  may  also  be  true  that  a  nonaxisymmetric  flow  dis¬ 
tribution  at  the  coupler  that  varies  with  supply  pressure  can 
lead  to  a  null  shift,  for  under  those  circumstances  the  features 
of  the  flow  field  at  zero  rate  could  depend  upon  the  manner  in 
which  the  power  supply  is  brought  up  to  working  level.  However, 
this  possible  noise  source  then  becomes  closely  linked  with 
one  of  the  speculated  sources  of  low  frequency  noise,  and  thus 
would  be  attacked  in  any  pregram  planned  to  eliminate  the 
latter. 


1.3.3  Outout  Noise 


The  problem  of  noise  in  the  output  signal  can  now  be 
put  into  perspective.  It  is  probably  surprising  to  most  ob¬ 
servers  that  an  angle  of  0.001  deg  can  be  resolved  in  a  turbu¬ 
lent  field  by  a  pickoff  that  may  be  sensitive  even  to  the 
thermal  stresses  induced  by  a  change  in  temperature  of  the 
working  fluid.  But  the  purported  causes  of  low  frequency  noise 
are  just  as  disquieting  when  referenced  to  this  small  angle. 


Suppose,  for  instance,  that  the  pickoff  were  mechan¬ 
ically  oriented  in  the  field  at  an  angle  of  0.1  deg  from  the 
angle  for  perfect  null.  Then  (3)  indicates  that  a  one  percent 


change  In  flowrate  would  be  sufficient  to  cause  an  output 
shift  equivalent  to  an  angular  change  of  0.001  deg.  Since 
both  of  these  requirements  are  realistic  if  not  optimistic 
estimates  of  the  accuracies  of  present  techniques,  it  would 
seem  that  this  effect  alone  could  be  sufficient  to  limit  the 
performance  of  the  rate  sensor.  However,  the  requirements  on 
axisymmetry  of  the  field  appear  at  first  glance  to  be  even 
more  stringent. 

The  continuity  equation  can  be  used  to  show  that  a 
nonaxisymmetric  upstream  disturbance  generates  perturbations 
with  respect  to  6  of  approximately  the  same  order  in  Vr  and 
Vg .  This  implies  that  the  present  noise  level  is  equivalent 
to  a  time-varying  perturbation  with  respect  to  0  in  Vr  at 
the  pickoff  that  is  about  four  orders  of  magnitude  smaller 
than  the  mean  value,  because 
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This  is  an  amazingly  low  disturbance  if  one  observes  that  little 
attention  has  been  given  to  the  designing  of  a  manifolding 
system  that  is  capable  of  supplying  the  outer  portion  of 
the  chamber  with  a  flow  whose  time-varying  distortion  with  re¬ 
spect  to  0  is  less  than  even  one  percent  of  the  mean  when  the 
sensor  is  powered  by  a  typical  power  supply.  This  implies 
that  the  noise  at  the  pickoff  due  to  entrance  effects  is  pro¬ 
bably  at  least  two  orders  of  magnitude  l^ss  than  that  at  the 
coupler . 


The  sensor  resolution  has  been  enhanced  somewhat  by 
the  design  of  the  push-pull  pickoff  discussed  previously, 
since  noise  at  the  ports  that  has  positive  correlation  tends 
to  cancel.  However,  the  only  fluid  mechanism  used  effectively 
to  condition  the  field  properly  for  precision  measurement  is 
the  fields'  natural  acceleration  as  it  flows  toward  the  sink, 
which  serves  to  decrease  upstream  nonaxisymmetry . 

It  should  not  be  surprising  that  the  present  vortex 
rate  sensor  does  not  have  the  desired  resolution  for  applica¬ 
tions  that  are  within  the  capabilities  of  other  gyro  technolo¬ 
gies  since  the  specifications  for  those  applications  would 
require  extremely  high  integrity  in  a  vortex  field.  In  an 
absolute  sense,  the  present  sensor  is  quite  effective  at 
measuring  very  small  streamline  angles,  due  at  least  as  much 
to  a  fortunate  choice  of  initial  sensor  concept  as  to  any 
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later  design  effort.  The  sink  field,  chosen  because  it  pro¬ 
vides  a  mechanism  for  amplifying  small  signals  at  the  coupler, 
has  a  side  benefit  in  its  effect  on  noise.  The  only  reason 
that  the  present  vortex  rate  sensor  works  as  even  a  moderately 
effective  sensor  seems  to  be  that  the  acceleration  of  the  sink 
field  leads  to  a  drastic  reduction  of  the  noise  due  to  entrance 
conditions . 

It  seems  clear  that  a  limit  has  been  reached  in  sensor 
performance,  and  that  it  would  be  unreasonable  to  expect  any 
improvement  unless  a  great  deal  of  care  is  taken  to  eliminate 
noise.  Of  the  three  types  of  noise  discussed,  signal  drift 
seems  to  be  the  most  troublesome.  Null  shift  can  be  elimina¬ 
ted  for  some  applications  by  using  a  high-pass  filter,  and 
high  frequency  noise  by  using  a  low-pass  filter.  However,  the 
low  frequency  noise  appears  in  any  bandwidth  chosen  for  the 
sensor  response.  If  this  latter  noise  is  related  in  part,  as 
suggested,  to  entrance  conditions,  then  special  care  should  be 
taken  to  properly  condition  the  flow  field  for  precision 
measurement . 


1.4  Proposed  Changes 

It  was  suggested  in  Section  1.2  that  h  should  at 
least  be  made  so  small  that  the  boundary  layer  size  is  not 
negligible  when  compared  to  h.  This  change  would  seem  to  be 
compatible  with  the  discussion  in  Section  1.3  concerning  the 
]ow  frequency  noise,  because  the  most  logical  methods  of  con¬ 
ditioning  the  flow  would  lead  to  significant  viscous  effects 
in  the  chamber. 

The  most  obvious  way  to  condition  the  flow  would  be 
to  make  the  chamber  radius  so  large  that  the  entrance  conditions 
dc  not  noticeably  influence  the  flow  near  the  pickoff.  This 
requires  at  the  same  time  that  the  coupler  have  no  effect  on 
the  signal  output.  The  signal  would  then  be  introduced  by 
means  of  wall  shear,  a  much  less  noisy  means  than  the  use  of  a 
coupler,  which  has  been  shown  (ref  8)  to  introduce  very  high 
distortion  in  the  flow  field.  If  the  chamber  is  sufficiently 
large,  the  only  output  noise  seen  by  a  rugged  pickoff  would 
be  that  due  to  changes  in  the  power  level  and  that  which  propa¬ 
gates  upstream  from  the  drain.  But  since  the  signal  propagates 
in  a  direction  normal  to  the  wall  or  in  a  direction  normal  to  the 
direction  of  propagation  of  these  noises,  it  may  be  possible 
by  proper  signal  processing  to  significantly  decrease  the  ef¬ 
fect  on  the  output  of  even  that  fluid  noise  which  careful  design 
does  not  eliminate. 

No  mention  has  been  made  of  the  effect  that  increasing 
the  chamber  radius  would  have  on  frequency  response  and  gain 
of  the  sensor.  It  might  help  at  this  point  to  try  to  establish 
intuitively  what  this  effect  would  be.  However,  it  is  doubtful 
that  such  a  discussion  would  encourage  the  reader  to  endure  the 
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lengthy  analysis  of  the  problem  to  follow  later  since  it 
would  be  necessary  for  an  adequate  discussion  to  redefine  old 
concepts  for  use  in  a  new  area  with  definitions  that  antici¬ 
pate  later  mathematical  difficulties ,  and  thus  possibly  only 
generate  confusion. 


Perhaps  some  reader  concern  can  be  allayed  if  it  is 
pointed  out  that  a  theory  is  developed  here  which  adequately 
predicts/  among  other  things,  the  gain  of  the  device  for 
all  values  of  h  for  which  the  boundary  layers  merge.  If 
the  pickoff  is  located  in  the  drain,  then  the  gain  decreases 
monotonically  to  zero  as  the  height  is  decreased.  For  suf¬ 
ficiently  large  h,  it  is  found  that  the  curve  of  gain  versus 
height  is  nearly  flat,  that  is,  that  gain  is  relatively  in¬ 
dependent  of  h.  At  small  h,  it  is  a  linear  function  of  it. 
In  this  later  region,  the  theory  indicates  that  the  coupler 
ha  no  influence  on  the  output  signal.  The  gain  is  reduced 
by  a  factor  of  about  two  in  reducing  h  from  a  value  in 
the  "flat"  range  to  one  in  the  upper  portion  of  this  linear 
rar.ge,  but  the  frequency  response  is  increased  by  a  factor 
of  about  twenty. 


It  can  be  then  shown  that  if  one  were  to  increase 
Rq  and  decrease  h  in  such  a  fashion  that  some  characteristic 
frequency  (such  as  that  at  which  the  phase  lag  is  90  deg)  re¬ 
mains  constant,  the  gain  would  increase.  It  would  increase 
by  about  an  order  of  magnitude  while  the  configuration  is 
changed  from  one  in  which  the  flow  in  the  chamber  is  nearly 
inviscid  to  one  in  which  the  entrance  conditions  do  not  appre¬ 
ciably  affect  the  signal. 


If  the  gain  and  frequency  response  are  independent  of 


Rn,  as  they  would  be  under  this  latter  condition,  the  coupler 


could  then  be  made  large  enough  to  properly  condition  the  flow. 
The  limit  on  size  would  be  fixed  by  physical  restrictions.  If 
these  do  not  allow  a  sufficiently  large  R0  for  the  desired 
noise  reduction,  then  several  schemes  other  than  a  single  cylin¬ 
drical  cavity  could  be  tried.  One  would  be  a  stack  of  thin 
annular  chamber*  mounted  coaxially  above  the  rate  sensor  cham¬ 
ber.  By  proper  interconnections,  the  fluid  could  alternately 
flow  inward  and  outward  through  successive  chambers  as  it  pro¬ 
gressed  through  the  stack  before  it  enters  the  working  chamber 
and  flows  out  the  sink.  If  the  deceleration  the  fluid  exper¬ 
iences  in  each  alternate  outward  stage  does  not  permit  stable 
flow,  either  the  chamber  height  could  be  changed  with  radius 
to  ensure  stability  or  other  schemes  should  be  tried.  It  is 
not  important  to  discuss  various  schemes  at  this  time,  partic¬ 
ularly  since  it  appears  that  their  number  is  limited  only  by 
the  imagination  of  the  designer  for  a  particular  set  of  physical 
restrictions. 
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In  summary,  it  seems  evident  that  the  operation  of  a 
rate  sensor  should  be  considered  for  the  general  regime  in 
which  the  viscous  effects  are  not  negligible,  since  the  changes 
required  to  increase  frequency  response  and  reduce  noise  lead 
either  directly  or  indirectly  to  the  consideration  of  signifi¬ 
cant  viscous  effects. 


Boundary  Layer  Effects 


1.5.1  Condition  for  Fully  Developed  Flow 


In  order  to  discuss  the  effect  of  viscous  losses  on 
the  sensor  characteristics,  it  is  necessary  first  to  establish 
which  dimensionless  viscous  parameter  specifies  the  relative 
size  of  the  boundary  layers.  Approximate  values  can  then  be 
sought  for  this  parameter  that  divide  its  range  into  regions 
of  negligible,  partial,  and  complete  viscous  interaction  of 
the  top  and  bottom  plates . 


The  boundary  layer  problem  for  sink  flow  over  a  single 
flat  disc  has  been  attacked  by  Sarpkaya,  Goto,  and  Kirshner 
(ref  9)  through  the  use  of  an  integral-momentum  method.  The 
results  obtained  with  their  method  are  shown  in  figure  2.  The 
boundary  layer  is  thickest  at  r  =  0.55  RQ  with 
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RADIAL  DISTANCE  r/R0 

Figure  2.  Boundary  layer  problem  for  sink  flow  over  a 
single  flat  disc  (method  of  Sarpkaya,  Goto 
and  Kirshner) . 
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There  has  been  no  known  published  work  that  examines 
in  detail  the  problem  of  profile  development  for  sink  flow  be¬ 
tween  two  flat  discs.  However,  an  approximate  value  for  the 
pertinent  viscous  parameter  at  which  the  radial  velocity  pro¬ 
file  first  becomes  fully  developed  can  be  obtained  by  usinc 
(8)  together  with  some  physical  arguments  based  on  the  profile 
development  for  plane  flow  between  two  flat  plates.  The 
phrase  "fully  developed"  is  used  hare  to  describe  a  profile  at 
a  radius  where  the  influence  of  the  entrance  conditions  is  no 
longer  significant,  even  though  the  velocity  profile  might 
change  with  decreasing  radius  due  to  acceleration  effects. 

The  laminar  boundary  layer  thickness  for  plane  flow 
over  a  flat  plate  at  zero  angle  of  incidence  to  a  uniform  ex¬ 
ternal  field  is  given  (ref  10)  by 
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where  Xi  is  the  distance  from  the  leading  edge  and  U  is  the 
free  stream  velocity.  The  entrance  length  £e  /  for  fully 
developed  plane  flow  between  parallel  flat  plates  separated 
by  a  distance  h  is  given  (ref  11)  by 


£e  =  0.04^ 


(10) 


Equations  (9)  and  (10)  can  be  used  to  establish  an 
approximate  criterion  in  terms  of  the  boundary  layer  thickness 
over  a  single  plate  at  xi  =  ie  and  the  plate  separation  dis¬ 
tance  h.  This  approximate  criterion  is  that  the  boundary 
layer  on  a  single  plate  (in  the  absence  of  the  other  plate) 
has  a  thickness  at  Ze  of 


UOh 


(11) 


Now  consider  sink  flow  between  two  flat  discs  as  h  is 
decreased  from  a  large  value  while  Q^,  v  and  RQ  are  fixed.  At 
some  critical  value  of  h,  the  flow  will  become  fully  developed 
at  some  radius  in  the  chamber.  If  this  radius  is  just  that  at 
which  the  boundary  layer  over  a  single  disc  is  at  its  thickest 
and  if  a  relationship  similar  to  (11)  holds,  then  an  approxi¬ 
mate  value  for  the  modified  Reynolds  number,  Re^,  for  this  case 
can  be  found  from  (8)  given  by 
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It  should  be  noted  that  (9)  refers  to  a  boundary  layer 
in  the  absence  of  an  external  pressure  gradient  and  that  an 
accelerating  radial  pressure  gradient  is  present  in  sink  flow. 
Since  the  boundary  layer  profile  in  an  accelerating  pressure 
gradient  is  blunter  in  the  region  near  the  free  stream,  it 
could  be  argued  that  the  boundary  layer  thickness  from  criter¬ 
ion  (11)  is  not  suf f iciently  large  to  ensure  profile  develop¬ 
ment  between  parallel  discs.  However,  the  velocity  profile 
that  the  fluid  attempts  to  attain  is  made  blunter  by  the  accel¬ 
eration  effects  than  the  profile  that  would  exist  in  the  absence 
of  an  accelerating  pressure  gradient,  so  the  degree  by  which 
the  viscous  forces  must  alter  a  uniform  entrance  profile  is 
not  so  severe  for  flow  in  the  presence  of  an  accelerating  pres¬ 
sure  gradient.  The  sign  of  the  error  in  the  approximation 
used  in  obtaining  (11)  must  remain  unspecified  at  this  time; 
however,  the  magnitude  of  the  error  should  not  be  large. 


1.5.2  Features  of  Flow  Field  for  Re^  <  Re^ 

crit 

If  the  free  stream  thickness  is  small,  the  additional 
acceleration  effects  caused  by  the  displacement  of  the  stream¬ 
lines  toward  the  centerplane  should  cause  the  boundary  layers 
to  become  relatively  flat  midway  in  the  chamber.  Consequently, 
if  the  point  in  the  chamber  at  which  the  flow  may  be  considered 
fully  developed  is  much  different  from  R0  at  Re^  ,  the  point 

crit 

rapidly  moves  toward  the  coupler  (R0)  as  Re^  is  decreased  below 
the  critical  value.  Then  the  flow  can  be  considered  fully  de¬ 
veloped  throughout  most  of  the  chamber  for  values  of  Re^  less 
than  the  critical  value.  ' 


1.5.3 


Features  of  Flow  Field  for  Re^ 
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Now  consider  the  flow  when  Re^  is  large.  The  shape 
of  the  velocity  profile  in  the  boundary  layers  is  specified 
by  the  parameter  A  where 

.  .  s£  t<(Vr  free  stream)  ( 
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If  the  boundary  layer  effects  on  the  free  stream  are  small  and 
the  flow  is  incompressible,  (5)  can  be  used  for  the  free  stream 
velocity  in  (13)  to  show  that 
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At  a  given  r/RQ,  the  quantity^L-  Re^jis  a  constant  to  the  ap 
proximation  of  small  boundary  layers,  since  figure  2  indicates 


is  a  function  of  ^ —  only.  Then  A  is  a  function  of 


that/— 2 
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•£—  only,  and,  in  particular,  is  independent  of  the  relative 

size  of  the  boundary  layer.  Thus,  the  boundary  layer  features 
can  be  treated  simultaneously  for  all  Re'  when  the  relative 
size  is  small.  Such  a  discussion  could  Be  extended  to  a  good 
approximation  into  the  range  of  Re^  in  which  the  boundary 
layers  interact  simply  by  observing  what  effect  the  asso¬ 
ciated  change  in  A  would  have. 


From  figure  2, 


Alr=RmaxS22  (15> 

where  Rmax  is  the  radius  at  which  the  boundary  layer  thickness 
is  a  maximum  (0.55Rq) .  The  value  of  A  at  Rmax  increases  as 
the  thickness  of  the  boundary  layers  approaches  the  chamber 
height,  due  to  the  added  acceleration  of  the  free  stream. 


When  the  boundary  layers  first  touch,  the  interactions 
near  the  center  plane  are  extremely  small  due  to  the  bluntness 
of  the  velocity  profile  in  this  region.  The  most  severe  inter¬ 
actions  will  occur  near  For  the  value  of  A  in  (15)  ,the 

velocity  profile  for  a  single  boundary  layer  is  so  blunt  (ref  12) 
that  the  boundary  layers  on  the  top  and  bottom  plates  could 
overlap  by  as  much  as  30  percent  without  causing  more  than  a 
few  percent  variation  in  the  velocity  throughout  the  overlapping 
region.  When  the  overlap  is  greater,  the  profile  obtained  near 


the  centerplane  by  a  simple  superposition  of  the  boundary 
layer  profiles  does  not  differ  significantly  from  a  good  es¬ 
timate  of  the  actual  solution  in  that  region.  Since  the 
greatest  signal  contribution  of  the  pickoff  comes  from  the 
fluid  passing  through  the  region  near  the  centerplane  of  the. 
chamber,  the  output  or  tne  sensor  can  be  described  fairly  well 
from  boundary  layer  analysis  even  if  the  overlapping  of  the 
boundary  layers  is  significant. 

The  increase  in  A  that  occurs  as  the  relative  size 
of  the  boundary  layers  grows  causes  the  viscous  effects  near 
the  centerplane  to  be  weaker  than  the  previous  discussion  an¬ 
ticipates.  Since  it  is  more  valid  to  approximate  viscous 
effects  by  a  superposition  of  boundary  layers  if  these  effects 
are  weaker,  the  previous  discussion  still  holds. 

It  would  seem,  then,  that  the  interactions  of  the 
boundary  layers  are  not  significant  in  the  case  of  slight  over¬ 
lapping.  Furthermore,  even  in  the  case  of  severe  overlapping, 
the  sensor  output  can  be  described  adequately  through  boun¬ 
dary  layer  analysis. 


1.5.4 


General  Features  of  the  Flow  Field 


In  summary,  the  range  of  Re^  for  rate  sensors  having 
a  cylindrical  chamber  can  be  divided  rather  sharply  into  two 
regions,  one  in  which  boundary  layer  analysis  applies  and  one 
in  which  the  radia.l  velocity  profile  can  be  considered  fu3  ly 
developed  throughout  the  entire  chamber.  The  value  of  Re^ 
which  divides  these  two  regions  is  about  7. 


Purpose  of  Investigation 


It  was  indicated  in  Section  1.4  that  the  operation  of 
the  sensor  should  be  considered  when  the  relative  size  of  the 
boundary  layers  is  not  negligible.  Section  1.5.4  then  indicates 
that  the  characteristics  of  the  sensor  should  be  examined  in 
the  fully  developed  range.  It  cannot  be  established  a  priori 
whether  some  gain  needs  to  be  sacrificed  to  both  increase  the 
frequency  response  and  improve  the  rate  threshold  to  the  de¬ 
sired  level;  however,  that  is  not  of  great  consequence  at  the 
present,  since  the  shortcomings  cf  the  present  vortex  rate 
sensors  are  its  poor  rate  threshold  and  poor  frequency  response, 
not  its  gain  characteristics.  Since  the  sensor  feature  for 
which  the  rest  of  the  control  system  usually  can  most  easily 
compensate  is  a  modest  deficiency  of  gain,  a  limited  sacrifice 
in  gain  is  not  an  obstacle. 
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An  investigation  of  the  characteristics  of  a  vortex 
rate  sensor  for  that  range  in  which  the  entrance  conditions 
do  not  appreciably  influence  the  radial  velocity  is  then  of 
paramount  interest. 


2.  ANALYSIS 

2.1  Formulation  of  ProbLem 

2.1.1  Analytical  Model 

In  what  follows  an  approximate  solution  will  be 
found  for  the  governing  equations  for  inward  flow  in  a  thin 
cylindrical  chamber  when  the  angular  rate  has  been  a  sinu¬ 
soidal  function  of  time  for  an  infinite  length  of  time.  The 
no-slip  condition  will  be  imposed  on  both  the  radial  and  tan¬ 
gential  velocities  at  the  top  and  bottom  plates  and  on  the 
tangential  velocity  at  the  coupler  wall.  The  entrance  effects 
on  the  radial  profile  will  be  neglected,  in  keeping  with  the 
purpose  of  this  study  as  stated  in  Section  1.6.  The  influence 
of  the  drain  on  the  flow  in  the  chamber  will  be  ignored;  the 
reasons  for  this  restriction  will  be  discussed  briefly. 


In  attempting  to  attack  the  problem  analytically,  it 
is  best  to  divide  the  physical  problem  into  three  regions; 
the  chamber,  the  drain  and  the  corner  region  between  the  two. 
The  chamber  is  that  region  in  which  the  velocity  normal  to  the 
top  and  bottom  plates  is  very  small.  The  drain  is  that  region 
in  which  the  streamlines  are  nearly  parallel  to  the  drain 
walls.  Since  a  different  set  of  terms  in  the  governing  equa¬ 
tions  can  be  ignored  in  each  region,  it  is  logical  to  obtain 
the  general  solutions  for  each  region  separately  and  then  to 
match  the  solutions  at  the  dividing  surfaces. 


Only  the  chamber  problem  is  discussed  in  detail  in 
this  paper,  although  the  qualitative  effect  of  the  corner  and 
the  drain  on  the  rate  sensor  characteristics  is  attacked  analy¬ 
tically.  There  are  two  reasons  for  considering  only  the  chamber 
in  detail.  First,  since  the  corner  problem  does  not  lend  itself 
readily  to  analytical  methods,  the  importance  of  the  corner 
effects  should  be  established  before  one  undertakes  the  pro¬ 
digious  task  of  solving  the  differential  equations  that  apply  in 
that  region  by  computer  techniques.  This  importance  can  be 
established  if  the  corner  is  the  only  significant  source  of 
error  that  is  ignored  in  attempting  to  describe  the  output 
characteristics  of  the  sensor,  as  will  be  done  here.  Second, 
if  the  corner  region  represents  a  noise  source  or  if  the 
viscous  losses  are  significant  in  that  reqion,  the  rate 
sensor  probably  should  be  designed  with  a  pickoff  just  up¬ 
stream  of  the  corner.  Intuition  indicates  that  the  corner 
does  have  just  such  features,  so  that  the  sensor 
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which  finally  results  front  this  effort  can  probably  be  de¬ 
signed  on  the  basis  of  the  chamber  solution  alone. 


In  keeping  with  this  restriction  on  the  scope  of  the 
problem  to  be  attacked,  it  is  assumed  (1)  that  the  drain  is 
a  small  tubular  sink  concentric  with  the  chamber  axis,  (2) 
that  the  pickoff  extends  the  length  of  the  chamber  near 
the  sink,  (3)  that  the  pickoff  responds  as  an  angle  of  at¬ 
tack  sensor,  and  (4)  that  the  pickoff' s  pressure  output 
is  an  average  over  the  entire  length  of  the  pressure  at 
the  angular  location  of  the  sensing  port.  Figure  3  is  an 
illustration  of  this  model. 

R=0 


-Z=0 


ANGLE  OF 

ATTACK 

SENSOR 


TUBULAR 

SINK 


A-A 

(ENLARGED) 


Figure  3.  Vortex  rate  sensor  analytical  model. 


The  radial  location  of  the  pickoff  in  the  mathematical 
model  is  set  equal  to  the  drain  radius  Rj  of  the  experimental 
model.  The  output  of  the  pickoff  located  in  the  drain  of  the 
real  sensor  is  assumed  to  be  a  linear  function  of  the  output 
of  this  hypothetical  pickoff  in  the  mathematical  model  of  the 
chamber.  The  implications  of  this  assumption  will  be  explored 
as  the  analysis  progresses. 


The  radius  of  the  sink  tube  should  be  small  enough 
that  the  shear  generated  on  its  surface  does  not  affect  the 
flow  near  the  pickoff  because  the  corner-drain  effects  are  to 
be  ignored.  It  will  be  shown  later  that,  for  the  values  of 
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flow  paramaters  in  this  experiment,  the  sink  tube  radius  need 
only  be  very  slightly  smaller  than  Rj  .  This  requires  that 
in  the  mathematical  model  the  pressure  at  the  sink  tube  wall 
must  be  slightly  lower  than  that  in  the  larger  drain  of  the 
real  model  so  that  the  flowrate  through  both  models  will  be 
the  same  for  the  same  supply  pressures. 

The  statement  that  the  drain  tube  need  only  be  slightly 
smaller  than  R,  may  be  surprising  if  one  expects  the  center 
of  the  field  to  behave  in  rigid  body  fashion,  which  is  a  fea¬ 
ture  of  three-dimensional  sink  flow.  It  can  be  shown  (App  A) 
that  no  such  rigid-body  core  exists  for  steady,  inward  flow 
between  two  concentric,  porous,  rotating  cylinders  that  are 
very  long  in  the  axial  direction.  If  the  cylindrical  cavity 
is  short,  the  viscous  effects  from  the  top  and  bottom  plates 
complicate  the  mathematics  and  the  physical  understanding  of 
the  problem,  but  the  essential  features  of  the  flow  fields  are 
quite  similar. 

The  drain  tube  cannot  be  made  arbitrarily  small  for 
a  given  pressure  and  flow  rate  in  the  exterior  part  of  the 
chamber  since  the  pressure  at  the  drain  wall  would  have  to  be 
zero  at  some  finite  radius  in  order  that  such  a  field  exist. 

This  is  the  so-called  "limit  circle"  effect.  But  in  reality, 
the  exterior  field  cannot  be  established  unless  the  pressure 
at  the  drain  tube  is  set  at  some  meaningful  value,  nor  does 
a  conscientious  designer  ever  specify  the  flow  rate  and  supply 
pressure  at  which  a  model  is  to  work  without  checking  the  pres¬ 
sure  drop  across  the  model  to  see  if  the  exit  pressure  is 
physically  realizable.  To  restrict  the  range  of  r  because 
of  the  artificial  concept  of  a  "limit  circle"  for  some  situa¬ 
tions  that  cannot  be  achieved  in  reality  is  to  needlessly 
complicate  the  resulting  design  equations.  With  this  under¬ 
standing  of  the  problem  the  radius  of  the  drain  tube  is  al¬ 
lowed  to  be  a  small  value,  6j . 

Figure  4  shows  a  schematic  of  the  drain  region  of  the 
rate  sensor  that  was  tested  in  this  study.  The  drain  dimen¬ 
sions  are  so  small  that  the  size  of  the  pickoff  output  port 
must  be  of  the  same  order  as  that  of  the  drain  to  make  the 
power  output  of  the  pickoff  readily  compatible  with  present 
flueric  amplifiers.  In  addition,  since  the  Reynolds  number 
based  on  pickoff  tube  diameter  and  velocity  in  the  drain  tube 
is  of  the  order  of  104  in  this  experiment,  the  local  distur¬ 
bance  due  to  the  picxoff  probably  causes  any  variation  in 
signal  level  across  the  stream  to  be  reduced,  if  not  washed 
out.  Because  of  this,  one  of  the  basic  assumptions  of  this 
analysis  is  that  the  pickoff  integrates  across  the  entire 
flow  stream,  as  indicated  previously  in  this  section. 
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Figure  4.  Rate  sensor  drain  region. 


2.1.2  Assumptions 

The  following  assumptions  are  made  about  the  flow  in 
the  chamber: 

(a  )  The  flow  is  laminar  and  incompressible. 

(b  )  The  flow  is  axisymmetric  with  respect  to  the  longi¬ 

tudinal  axis  of  the  chamber. 

(c  )  The  vortex  is  weak,  (that  is  I  V#!  «  I  Vrl  ) . 

(d  )  The  displacement  of  a  fluid  particle  in  a  direction 

normal  to  the  end  walls  as  it  flows  toward  the  sink 
is  small  compared  with  the  chamber  height. 

(e  )  The  body  forces  are  small. 

Assumption  (d.)  and  the  conditions  that  the  chamber 
is  thin  can  be  used  to  show  that 
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If  the  assumptions  (b) ,  (c) ,  and  (e) ,  and  (16)  are 
used,  these  equations  become 


/dv. 

dv,  \ 

dp  . 

d2v,  .  : 

v*-+ 

~3tH 

^dr2  r 

!  (dve+  avg+  V|A  +  + 

y  dt  r  dr  r  r  J  ydr^  r  dr  r2  dz% ) 


4^  +  ^  +  4^  = 


2 . 1 . 3 . 1  Simpli  :  if.::-  Radial  Momentum  Equation 

Notice  that  Vg  does  not  appear  in  (21) .  Since  the 
time-dependent  boundary  conditions  apply  only  to  Vfl ,  under  the 
above  assumptions  and  boundary  conditions  Vr  is  independent 
of  time.  Thus,  the  time  derivative  in  equation  (21)  can  Tae 
dropped.  Before  the  rest  of  the  equation  can  be  reduced,  it 
is  necessary  first  to. examine  the  relationship  between 
Vr/r  and  3Vr/3r. 


If  the  radial  velocity  profile  at  RQ  has  only  one 
maximum,  then  at  all  r  it  will  have  only  one  maximum.  This 
implies  that  Vz  can  have  only  three  zeroes,  one  at  each  wall, 
and  one  at  the  point  where  Vr  is  a  maximum.  In  addition,  if 
the  radial  profile  at  Rq  is  not  sharp,  V2  must  be  a  rather 
smooth  function  of  z  for  all  r.  These  facts,  together  with 
the  observation  that  the  magnitude  of  Vr  increases  smoothly 
as  r  decreases,  can  be  used  with  inequality  (16)  to  show  that 


dr  dz 


This  has  the  same  consequence  in  the  continuity  equation  as 
saying  that  the  mass  flux  through  an  elemental  surface  lying 
in  the  z-plane  can  be  ignored  in  comparison  to  that  uhrouah 
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an  elemental  cylindrical  surface  perpendicular  to  r. 


Equation  (24)  and  inequality  (25)  can  then  be  used 
to  show  the  approximate  relationship 

..£=-*£-***  *£  (26) 
r  dr  dz  —  dr 

This  expression  implies  that  (21)  can  be  closely  approximated 
by 


V,2  dp  d2vr 


If  one  lets 


-A(r)k  (r,i) 


and  integrates  (24)  with  respect  to  z  from,  -h/2  to  h/2,  it 
can  be  shown  that 


f h/  k(r,  z)dz 


where  k(r,o)  =  1,  k(r,  ±  h/2)  =  0. 

The  condition  that  k(r,o)  =  1  has  been  artificially 
imposed  in  an  effort  to  clarify  the  problem.  The  function  k 
represents  a  normalized  velocity  profile,  and  the  function 
A  defines  the  magnitude  of  the  velocity  in  the  centerplane. 

If  (23)  and  (28)  are  used,  (27)  becomes 


A  dZk  _  A2/>k2  _  dp 
**  r  3i2  r5  dr 


Under  assumption  (d) ,  A  is  a  weak  function  of  r,  be¬ 
cause  A  would  be  constant  if  the  displacement  were  zero.  (If 
it  is  zero,  then  the  profile  is  fixed  and  k(r,z)  =  k(z)  in  (29) 
so  that  A(r)  would  be  independent  of  r) .  Since  k  is  monotonic 
from  z  =  0  to  z  =  h/2,  the  ratio  of  the  two  terms  involving  k 


in  (30)  can  then  be  approximated  by 


a  ifiL 
*  r  di2 


where  R{r)  has  a  dependence  on  r  of  the  order  of  r2  and  the 
bar  refers  to  an  average  over  z.  Thus,  at  sufficiently  large 
r  'ajsordng  that  r  were  not  restricted  to  values  less  than 
Rq) ,  the  numerator  of  the  left  side  of  (31)  is  much  larger 
than  the  denominator,  and  (30)  becomes 

d*ka>  rdp 

-/*A«  -72  -  —  (32) 


Equation  (32)  implies  that  32k„/3z2  is  independent 
of  z.  If  this  result  is  integrated  twice  and  the  boundary 
conditions  are  imposed,  one  has 


•(•-#) 


Furthermore ,  from  ( 29 ) 


At  small  r  (except  near  the  points  r 
the  k2  term  predominates,  and  (30)  becomes 


=  0,  z  =  ±  h/2) , 


.  2 . 2_  rJdp 
A0p*o  -  ~ 


This  implies  that  k0  is  independent  of  z  for  a  point  away  from 
the  wall.  Therefore, 


for 
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II 


As  the  fluid  flows  toward  the  origin  from  a  large 
radius,  the  radial  velocity  profile  changes  smoothly  from 
parabolic  to  uniform,  and  A  decreases  to  two-thirds  of  its 
value  at  a  large  radius.  Therefore,  the  r  dependence  of 
A(r)*k  (r,z)  is  much  weaker  than  that  of  the  1/r  term  except 
possibly  near  the  points  (0,  ±  n/2) .  Thus,  the  radial  velocity 
in  the  vicinity  of  some  finite  radius  rj  can  be  approximated 
closely  by 


v. 


r 


(38) 


where  At  =  A(rt  )  and  kt  (z)  =  k(rj  ,  z) . 
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2. 1.3. 2  Simplified  Tangential  Momentum  Equation 

The  simplified  tangential  momentum  equation  is 


■5 
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dvg 

at 


+ 


(39) 


The  boundary  conditions  for  Vq  at  the  solid  boundaries  are 
given  by 

0 

Vg(  Ro.i)  =Roft0«*"» 


r 


(41) 


Vg(;,z)  = 


(42) 


-a 


where  Cj  =  sink  tube  radius. 

Although  the  effect  of  the  entrance  conditions  on  the 
radial  velocity  profile  has  been  ignored,  their  effect  on  the 
tangential  profile  is  considered,  as  defined  by  (40) .  Since 
the  type  of  equation  and  the  boundary  conditions  satisfied  by  each 
of  the  two  components  of  velocity  are  different,  the  "entrance 
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lengths"  in  general  will  not  be  the  same.  Therefore,  the 
effect  of  the  entrance  conditions  on  the  tangential  profile 
is  retained  to  discover  a  posteriori  whether  they  significantly 
influence  the  flow  field  after  entrance  effects  on  the  radial 
profile  have  been  assumed  to  disappear. 

Since  (39)  is  linear,  it  may  lend  itself  readily  to 
analysis  if  V  is  found  first.  In  this  paper,  a  method  is 
shown  for  finding  Vg  if  Vr  is  given  by  a  set  of  functions  A^ 
and  that  apply  in  appropriately  thin  annular  regions  that 
together  fill  the  chamber.  A  less  complicated  method  is 
used  to  obtain  bounding  expressions  for  the  gain  and  frequency 
response  of  the  rate  sensor. 


2.2  Solution  of  Tangential  Momentum  Equation 

2.2.1  Solution  in  a  Given  Annulus 

Assume  that  the  solution  for  Vg  in  the  region  about 
in  which  At  and  k(  apply  is  of  the  form 

+  u2(r,*)]ft0««“*  (43) 


The  function  ol  is  chosen  such  that  it  satisfies  the  boundary 
conditions  on  the  top  and  bottom  plate  for  all  r,  f>  jsuming 
that  k,  applies  over  all  r,  that  Rq  is  infinite,  c  1  that  the 
drain  effects  are  ignored.  The  function  u2  is  then  used  to 
satisfy  the  upstream  and  downstream  boundary  conditions  on  the 
annulus  in  which  kt  applies. 

The  boundary  conditions  at  the  top  and  bottom  plates 
are  therefore  given  by 


(44) 


and 


u2 


(45) 


Assume  that 


~P“ 


9ZU) 


(46) 
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where 


*H)“ 


*(**)■■ 


From  (38),  (43),  and  (46),  it  can  be  shown  that 
(39)  will  be  satisfied  if 

2 

d  92  _  -iw92 
d22  =  v 


-  d  1 

f  dug  uj>  \ 

i  +  i!k 

\dr  r  ) 

dr  \ 

r  ) 

1  diZ 

Equations  (49)  and  (50)  can  be  integrated  directly 
and  the  boundary  conditions  in  (47)  and  (48)  can  be  imposed 
to  obtain 


.  coth(oz) 

coitifail) 


2Ai  ,1^2  I  f i) 

=  k'(e,CO,h2‘oC,d^ 


where  a  -  / 

Now  consider  (51)  for  the  case  in  which  k1  =  kQ  =  1. 
Because  this  is  a  linear  partial  differential  equation  which 
has  coefficients  constant  with  respect  to  z  and  which  has  only 
one  term  that  involves  differentiation  with  respect  to  z,  it 
can  be  separated  into  two  ordinary  differential  equations  if 
u3  (r,z)  is  chosen  to  be  a  product  of  a  function  f (r)  and  a 
ftmction  'I'(z) ,  or  more  generally,  a  summation  of  such  products. 


Since  boundary  conditions  that  are  symmetric  with  respect  to 
the  centerplane  will  be  imposed  at  a  given  radius,  the  proper 
form  of  u3  io 


u2  ”  £Cn(Wft  + 


(54) 


where  (ft)n  and(f2)n  are  the  two  independent  solutions  of  the 
equation  governing  fn  and  is  the  solution  of  the  equa¬ 
tion  that  i.;  symmetric  with  respect  to  the  centerplane.  The 
function  (f,)n  is  chosen  as  the  solution  that  specifies  the 
manner  in  which  information  propagates  in  the  downstream  di¬ 
rection,  while  (f2)n  specifies  the  upstream  propagation. 

Equation  (54)  does  not  apply  exactly  in  the  region 
near  the  drain,  because  the  drain  radius  is  finite  (implying 
that  the  radial  velocity  profile  is  not  exactly  uniform) . 
However,  this  solution  for  (f2)  should  accurately  describe 
how  information  propagates  upstream  in  this  region.  It  can 
be  used  to  show  approximately  how  small  the  drain  tube  should 
be  to  eliminate  its  effect  on  the  pickoff,  but  will  otherwise 
be  ignored  in  keeping  with  the  initial  description  of  the 
problem. 


Consider  then  only  the  effect  of  information  that  pro 
pagates  in  the  downstream  direction  from  the  upstream  solid 
boundary,  the  coupler.  If  only  the  effect  of  the  '"oupler  (as 
well  as  the  top  and  bottom  plates)  on  the  flow  in  the  chamber 

is  considered,  then  the  term  (•— - -•  +  -^2-/  in  equation  (51) 

3r  3r  r 

can  be  ignored.  This  will  be  verified  later  for  the  case  in 
which  ki  =  k_  (that  is,  when  the  radial  profile  is  uniform) 
by  shovjihg  that  the  solution  for  (fj)n  obtained  by  ignoring 
this  term  is  very  nearly  equal  to  the  solution  obtained  when 
it  is  included.  In  addition,  since  it  can  also  be  shown  that 
the  ratio 


a  /&*2  ,  uzA  j 

dr  U  Jll  = 

of  1 

is  a  maximum  for  the  case  in  which  kj  =  kD,  this  term  can  be 
neglected  in  the  general  case.  If  tuis  term  is  neglected, 
then  (54) ,  with  Dn  =  0,  represents  the  general  form  of  the 
solution  of  (51)  for  all  k,  .  The  governing  equations  are 


III  POnl 
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(55) 
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where  Pn  is  a  separation  constant.  The  boundary  conditions 
on  4>n  are 


It  should  be  noted  that  the  ordinary  differential 
equations  (49),  (50),  (55)  and  (56)  have  been  obtained  just 
by  assuming  the  form  of  the  solution  of  (39),  that  is,  by 
using  (43),  (46)  and  (54).  However,  the  partial  differen¬ 
tial  equation  (39)  is  linear  and  the  boundary  conditions 
specified  by  (40) ,  (41)  and  (42)  are  such  that  the  problem 
is  well-posed.  Then  the  complete  solution  thus  obtained  is 
unique  if  it  satisfies  the  appropriate  boundary  conditions, 
because  existence  is  a  necessary  and  sufficient  condition 
for  uniqueness  in  this  class  of  problems. 


Equation  (55)  has  the  solution 


where  an  arbitrary  multiplicative  constant  of  this  solution 
has  been  absorbed  into  ^n. 


tion 


It  is  clear  that  in  a  region  near  r1#  u2  has  a  solu- 

§  v^n'Z 
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and  that  this  is  an  eigenvalue  problem.  The  eigenfunctions 
*n  depend  on  kt  in  that  region.  The  eiaenvalues  can  be  found 
by  imposing  the  boundary  conditions  given  by  (57) .  The  co¬ 
efficients  cn  can  be  found  by  requiring  the  solution  to 
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satisfy  a  boundary  condition  at  the  upstream  limit  of  the 
region  in  which  k,  holds. 

It  should  be  pointed  out  that  this  eigenvalue  problem 
is  not  Sturm-Liouville.  The  eigenvalues  are  complex  and  the 
differential  operator  is  not  self-adjoint  since  a2  in  (56)  is 
imaginary.  Nevertheless,  it  is  still  possible  (App  B)  to 
use  typical  orthogonal  expansion  techniques  in  finding 
Cn.  However,  before  the  integration  is  performed  to  find 
C_,  both  sides  of  the  boundary  value  equation  (obtained  from 
(43),  (59)  and  the  upstream  boundary  condition)  must  be  multi¬ 
plied  by  k  instead  of  the  complex  conjugate  of  this  quantity, 
which  is  the  usual  procedure  with  a  Sturm-Liouville  problem. 


3 


-a 


2.2.2  .  Solution  for  the  Entire  Chamber 

An  approximate  solution  for  Vg  at  the  pickoff  can  be 
found  through  a  stepwise  procedure.  The  procedure  is  started 
by  specifying  the  functions  A,  and  k,near  the  coupler  and 
the  range  over  which  these  functions  will  be  applied.  The 
boundary  condition  at  RQ  given  by  (40)  is  used  to  find  the  co¬ 
efficients  of  the  series  solution  for  u2  valid  in  that  region. 

The  solution  is  then  found  for  the  lower  limit  of  r  at  which 
these  expressions  for  A,  and  k, apply.  This  function  is  taken 
as  the  upstream  boundary  condition  for  the  next  annulus  in 
which  a  new  pair  of  expressions  for  Aj  and  kj  apply,  and  the 
process  continues  until  a  solution  is  obtained  in  the  vicinity 
of  the  hypothetical  pickoff  location.  When  this  solution  is 
found,  the  characteristics  of  the  rate  sensor  can  be  described. 

It  should  be  pointed  out  that  the  tangential  shear 
is  not  continuous  across  the  boundaries  of  the  annuli.  Since 
(f2)n  has  been  ignored,  the  shear  is  not  continuous  if  the 
velocities  are  matched.  However,  the  error  can  be  made  ap¬ 
propriately  small  by  choosing  annuli  that  are  sufficiently  thin 
in  the  r  direction. 

According  to  assumptions  (3)  and  (4)  made  in  Section  2.1, 
the  output,  4> ,  of  the  hypothetical  pickoff  is  given  by 


A  |{t^/0  '  [(V®r«l  +  Vr2)  *in{orcton  ^f}J  drjj- 


(60) 


where  Vgre^  is  the  tangential  velocity  relative  to  the  wall 
and  Kj  is  some  constant.  Since  |Vfll«  I  V^l  and  Vg  is  a  linear 
function  of  «0f  K  h/2r  1 

Yp  ,v,  a  (61) 

1  rW  JrrR. 


33 


^jj 


9 

| 


1 

•3 


I 

1 

3 


I 


s 

A 

m 

% 


where  zt  is  the  pickoff  location  in  the  drain. 

If  experimental  data  obtained  from  a  sensor  with  a 
pickoff  in  the  drain  is  compared  to  expressions  for  the  out¬ 
put  of  the  hypothetical  pickoff  in  the  mathematical  model, 
then  changes  in  K2  during  an  experiment  can  lead  to  an  error 
that  is  not  related  to  the  accuracy  of  the  analysis  in  pre¬ 
dicting  the  features  of  the  flow  in  the  chamber,  the  domain 
in  which  this  solution  applies. 

In  the  real  sensor,  K2  will  be  influenced  by  changes 
in  the  fluid  losses  in  the  corner-drain  region  brought  on  by 
changes  in  the  fluid  parameters,  and  by  any  dependence  that 
the  streamline  pattern  may  have  on  the  Mach  number.  It  should 
be  remembered  that  one  purpose  of  this  approach  is  to  estab¬ 
lish  the  significance  of  these  effects  by  ignoring  them  in  the 
analysis.  Thus,  in  the  theoretical  results  it  is  necessary 
to  assume  that  K  is  some  appropriate  constant  without  any 
attempt  at  justification.  Because  of  the  assumption  used  to 
obtain  (62) ,  this  is  equivalent  to  assuming  that  K2  is  constant. 

By  comparing  the  output  at  u  =  0  to  that  vhich  would 
exist  if  the  flow  in  the  chamber  were  inviscid  (and  u=0) ,  Ki 
can  be  eliminated  from  the  expression  for  d-c  gain.  From  (4) , 
(5),  and  (60),  the  expression,  Gq,  for  d-c  gain  would  be  given 
by 
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(63) 


Even  though  this  expression  is  1  for  large  Re^  (assuming  that  a 
solution  over  the  entire  range  of  Re^  were  known) ,  this  does 
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not  imply  that  the  efficiency  of  the  sensor  is  100  percent  in 
that  range  of  Re0',  for  the  losses  in  the  drain  region  could 
be  significant.  The  value  of  (63)  is  independent  of  such  losses 

If  K]  is  chosen  so  that  is  1  at  «  =  0,  another  func¬ 
tion,  ,  is  obtained  that  describes  the  frequency  response 

of  the  sensor. 


2.3  Bounding  Expressions  for  the  Stepwise  Solution 

The  method  described  here  will  involve  a  great  deal 
of  effort,  principally  because  of  the  stepwise  procedure  out¬ 
lined.  Before  this  task  is  undertaken,  the  importance  of  the 
changes  in  kj  with  respect  to  r  should  be  established.  This 
can  be  done  by  solving  tv;o  hypothetical  problems,  one  in  which 
the  radial  velocity  profile  is  assumbed  parabolic  throughout 
the  chamber,  and  one  in  which  it  is  assumed  uniform. 


It  can  be  shown  (App  C)  that  these  two  solutions  will 
provide  upper  and  lower  bounds  for  the  results  obtained  by  the 
lengthy  stepwise  procedure  for  at  least  a  portion  of  the 
ranges  of  the  governing  parameters.  The  parabolic  solution 
can  be  shown  to  provide  an  upper  bound  for  G0,  and  the  uniform 
one,  a  lower  hound.  To  discuss  the  nature  of  the  bounds 
these  solutions  provided  for  G^u),  the  contributions  of  ui 
and  U2.  tq  ..the  output  signal  must  be  discussed  separately.. 

For  all  w'  (where  u-'=u)h2/v) ,  the  parabolic  solution 
provides  a  maximum  expression  (and  the  uniform  a  minimum',  for 
phase  lag  in  the  signal  due  to  ui.  At  small  w' ,  the  amplitude 
of  Gio)  due  to  ui  depends  only  weakly  on  the  shape  of  the  rad¬ 
ial  profile,  although  it  does  appear  that  the  parabolic  solu¬ 
tion  provides  a  maximum  expression  for  amplitude,  and  the 
uniform,  a  minimum.  As  id'  increases,  the  effect  on  the  am¬ 
plitude  caused  by  the  variation  of  the  radial  profile  eventual¬ 
ly  becomes  significant.  For  large  oy' ,  the  uniform  solution 
provides  a  maximum,  and  the  parabolic,  a  minimum  for  the 
amplitude  of  G^u  due  to  uj. 

For  the  portion  of  G^w  due  to  U2,  the  uniform  solu¬ 
tion  provides  a  maximum  expression  for  phase  shift,  and  a 
minimum  for  amplitude.  The  parabolic  provides  a  minimum  for 
phase  shift  and  a  maximum  for  amplitude. 

Several  interesting  features  of  G^  should  be  pointed 
out.  The  amplitude  of  the  contribution  of  Uj  to  G^  diminishes 
to  zero, and  the  phase  lag  approaches  some  constant  as  the  fre- 
cruency  increases.  The  amplitude  of  the  contribution  of  U2  to 
G  •  u  approaches  a  constant  roughly  4/Rep  times  that  at  dc,  and 
tne  phase  lag  increases  without  limit. 
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If  Re£  is  sufficiently  small,  the  dc  signal  due  to 
U2  is  attenuated  so  much  as  it  is  transported  from  the  coup¬ 
ler  to  the  pickoff  that  its  contribution  is  below  the  noise 
level  for  all  to' .  Then  only  the  contribution  dv>e  to  ui  need 
be  considered  and  the  hypothetical  solutions  provide  bounds 
as  specif led  .above. 

At  small  (ji' ,  if  Red,  is  not  sufficiently  small  to 
wipe  out  the  contribution  due  to  U2  but  small  enough  to  ensure 
that  the  contribution  due  to  ui  predominates,  then  either  ex¬ 
pression  for  the  amplitude  may  be  used  because  of  their  slight 
difference,  while  a  maximum  value  for  phase  shift  can  be 
found  from  the  parabolic  solution  and  a  minimum  from  the 
uniform. 

At  sufficiently  large  Re^ ,  the  U2  contribution  pre¬ 
dominates  and  specifies  the  nature  of  the  bounds  these  ex¬ 
pressions  provide.  At  intermediate  values  of  Re£  when  u'  is 
small,  both  contributions  are  of  nearly  equal  magnitude,  and 
the  behavior  is  mixed. 

As  id'  increases,  the  value  of  Reo  at  which  the  be¬ 
havior  becomes  mixed  decreases,  as  does  the  value  at  which 
the  contribution  due  to  U2  predominates.  At  sufficiently 
large  w',  u2  predominates  for  all  values  of  Re0  sufficiently 
large  to  insure  that  its  amplitude  is  not  attenuated  below 
the  noise  level. 

This  information  has  been  summarized  in  the  repre¬ 
sentation  in  figure  5,  with  some  added  information  based  on 
a  number  of  comparisons  of  the  solutions  to  be  obtained  in 
2.3.1  and  2.3.2. 


Zone  i  represents  the  region  in  which  ui  predominates. 
Either  solution  can  be  used  to  express  the  amplitude  of  Giw  at 
small  a)'.  At  larger  w',  the  difference  between  these  two  ex¬ 
pressions  for  amplitude  is  not  su.all .  In  that  range  of 
the  parabolic  solution  expresses  a  lower  bound  for  amplitude, 
and  the  uniform,  an  upper  bound.  The  parabolic  solution  ex¬ 
presses  a  maximum  for  phase  shift, and  the  uniform, a  minimum 
for  all  w'. 

Zone  2  is  the  region  in  which  the  contribution  due  to 
ui  and  U2  must  both  be  considered,  but  in  which  the  phase  lag 
of  either  is  small.  If  the  difference  in  amplitude  or  phase 
of  the  two  solutions  is  significant,  they  may  be  regarded  as 
bounding  expressions.  Otherwise,  either  expression  may  be 
used. 

Zone  3  represents  the  only  region  in  which  these  two 
solutions  cannot  be  used  directly  to  bound  the  amplitude  and 
phase  of  G£u.  In  this  region,  the  contributions  due  to  Uj  and 
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Figure  5.  Boundary  zones  of  hypothetical  solutions. 


U2  can  interfere  destructively,  since  they  are  nearly  of  equal 
amplitude  but  out  of  phase  by  large  margins.  The  amplitude 
and  phase  can  change  radically  with  to'  in  this  zone.  This 
behavior  is  reflected  in  both  solutions ,  but  the  values  of 
w'  and  Reo  a  -  hich  similar  phenomena  occur  are  different  in 
the  two  soli  I  ns.  This  implies  that  the  two  solutions  can 
be  used  to  isolate  the  range  of  Re^  at  which  a  specific  fre¬ 
quency  response  feature  occurs  (such  as  sudden  reduction  and 
subsequent  recover/  of  amplitude  of  as  in'  increases) ,  as 
well  as  the  range  of  to'  in  which  this  behavior  occurs.  How¬ 
ever,  the  two  solutions  at  a  specific  Re,!,  and  w'  cannot  be 
used  in  a  direct  manner  to  estimate  the  actual  behavior. 

The  solution  at  neighboring  values  of  Re^  and  oj"  is  required 
for  this. 


The  exact  range  of  zone  3  has  not  been  accurately 
specified  at  this  time,  but  it  appears  to  extend  from  approx¬ 
imately  .5  to  3  on  the  Re^  axis  when  w'  =  15. 


Zone  4  represents  the  region  in  which  u2  predomi¬ 
nates.  In  this  region,  the  parabolic  solution  provides  an 
estimate  for  the  maximum  amplitude  of  and  the  minimum 
phase  lag.  The  uniform  solution  has  the  complimentary 
features. 
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2.3.1  Solution  for  a  Parabolic  Profile 


If  (33)  and  (34)  are  assumed  to  apply  throughout  the 
chamber,  it  can  be  shown  (App  D)  that 
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M  [at  ,  a2  ,x].  is  the  confluent  hypergeometric  function,  and  Sn 
is  a  member  of  the  ordered  set  of  values  which  satisfy  the 
equation 
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The  prescript  "o"  used  in  (64)  is  used  to  denote  the  values 
of  1^,  ,  and  Xn  when  is  zero. 

The  computer  program  used  to  numerically  evaluate 
(64)  and  (65)  is  presented  in  Appendix  E.  The  series  in  both 
(64)  and  (65)  have  been  truncated  after  the  first  term  in  the 
evaluation  presented.  The  reason  for  this  can  be  discussed 
more  clearly  when  the  solution  for  a  uniform  radial  profile 
has  been  shown. 


2.3.2 


Solution  for  a  Uniform  Profile 


If  (36)  and  (37)  are  used/  it  can  be  shown  (app  F) 
that  the  functions  that  describe  the  gain  and  frequency  re- 
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In  these  equations  both  of  the  infinite  series  expressions,  which 
have  been  generated  by  u2,  can  be  truncated  after  the  first  term. 
The  reasons  for  this  deserve  elaboration. 

When  R€q  is  of  the  order  of  1,  the  magnitude  of  the 
exponential  term  diminishes  rapidly  as  n  increases.  This  re¬ 
sults  from  the  fact  that  the  higher  harmonics  of  are  atten¬ 
uated  more  than  the  fundamental  as  the  signal  propagates  from 
the  coupler  to  the  pickoff  due  to  the  higher  viscous  shear 
they  generate. 

In  addition,  the  contribution  of  each  mode  due  to  the 
boundary  condition  at  the  coupler  varies  with  the  order  of 
the  mode.  Since  the  function  that  the  series  must  satisfy  at 
the  coupler  when  to"  is  small  has  a  single  maximum  and  no  in¬ 
flection  points,  the  fundamental  mode  contributes  most  strongly 
to  the  series  representation  at  small  to".  In  (78),  the  co¬ 
efficient  of  the  exponential  for  n=0  is  81  times  that  for  n=l. 
This  is  nearly  the  same  as  in  the  series  in  (64).  Thus,  for 
small  to",  even  for  that  range  of  Reo  in  which  the  higher  order 
modes  are  not  greatly  attenuated  between  the  coupler  and  the 
pickoff,  their  contributions  are  not  significant. 

As  to"  increases,  the  significance  of  the  higher  order 
modes  increase,  but  never  predominates.  For  very  large  to", 
the  coefficient  in  (79)  for  n=0  is  9  times  that  for  n=l. 

Again,  the  behavior  is  nearly  the  same  In  (65).  It  might  be 
argued  that  this  significance  is  not  sufficiently  weak  to 
warrant  truncation.  However,  in  that  r*xnge  of  to",  the 
amplitude  '  is  so  much  less  than  the  d-c  value  that  the 
exact  features  are  not  of  great  interest  to  a  control  system 
designer.  Therefore,  the  series  can  be  truncated  after  the 
first  term  for  all  Re^  and  u"  within  the  intent  of  this  in¬ 
vestigation. 

No  program  is  presented  for  the  solution  of  (78)  and 
(79)  because  of  their  rather  simple  nature.  It  should  be 
noted  however,  that  double  precision  is  required  in  evaluating 
the  bracket  containing  the  hyperbolic  tangent  terms  when  w" 
is  very  small. 
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The  error  in  the  solution  for  u2  caused  by  neglecting 
this  viscous  term  will  be  larger  if  Re^  is  very  small,  because 
more  terms  of  the  series  in  (80)  must  then  be  considered, 
making  the  approximation  for  the  ratio  of  the  gamma  functions 
less  accurate.  However,  when  Re^  is  very  small,  the  contri¬ 
bution  of  Uj  to  the  solution  for  Vg  near  the  pickoff  is  negli¬ 
gible.  Thus  the  error  generated  in  GQ  and  G-  by  neglecting 
the  term  3/3r(3u2/3r  +  u2/r)  is  small  ror  all  Re^  in  the  fully 
developed  range  when  the  radial  profile  is  uniform.  As  in¬ 
dicated  in  Section  2.2.1,  this  implies  that  the  error  is  small 
in  general. 


2.4.2  Error  Due  to  Neglecting  Entrance  Effects  on  the 

Radial  Velocity*  Profile 

It  was  assumed  that  the  entrance  effects  did  not  in¬ 
fluence  the  radial  profile  throughout  the  chamber.  In  the 
actual  sensor  the  radial  profile  should  be  significantly  in¬ 
fluenced  hy  the  upstream  boundary  conditions  in  the  outer 
portion  of  the  chamber  for  the  values  of  Reo  near  7.  When 
Re©  is  decreased  by  changing  h  in  the  upper  portion  of 
the  fully  developed  range,  the  decrease  in  Vq  in  the 
actual  sensor  will  not  be  as  large  as  that  predicted  by  this 
analysis,  since  the  viscous  losses  in  the  actual  model  are  not 
as  large.  When  Re^  is  decreased  in  the  lower  portion  of  the 
fully  developed  range,  however,  the  proportionate  decrease  of 
Vq  should  be  predicted  accurately,  since  the  range  of  r  in 
wnich  the  profile  is  not  fully  developed  is  quite  small.  Con¬ 
sequently,  if  the  theoretical  curve  and  the  data  are  normalized 
so  that  they  intersect  when  Re^  is  7,  the  slope  of  the  theoret¬ 
ical  curve  of  G0  as  a  function  of  Re,!,  in  the  upper  portion  of 
the  fully  developed  range  should  not  be  as  flat  as  that  of  the 
empirical  curve. 


2.4.3  Effects  of  the  Corner  and  Drain  on  Sensor  Characteristics 

It  is  shown  in  Appendix  C  that  if  k,  =  kD,  the  func¬ 
tion  (f2)  which  describes  the  propagation  of  the  information 
in  the  upStream  direction  is  given  by 


~R|q 

(*2)n=  f  ^  kR«0  (•'/iu'  +  (2n+  -jj-) 

~T~' 


(83) 
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where  Re0  =  UqRqA  =  Q/2irhi>,  and  K^Re  /2)-i  is  the  modified  | 

Bessel  function  of  the  second  kind  of  order  (Re0/2)-l.  | 

Just  upstream  of  the  drain,  (f2)n  behaves  nearly  as  | 

l"Re  -  S 

r  o.  since  ReQ  is  on  the  order  of  104  in  this  experiment,  J 

drain  effects  are  extremely  localized.  The  corner  in  the  | 

real  sensor  can  then  be  defined  as  the  region  extending  from  J 

a  radius  just  upstream  of  the  drain  to  a  plane  sufficiently  I 

far  downstream  in  the  drain  to  ensure  that  the  streamlines  lie  i 

in  a  cylindrical  surface  concentric  with  the  drain  wall.  J 

Consider  what  would  happen  if  (f2)a  were  retained  J 

in  the  chamber  solution,  and  the  corner  region  of  a  two  drain  f 

sensor  were  instead  described  simply  as  a  curve  F(r,z)  as  J 

shown  in  figure  6.  The  curve  F(r,z)  is  chosen  so  that  con-  | 

tinuity  is  satisfied  locally  when  the  two  transport  velocity  f 

profiles  are  matched.  § 


-I 

1 
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Figure  6 .  Corner  Region  Model 


Since  the  fluid  particles  approach  F(r,z)  in  a  plane 
parallel  to  the  chamber  walls  and  leave  on  a  surface  parallel 
to  the  drain  wall,  this  description  of  the  acceleration  ef¬ 
fects  on  the  transport  velocity  is  clearly  nonphysical.  How¬ 
ever,  this  description  can  be  used  to  discuss  to  a  first-order 
approximation  the  total  viscous  losses  in  the  tangential  com¬ 
ponent  near  the  corner. 
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Note  first  that  in  the  tangential  momentum  equation 
the  acceleration  across  this  curve  is  zero  if  there  is  no  dis¬ 
sipation.  The  nonphysical  description  of  the  transport  com¬ 
ponent  of  velocity  influences  the  tangential  equation  only  in 
that  it  arbitrarily  sharpens  the  manner  in  which  the  transport 
velocity  changes  from  the  chamber  function  to  the  drain  func¬ 
tion.  If  the  change  is  smooth  and  monotonic  in  the  actual 
case,  as  it  is  near  the  points  (R  ,  ±  h/2) ,  then  the  error  is 
second  order.  If  the  fluid  stagnates  in  going  from  one  region 
to  the  next,  as  it  does  near  (0,  0) ,  then  the  error  may  not 
be  small.  This  approach  can  then  be  used  to  examine  the  vis¬ 
cous  loss  in  the  tangential  component  near  (Rw  .+  h/2). 


The  tangential  velocity^ in  the  chamber  in  the  region 
naar  F(r,z)  is  given  approximately  by 


ve  =  [(t  +  0  «2  -  £ cn(f,)n*n  + 1  Go*4'-'* 


(84) 


where  g, ,  g2  ,  (f^n  and  *n  are  iound  by  setting  kt  =  kD. 
In  the  drain  the  tangential  velocity  is  given  (app  H)  by 
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and  m  =  mass  flow  rate.  Ii  is  the  modified  first-order  Bes¬ 
sel  function  of  the  first  kind  and  the  5n's  are  the  ordered 
set  of  values  that  satisfy  the  relationship 


The  coefficients  CQ  are  specified  by  the  boundary 
conditions  at  some  boundary  upstream  of  the  drain  (since  (f*)n 
is  effectively  zero  just  upstream  of  Ri)  while  and  Dn  are 
chosen  so  that  the  tangential  velocity  and  the  shear  are  con¬ 
tinuous  across  F(r,z) .  Some  general  features  of  this  matching 
can  be  discussed  without  going  into  explicit  detail. 

Consider  the  tangential  velocity  profile  in  the  chamber 
sufficiently  far  upstream  of  F(r,z)  to  eliminate  the  influences 
of  (f2)n.  If  the  tangential  profile  decreases  evenly  from  a 
maximum  at  the  centerline  to  a  minimum  at  the  wall,  then  the 
values  found  for  the  coefficients  Dn  are  small,  and  the  coef¬ 
ficients  Cn  are  nearly  equal  to  those  that  are  chosen  simply 
to  match  Vq  at  F(r,z).  However,  if  the  profile  is  significantly 
blunted,  the  shear  generated  by  this  matching  condition  will 
be  large. 

In  the  real  sensor,  the  matching  region  may  well  span 
most  of  the  region  I  z  I  <  h/2,  r  <  R,  ,  as  well  as  include  a  sig¬ 
nificant  stagnation  region  about  the  point  (0,0)  and  a  separa¬ 
tion  region  around  the  corner.  In  addition,  the  drain  may 
distort  the  radial  profile  upstream  of  the  corner  and 
cause  a  small  V2  to  exist.  However,  to  a  first  order  ap¬ 
proximation  this  distortion  will  not  change  the  axisymmetric 
propagation  of  tangential  shear  along  a  given  streamline.  The 
■  field  will  only  be  distorted.  Therefore,  the  net  effect  of  the 
viscous  dissipation  at  the  corner  on  the  output  characteristics, 
should  be  about  the  same  as  that  predicted  by  this  description, 
that  is,  that  the  losses  'ear  the  corner  will  be  relatively 
high  if  the  tangential  velocity  profile  is  somewhat  blunted. 

Consider  what  happens  for  the  case  when  'J  =  0  as  Re^ 
is  reduced  (by  decreasing  h,  for  instance)  from  a  large  value 
to  the  value  Reocrit'  the  point  at  which  the  flow  becomes 
fully  developed.  The  tangential  velocity  profile  near  the 
corner  changes  from  nearly  uniform  xfor  nearly  inviscid  flow) 
to  nearly  parabolic  (gj  at  u  =  0  is  the  parabolic  profile  if 
hi  *  k  ) .  The  bluntness  of  the  tangential  profile  is  then 
reduced  considerably.  Although  a  pickoff  mounted  in  the  chamber 
would  see  a  decrease  in  G0  as  Re’  is  reduced,  due  to  the  in¬ 
creased  chamber  losses,  one  mounted  in  the  drain  would  see  at 
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best  a  small  change.  Because  little  of  the  information  des¬ 
troyed  by  the  boundary  layer  can  fce  recovered  by  a  pickoff  in 
the  drain  due  to  the  bluntness  of  the  tangential  profile, 
changing  the  relative  boundary  layer  size  should  not  greatly 
change  the  sensor  gain  relative  to  the  inviscid  case. 

Consider  next  the  effect  of  the  corner  losses  ->n  the 
frequency  response  of  the  sensor.  Since  the  magnitude  of 
g2  (z)  near  the  wall  relative  to  the  value  at  the  center  plane 
increases  as  c 3  increases,  the  profile  is  reduced  more  at  the 
center  plane  with  increasing  cJ  .  Because  this  effect  leads 
directly  to  a  blunting  of  the  tangential  profile,  the  pickoff 
will  recover  a  smaller  portion  of  the  total  time-dependent  in¬ 
formation  as  <J  increases.  This  implies  th&t  the  attenuation 
predicted  by  this  analysis  will  not  be  as  large  as  that  measured 
with  the  experimental  model.  As  cJ  increases  the  pickoff  will 
also  recover  a  smaller  portion  of  the  information  near  the  wall, 
which  is  in  the  closest  phase  agreement  with  the  sensor  motion. 
This  implies  that  the  phase  shift  predicted  by  this  analysis 
will  not  be  sufficiently  large. 
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It  is  possible  to  make  two  other  observations  about 
the  losses  in  the  corner-drain  region.  The  magnitude  of  the 
lower  order  eigenvalues  of  (90)  is  of  the  order  of  10.  Since 
Re!  is  of  the  order  of  104  in  this  experiment,  ^  is  very 
nearly  one  for  the  lower  order  eigenfunctions  ifz  is  of  the 
order  of  Rj .  In  addition,  the  term  6^* /Re2,  in  the  first  co¬ 
efficient  of  the  hypergeometric  function  can  then  be  ignored. 
Only  a  few  terms  should  be  needed  to  describe  the  tangential 
profile  in  the  drain,  so  the  first  result  states  that  the 
losses  in  the  drain  tube  can  be  ignored.  The  second  result 
implies  that  the  losses  in  the  corner  region  are  independent 
of  Rej ,  for  this  parameter  does  not  appear  in  either  (84)  or 
in  the  simplified  version  of  (85) .  This  latter  result  is  im¬ 
portant  and  will  be  discussed  in  the  experimental  phase  of 
the  study. 


3.*  EXPERIMENT 

3.1  Approach 

It  is  possible  to  obtain  some  important  information 
about  the  viscous  effects  in  the  chamber  and  in  the  corner 
region  by  comparing  this  analysis  with  data  obtained  for  a 
properly  chosen  set  of  conditions. 


3.1.1 


Determination  of  G. 


The  accuracy  of  the  analysis  in  predicting  GQ  when 
only  the  chamber  losses  are  varied  is  first  examined.  Data 
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relating  G0  to  Re^>  was  obtained  by  changing  h,  since  this  is 
the  only  parameter  that  can  be  easily  varied  in  Re^  to  cause 
the  gain  to  approach  some  constant  at  large  Re^  (for  the  pur¬ 
pose  of  normalizing  the  data  to  compare  it  with  theory) . 

Viscosity-  could  also  be  varied,  but  such  an  experiment  would 
surely  be  unmanageable. 

Since  Rei  is  independent  of  h,  the  accuracy  of  pre¬ 
dicting  chamber  losses  can  be  tested  independently  of  testing 
the  prediction  that  corner  losses  are  not  a  function  of  Rei . 

If  the  previous  test  is  repeated  for  various  values  of 
A  (and  thus  various  values  of  Rei),  then  curves  can  be  obtained 
that  examine  ranges  of  Re^  that  have  differing  upper  limits. 

The  degree  of  correlation  of  these  curves  would  reveal  some 
information  about  the  dependence  on  Rei  of  the  losses  at  the 
corner,  as  will  be  discussed  later. 

3.1.2  Determination  of  G,  G^w,  and  the  Corner  Losses  due 

to  Ret  ~~ 

3. 1.2.1  Influence  of  Re,  on  G  and  Gj(t} 

It  is  of  great  interest  to  examine  the  accuracy  of  this 
analysis  in  predicting  G^w  as  well  as  the  dimensional  d-c  sensor 
gain  (that  is,  not  divided  by  inviscid  gain) ,  G,  if  Q0  and  v 
are  varied.  But  since  these  latter  variations  would  lead  to  a 
change  in  Rei,  it  is  necessary  to  carefully  choose  the  experi¬ 
mental  conditions  in  order  that  accurate  conclusions  can  be  drawn 


In  the  range  Reo<  1,  even  the  first  terms  in  the 
series  portions  of  the  solutions  for  GQ  and  Gj.w  can  be 
ignored.  This  is  equivalent  to  stating  that  the  contributions 
due  to  the  entrance  conditions  are  negligible,  or  that  the 
features  of  Vq  near  the  pickoff  are  specified  by  u, . 

It  will  now  be  shown  that  in  this  range  of  Re'  not 
only  can  one  conduct  tests  that  predict  G  for  changing°Q0 
and  v  as  well  as  tests  that  predict  G^  ,  but  that  these  same 
tests  can  investigate  in  different  manners  whether  the  losses 
at  the  corner  depend  on  Re,  .  It  is  not  unreasonable  to  re¬ 
strict  the  range  of  Re^  to  values  less  than  one,  particularly 
in  the  first  planned  experimental  study  of  the  fully  developed 
range,  since  this  limited  range  is  just  that  which  was  pre¬ 
viously  specified  by  intuitive  means  to  the  most  advantageous 
range,  the  range  in  which  entrance  conditions  have  no  effect 
on  signal  output. 

3. 1.2. 2  G^w  and  the  Depertdence  of  the  Corner  Losses  on  Re, 

Suppose  that  an  experiment  were  conducted  in  this 
range  of  Re£  to  predict  G^w  for  various  values  of  Re^  obtained 


by  varying  ft.  Because  Re,  is  also  a  function  of  ft,  any  in¬ 
fluence  of  Re,  on  the  losses  at  the  corner  could  appear  and 
confuse  the  results. 


For  each  value  of  Re^,  data  would  be  reduced  for 
comparison  to  the  analytical  expression  for  G^w  by  comparing 
the  magnitude  of  the  signal  at  various  values  of  'J  to  that 
d-c.  Thus,  this  data  would  not  be  influenced  by  the  magnitude 
of  u,  at  d-c  or  any  dependence  on  Re,  of  the  manner  in  which 
the  magnitude  of  the  tangential  velocity  profile  is  altered 
at  the  corner,  since  those  dependences  are  Eliminated  when 
the  signal  is  compared  to  d-c. 


However,  the  data  could  be  influenced  by  any  de¬ 
pendence  on  Re,  of  the  manner  in  which  the  shape  of  the  tan¬ 
gential  profile  is  altered  at  the  corner.  It  has  already 
been  predicted  in  section  2.4.3  that  the  losses  at  the  corner 
depend  upon  the  bluntness  of  the  profile.  Because  the  shape 
of  u,  is  a  function  of  -J  ,  an  error  should  be  expected  in  the 
analytical  expression  for  that  is  a  function  of  <</  .  How¬ 
ever,  since  the  shape  of  u,  is  not  a  function  of  Re^,  the  shape 
of  the  velocity  profile  at  the  corner  is  the  same  for  a  given 
(J  for  all  values  of  Re,  tested.  As  a  result,  the  shape  ef¬ 
fects  at  the  corner  are  independent  of  Re,  if  the  error  in 
predicting  G-u  is  a  function  of  only  tJ  for  values  of  Re^  less 
than  one  that  are  obtained  by  changing  ft. 


3. 1.2. 3  G  and  the  Dependence  of  the  Corner  Losses  on  Re, 


For  Re'  <1,  Gq  is  proportional  to  Re^  in  both  the 
uniform  and  parabolic  radial  profile  solutions.  It  can  then 
be  shown  that  this  analysis  predicts  that,  for  a  pickoff  in 


the  drain  and  Re^  <1,  G  is  given  by 


where  K3  is  a  function  of  pickoff  design  and  location  as  well 
as  the  losses  at  the  corner. 


It  should  be  recalled  that  only  the  magnitude  of  u, 
varies  if  ft  and  n  are  changed  for  w  =  0.  The  outcome  of  an 
experiment  conducted  to  predict  G  will  then  be  influenced  by 
the  dependence  on  Re  of  the  manner  in  which  the  magnitude  of 
the  tangential  profile  is  altered  by  the  viscous  losses  at  the 
corner.  The  degree  to  which  K3  depends  on  Re,  can  be  demon¬ 
strated  by  tests  that  show  how  curves  for  G  versus  ft  correlate 
for  various  gases.  Since  Re, changes  for  a  given  m  when  u 
is  changed,  the  magnitude  effects  will  be  independent  of  Re, 
if  the  experimental  curves  correlate  when  (91)  is  used. 
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These  tests  should  show  whether  the  corner  losses 
are  independent  of  Re,  and  should  demonstrate  the  accuracy 
of  this  analysis  in  predicting  the  effects  of  any  design 
change  other  than  a  change  in  the  pi.koff-drain  configuration. 


>aratus 


Data  points  were  obtained  from  the  experimental  model 
shown  in  figure  4.  For  the  frequency  response  data,  the  cham¬ 
ber  height  was  fixed  at  0.035  in. 


Data  points  relating  G~  to  Re^  were  obtained  by  varying 
the  chamber  height  between  0.250  and  0.004  in.  Mass  flow  rate 
was  held  constant  in  obtaining  a  single  curve,  but  was  varied 
in  a  series  of  such  tests,  as  previously  discussed. 

'The  predictions  of  (91)  were  tested  by  varying  fh  at  a 
fixed  chamber  configuration.  The  gain  was  measured  using 
three  different  gases:  carbon  dioxide,  air,  and  helium. 


In  every  test  the  output  of  each  channel  of  the  pickoff 
could  be  monitored.  A  miniature  hot-wire  probe  was  used  to 
measure  flow  rate,  and  a  strain  gauge  transducer  was  used  to 
measure  pressure.  The  load  resistance  in  each  channel  could 
be  varied  by  changing  a  valve  downstream  of  the  transducer 
cavity. 


Only  the  data  points  for  an  infinite  load  resistance 
are  presented  in  this  paper,  because  it  was  decided  that  the 
characteristics  of  the  sensor  for  a  blocked  load  should  be 
analyzed  before  the  sensor's  load  sensitivity  is  discussed. 

With  an  infinite  load  resistance,  the  frequency  response  of 
the  readout  system  was  high  compared  with  that  of  the  rate  sen¬ 
sor.  Direct  measurement  with  the  pickoff  and  readout  assembly 
indicated  that  the  amplitude  was  down  about  1  dB  relative  to 
d-c  at  100  Hz. 
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4 .  RESULTS 


Data  and  Theory  for  G„ 


Data  points  relating  0o  to  Re,!,  are  shown  in  figure  7. 
Solutions  of  (64)  and  (78)  are  also  shown  in  this  same  figure. 
The  data  have  been  normalized  to  fall  midway  between  the  two 
theoretical  curves  so  that  the  shapes  of  these  curves  can  be 
compared.  The  actual  value  of  GQ  for  each  data  point  can  be 
obtained  by  multiplying  the  value  shown  by  1.49,  or  by  regarding 
the  value  of  the  data  point  at  Re'  =  7  as  having  a  value 
G0  =  1.00  (+  0.02) . 
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Data  obtained  at  a  very  high  Itej  (over  a  range  of  Re' 
up  to  20)  showed  that  the  curve  was  flat  within  an  experimental 
accuracy  of  several  percent  for  va.lues  of  Re^  above  7 .  The 
data  for  the  largest  Re,  shown  here  were  therefore  normalized 
with  respect  to  its  value  near  Re^  =  7,  the  largest  value  of 
R6q  that  could  be  obtained  at  that  value  of  Rej . 


Since  the  curves  for  the  smaller  values  of  Re,  did  not  have 
that  great  a  range  of  Re^,  they  were  normalized  to  fit  the 
curve  for  the  largest  Re,  at  their  outer  data  point. 


This  method  of  fitting  the  data  curves  together  is 
in  a  sense  artificial,  but  it  can  reveal  some  information. 

It  should  be  noted  that  the  shape  of  V#  does  depend  on  Re^ 
in  this  range,  since  u2  cannot  be  neglected.  Then  if  the 
manner  in  which  the  losses  near  the  corner  change  .  the  shape 
of  the  tangential  profile  were  to  depend  on  Re,  ,  the  shapes 
of  these  curves  would  be  different.  They  would  not  be  super¬ 
imposed  on  each  other  when  normalized  in  this  manner. 

The  information  that  is  lost  through  this  method  is 
the  dependence  on  Re  of  the  alteration  of  the  magnitude  of 
the  tangential  profile  near  the  corner, since  the  information 
about  the  magnitude  of  the  curve  is  lost.  However,  this  ques¬ 
tion  is  explored  in  the  next  comparison. 
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4 . 2  Data  and  Theory  for  G 

Two  series  of  tests  were  conducted  co  investigate 
the  dependence  of  K3  in  (91)  on  Rej  .  The  results  of  the 
first  series  are  shown  in  figure  8.  The  gain  of  the  sensor 
as  a  function  of  standard  volumetric  flow  rate  is  plotted 
for  two  different  gases:-  air  and  carbon  dioxide.  A  para¬ 
bolic  curve  was  drawn  to  intersect  the  empirical  curve  for 
air  somewhere  in  the  middle  of  its  range.  Deviations  of  the 
data  from  this  curve  can  be  considered  to  arise  front  a  change 
in  th»-  value  of  K3  .  This  parabolic  curve  and  (91)  were  then 
used  to  obtain  a  similar  curve  through  the  data  points  for 
carbon  dioxide  by  assuming  K3  is  constant. 


STANDARD  FLOWRATE  (cC/MIN  X I0*3) 


Figure  8.  Sensor  Gain  vs  Standard  Volumc-cric  Flow 
Rate  -  Air  and  Carbon  Dioxide. 
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Theory  indicates  that  K3  is  independent  of  Re, .  How¬ 
ever,  it  may  depend  on  some  other  parameter.  An  effort  is 
made  here  to  correlate  the  variations  in  K3  to  the  Mach  number 
at  the  exit  of  the  drain.  Since  the  pickoff  is  downstream  of 
a  sudden  enlargement  of  drain  diameter  in  the  model  tested 
(see  figure  4),  the  output  may  vary  as  the  Mach  number  changes 
due  to  the  associated  change  in  streamline  pattern. 

The  percentage  deviation  of  the  empirical  values  from 
the  parabolic  curves  is  plotted  for  both  gases  in  figure  9 
as  a  function  of  Mach  number  at  the  pickoff.  The  Mach  number 
was  calculated  by  assuming  the  flow  was  isentropic  and  the 
exit  pressure  was  one  atmosphere.  Therefore, 


Figure  9.  Percent  Deviation  of  Empirical  Values 
from  Parabolic  Curves  of  Figure  8. 
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where  7  *  ratio  of  specific  heats,  Pc  =  chamber  pressure, 
and  Pa  =  atomspheric  pressure.  Although  the  flow  in  the 
chamber  is  clearly  not  isentropic  in  this  case,  the  error  in 
(92)  should  not  be  large.  In  a  direct  test,  no  change  in  the 
chamber  pressure  for  fixed  flow  rate  was  detected  when  the 
chamber  height  was  varied  over  the  widest  possible  range,  in¬ 
dicating  that  the  acceleration  effects  dominate  in  vietermining 
the  pressure  distribution. 

This  test  was  repeated  with  helium  and  air,  as  shown 
in  figurelO  and  11 (note  scale  difference  between  figures  8 
and  10) .  The  shapes  of  the  curves  in  figures  9  and  11  are 
ditferent  because  the  pic)  off  position  was  changed  before  the 
latter  test  was  conducted,  causing  the  effect  brought  on  by 
a  changing  streamline  pattern  to  vary. 


FigurelO.  Sensor  Gain  versus  Standard  Volumetric 
Flow  Rate — Air  and  Helium. 
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Figure  11.  Percent  Deviation  of  Empirical  Values  from 
Parabolic  Curves  cf  Figure  10. 


4.3  Data  and  Theory  for  G.; 

Numerical  solutions  of  (65)  and  (79)  are  plotted 
with  the  data  in  figures  12  through  15.  The  difference  (in 
decibels)  between  the  magnitude  given  by  (65)  and  the  empirical 
amplitude  ratio  is  plotted  in  figure  16.  The  difference  in 
phase  is  plotted  in  figure  17.  The  magnitude  and  phase  of 
(65)  are  plotted  in  figures  18  and  19,  respectively,  for  a 
wider  range  of  Re^. 


Figure  12.  Numerical  Solutions  of  Equations  (65) 
and  (79) — Re^  =  0.273 
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Figure  19.  Phase  of  Equation  (65). 
DISCUSSION 


5.1  Prediction  of  GQ 

The  analysis  predicts  rather  well  the  behavior  of 
G0  as  a  function  of  Re^.  The  shapes  of  the. data  curve  and 
the  two  theoretical  curves  in  figure  7  are  quite  similar  ex¬ 
cept  that  the  empirical  curve  is  somewhat  flatter  in  the 
range  of  Re^  between  2  and  7,  as  anticipated  in  section  4.2. 

It  was  also  expected  that  the  empirical  curve  would  be 
nearly  flat  above  Reocrit,  even  though  the  chamber  losses 
change  in  that  region,  because  little  of  the  signal  destroyed 
by  the  boundary  layers  is  recovered  by  a  pickoff  in  the  drain. 
This  latter  effect  makes  the  empirical  value  of  GQ  at  ReQCrit 
greater  than  the  theoretical  value  when  both  the  theoretical 
and  empirical  curves  are  normalized  with  respect  to  their 
values  in  the  inviscid  range. 

No  significant  difference  in  the  shape  of  the  data 
curves  taken  for  the  three  different  values  of  Re;  is  evi¬ 
dent.  A  more  careful  experiment  might  disclose  differences 
in  the  range  of  Re^  between  3  and  5,  for  there  does  seem  to 
be  slight  difference  in  trend  of  the  three  data  curves  in 
this  region.  This  difference  is  so  small,  however,  that  it 
is  unimportant  with  regard  to  the  immediate  task  of  designing 
a  greatly  improved  rate  sensor. 


5.2 


Prediction  of  G 


It.  is  obvious  from  figures  8  and  10  that  K,  does 
change.  However,  figures  9  and  11  indicate  that  this  is  due 
in  great  part  to  a  change  in  the  Mach  number  at  the  pickoff. 

A  slight  effect  is  present  that  cannot  be  correlated  with  a 
change  in  Mach  number,  because  the  error  in  the  data  points 
does  not  completely  coincide  in  figures  9  and  11.  This  may 
have  been  caused  by  an  inaccuracy  in  (92)  or  by  the  way  in 
which  gain  was  measured  from  the  rough  data.  Gain  was  de¬ 
fined  as  the  slope  in  the  linear  portion  of  the  plot  of  pres¬ 
sure  output  versus  angular  rate.  Since  the  shape  of  the  output 
curve  was  at  times  somewhat  distorted,  and  was  generally  not 
identical  for  two  different  gases  at  the  same  Mach  number,  it 
is  possible  that  the  manner  in  which  the  linear  range  was 
selected  graphically  could  introduce  a  small  error. 

It  is  extremely  unlikely  that  the  lack  of  correla¬ 
tion  in  figures  9  and  11  is  caused  by  a  dependence  of  the 
drain  losses  upon  Re,,  .  In  both  figures,  the  air  data  is  about 
10  percent  lower  than  that  of  the  other  gas  near  M  =  0.9,  for 
example,  even  though  Re,  for  carbon  dioxide  was  about  twice 
that  for  air  in  the  region  and  Ret  for  helium  was  an  order 
of  magnitude  smaller.  Any  dependence  on  Re  would  be  ex¬ 
pected  to  be  monotonic. 

5 . 3  Prediction  of  G^ 

Figures  12  through  15  indicate  that  the  differences 
between  (65)  and  (79)  are  slight  and  that  qualitatively  both 
predict  the  frequency  response.  Figures  16  and  17  indicate 
that  the  error  in  this  data  is  independent  of  Re, .  These 
figures  imply  that  the  viscous  loss  near  the  corner  is  a 
function  of  the  shape  of  the  tangential  velocity  profile  up¬ 
stream  of  the  corner.  In  particular,  the  blunter  the  tan¬ 
gential  profile,  the  greater  the  loss  will  be, as  discussed 
in  section  2.4.3. 

5 . 4  Use  of  Semiempirical  Equations 

Some  very  simple  design  equations  that  are  valid  in 
the  range  Re'  <  1  can  be  formulated  on  the  basis  of  the  dis¬ 
cussion  in  tne  two  previous  sections.  In  section  5.2  it  was 
shown  that  (91)  is  valid  to  within  about  10  percent  if  a 
single  experiment  is  first  conducted  to  establish  the  func¬ 
tional  dependence  of  K3  on  M.  In  section  5.3,  it  was  shown 
that  the  error  in  predicting  is  a  function  only  of  w  in 

this  range  of  Re^.  Again,  only  a  single  experiment  need  to 
be  conducted  to  establish  this  dependence. 
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corner  has  a  large  radius  of  curvature  (to  reduce  viscous 
losses) ,  the  attenuation  due  to  the  corner  will  not  be 
as  large  as  that  indicated  by  figure  16.  This  could  lead 
to  a  slender  stability  margin  in  a  control  system,  since 
the  phase  shift  is  not  accompanied  by  significant  attenua¬ 
tion.  It  would  be  easier  to  design  a  stable  system  using 
a  sensor  with  fully  developed  flow  because,  for  example, 
the  amplitude  in  such  a  design  would  be  down  about  20  dB 
when  the  phase  lag  is  180  degrees. 

5.5  Variation  with  r  of  the  Shape  of  the  Radial 

Profile 

The  change  in  shape  of  the  radial  velocity  profile 
as  the  fluid  flows  toward  the  drain  does  not  appear  to  af¬ 
fect  the  sensor  characteristics  significantly.  The  most 
significant  variation  in  the  two  solutions  obtained  is  in  the 
magnitude  of  the  curves  for  GQ  ,  The  shape  of  the  two  theo¬ 
retical  curves  for  G  agree  quite  closely  as  do  both  the  magni¬ 
tude  and  shape  of  th@  curves  for  G^w  . 

It  appears  on  first  reflection  that  the  suape  of 
the  curve  of  GQ  versus  Re^  obtained  from  the  stepwise  pro¬ 
cedure  need  not  be  similar  to  that  obtained  from  either  hy¬ 
pothetical  solution.  Since  the  ratio  of  the  acceleration 
term  to  the  viscous  term  in  (30)  is  proportional  to  Re^,  as 
indicated  by  (31),  the  acceleration  effects  will  decrease  as 
Re«  decreases  and  the  profile  will  become  more  nearly  para¬ 
bolic  throughout  the  chamber.  The  results  of  the  stepwise 
procedure  might  then  be  expected  to  yield  a  solution  for  G 
that  agrees  more  closely  with  the  solution  for  a  parabolic 
radial  profile  as  Re^  decreases.  However,  this  tendency 
is  opposed  by  another  effect. 


J3 


The  opposing  effect  can  be  shown  by  first  exam¬ 
ining  some  features  of  Uj  and  u2  where  =  0.  As  indicated 
by  (53),  g,  is  proportional  to  Re*.  As  k,  becomes  more 
nearly  uniform,  the  proportionality  constant  decreases 
slightly.  Consequently,  in  the  stepwise  procedure  outlined, 
the  Uj  in  each  annular  region  contributes  a  positive  portion 
to  the  U2  in  the  next  inner  region. 

This  implies  that  if  the  stepwise  solution  were 
known  for  infinite  RQ,  the  magnitude  of  the  tangential 
velocity  relative  to  the  wall  predicted  by  the  stepwise 
solution  would  be  greater  everywhere  than  the  solution  ob¬ 
tained  by  assuming  kj  =  kQ  everywhere.  But  most  importantly, 
the  ratio  of  the  former  to  the  latter  at  Rj  increases  as  the 
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changes  in  k,  upstream  become  more  significant.  Now  since 
the  change  in  the  radial  profile  between  R  and  R,  is  greater 
with  increasing  Re^,  this  implies  that  the  ratio  of  GQ  pre¬ 
dicted  by  the  stepwise  analysis  to  the  G  predicted  by  every¬ 
where  setting  kj  =  k  will  increase  as  Re^  increases.  This 
effect  then  causes  tne  stepwise  curve  to  attempt  to  approach 
the  parabolic  curve  as  Re^  increases,  in  direct  opposition 
to  the  former  effect. 

Since  these  effects  oppose  each  other,  the  shape  of 
the  curve  obtained  from  the  stepwise  procedure  is  probably 
close  to  that  obtained  for  either  of  the  theoretical  cases; 
only  its  magnitude  is  not  well  defined.  This  poses  no  prob¬ 
lem  for  the  use  of  these  results  as  design  equations,  for  it 
contributes  only  to  an  uncertainty  in  K3  in  (91) ,  a  constant 
which  must  be  measured  empirically  anyway. 

5.6  Temperature  Effects 

No  effort  was  made  in  this  experiment  to  investigate 
whether  the  temperature  dependence  of  the  sensor  is  influenced 
if  the  suggested  change  in  viscous  effects  is  pursued.  How¬ 
ever,  this  problem  can  now  be  discussed  from  an  analytical 
viewpoint. 

Since  the  shape  of  the  stepwise  curve  for  G0  will  be 
very  nearly  similar  to  that  of  either  (64)  or  (78),  (61), 

(63)  and  (78)  can  be  used  to  show  that  for  Re^  <7,  the  di¬ 
mensional  gain,  G  ,  is  given  by 


m  (. 

3  =  K4 m  ft 


where  the  product  of  K4  and  Ks  is  a  function  of  h,Rlf  the 
pickoff  design,  and  the  corner  losses.  The  portion  of  the 
right  of  (93)  contained  in  the  brackets  is  just  the  contribu¬ 
tion  due  to  G0 •  Clearly  m/y  goes  as  Reo  if  temperature 
changes.  Then  from  figure  7,  the  change  in  G0  with  tempera¬ 
ture  is  a  minimum  near  Re<}  =  7,  and  a  maximum  at  small  Re£  . 

Now  consider  the  full  expression  for  G.  If  the  change 
in  m  with  respect  to  temperature  has  the  same  sign  as  that  of 
m/y,  then  decreasing  Re^  below  7  increases  the  temperature 
sensitivity  of  G.  If  the  sign  of  the  changes  is  opposite, 
however,  this  can  lead  to  a  reduction  of  temperature  effects, 
and  possibly  even  temperature  insensitivity.  This  point  de¬ 
serves  further  elaboration. 

If  the  temperature  changes  by  some  small  amount  AT, 
then  it  can  be  shown  that  the  change  in  gain,  AG,  is  given  by 
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where  the  subscripts  on  m  and  y  refer  to  the  respective  values 
at  the  initial  temperature. 

If  the  flow  is  inviscid  (that  is,  if  Re£>  is  large),  ' 

then  ^  (95) 

G  m0 


If  ReQ  is' very  small,  (94)  can  be  used  to  show  that 

A_G  „  Am  _ 

G  c  m0  "  ^0 

Since  (96)  describes  the  greatest  possible  contribution  of 
G0  to  the  temperature  effects,  it  can  then  be  sbov’n  from 
(95)  and  (96)  that  temperature  sensitivity  can  be  reduced 
below  that  for  inviscid  flow  if  Am  Au 


Let  us  then  define  the  quantity 


A  m  A/i 
m0  Mo 


N  =  Aii 

MO 


It  can  then  be  shown  from  (97)  that  the  temperature 
sensitivity  increases  as  Reo  is  decreased  if  N<1.  (96)  can 
be  used  to  show  that  if  1<N<2,  the  temperature  sensitivity 
decreases  as  Re^  is  decreased.  If  N>2,  (94)  can  be  used  to 
show  that  gain  is  insensitive  to  small  changes  in  tempera¬ 
ture  at  some  value  of  Re^  .  In  the  fully  developed  range 
this  value  can  be  approximated  by  solving  the  expression 

Re°  (fe"  e  2R*°  “  n-2 

Two  typical  cases  can  be  considered  to  discuss  these 
effects.  Suppose  first  that  a  sensor  using  a  hydraulic  fluid 
is  powered  by  a  constant  pressure  source.  If  the  sensor  de¬ 
sign  is  that  of  a  plenum,  chamber,  and  a  drain  in  series,  then 
the  source  sees  the  sensor  as  a  series  combination  of  an 
orifice  and  a  laminar  resisto1'  The  resistance  of  the  orifice 
does  not  change  appreciably  with  temperature,  but  a  decrease 
in  viscosity  causes  a  decrease  in  laminar  resistance,  and  thus 
an  increase  in  mass  flowrate.  Therefore,  N  is  negative,  and 
temperature  sensitivity  increases  as  Re^  decreases. 

Now  suppose  a  pneumatic  sensor  is  powered  by  a  gas 
bottle  that  supplies  fluid  to  the  sensor  through  a  choked 
orifice.  Since  the  Mach  number  at  this  orifice  is  one  for 
all  temperatures,  the  velocity  through  the  orifice  goes  as  the 
speed  of  sound;  that  is,  proportional  to  the  square  root  of 
temperature.  The  density  in  the  bottle  does  not  change  with 
temperature  if  the  volume  of  the  bottle  can  be  considered  con¬ 
stant,  so  that  the  density  at  the  minimum  cross  section  of 
the  orifice  is  also  independent  of  temperature  if  the  ratio 
of  specific  heats  can  be  considered  constant.  Therefore, 
mass  flowrate  through  the  sensor  increases  in  proportion  to 
the  square  root  of  temperature-  If  y  is  expressed  in  terms 
of  temperature  by  a  relationship  such  as  Sutherland's  for¬ 
mula14,  it  can  be  shown 
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that  for  the  most  commonly  used  gases  operating  in  the  standard 
military  temperature  range,  N  falls  in  the  range  1<N<2. 
Therefore,  the  temperature  sensitivity  of  the  sensor  de¬ 
creases  as  Reg  decreases. 

The  slope  of  m  as  a  function  of  temperature  can  be 
made  to  decrease,  thereby  increasing  N,  if  the  orifice  is 
not  choked,  but  loaded  by  a  laminar  resistance.  In  general, 
it  should  be  possible  to  choose  a  particular  combination  of 
orifice  and  laminar  resistance  for  a  given  mean  operating 
temperature  such  that  N  is  equal  to  2  'for  small  changes  of 
temperature.  At  that  value  of  N,  the  sensor  gain  is  insensi¬ 
tive  to  temperature  for  values  of  Keo  that  are  sufficiently 
small  to  cause  the  contributions  of  U2  to  G0  to  be  negligible. 
If  a  value  of  Reo  less  than  0-.  1  were  chosen,  for  instance, 
the  temperature  could  increase*  without  limit  and  decrease  by 
at  least  two  orders  of  magnitude  without  causing  any  sig¬ 
nificant  contribution  from  U2  to  be  felt.  The  degree  of 
temperature  insensitivity  would  therefore  be  specified  by 
the  degree  to  which  N  could  be  held  constant  over  the  tem¬ 
perature  range. 


5.7  Other  Observations 


It  is  apparent  that  u2  affects  G0  and  G^  at  values 
of  Re'  significantly  less  than  7,  the  value  at  which,  it  was 
argued,  the  radial  velocity  profile  becomes  fully  developed. 
However,  deviations  in  G0  caused  by  the  entrance  effects  on 
the  radial  profile,  as  discussed  in  sections  2.4.2  and  5.1, 
appear  slight.  The  entrance  lengths  for  the  two  components 
of  velocity  thus  appear  different,  but  since  G0  is  not  very 
sensitive  to  the  shape  of  the  radial  profile,  a  more  careful 
analysis  of  the  development  of  the  radial  velocity  profile 
is  needed  to  establish  formally  the  influence  of  the  entrance 
conditions  on  that  profile  for  this  range  of  Re^. 


It  should  be  noted  that  the  chamber  pressure  in  the 
frequency  response  test  for  which  Re^  =  0.806  was  22  psig, 
while  for  Re^  =  0.273,  it  was  4  psig,  causing  a  significant 
change  in  compressible  effects  in  this  experiment.  Since 
these  effects  would  necessarily  be  localized  in  or  very  near 
the  drain,  they  could  be  lumped  with  the  other  drain  effects 
as  a  matter  to  be  ignored  in  this  analysis  in  order  to  es¬ 
tablish  their  significance.  They  clearly  have  no  effect  on 
the  dynamic  characteristics,  as  shewn  by  figures  16  and  17. 
Their  effect  on  dc  gain  has  been  discussed  iii  section  4.2. 


It  is  also  interesting  to  note  that  the  signs  of  the 
contribution  in  G-j^  due  to  U2  is  negative,  while  that  due  to 
ui  is  positive.  It  may  then  be  possible  to  choose  conditions 
in  which  the  two  contributions  are  nearly  180  degrees  out  of 
phase  and  of  equal  amplitude,  so  that  the  combined  signal  is 
very  small.  (65 )•  can  be  used,  for  instance,  to  show  that  for 
Red,  =  1.3,  as  u'  increases  from  15  to  21  the  amplitude  drops 
from  -9db  to  -28  db  and  then  recovers  to  -12db  at  to  *  =  30. 

A  similar  phenomenon  occurs  at  slightly  larger  values  of  Re^ 
and  to'  if  (79)  is  used.  It  would  therefore  appear  that  a 
sensor  with  a  very  simple  design  could  be  made  to  have  a 
notched  frequency  response. 


6.  CONCLUSIONS  ABOUT  PRESENT  MODEL 

The  nature  of  the  gain  and  frequency  response  of 
the  present  rate  sensor  design  in  the  fully  developed  range 
can  be  closely  predicted  through  the  use  of  this  chamber 
analysis  and  a  single  test.  Tests  carried  out  in  the  range 
of  Re^  <  1  demonstrate  that  the  error  in  this  analysis  is 
not  a  function  of  Re! ,  which  implies  that  losses  at  the  cor¬ 
ner  and  in  the  drain  are  independent  of  Re, .  This  also  in¬ 
dicates  that  once  a  single  test  is  made  on  the  present  sen¬ 
sor  to  specify  those,  losses,  the  results  of  making  any  de¬ 
sign  change  other  than  one  in  the  pickoff-drain  configura¬ 
tion  can  be  predicted  to  within  at  least  10  percent. 

The  change  in  shape  of  the  radial  velocity  profile 
as  the  fluid  flows  toward  the  drain  does  not  significantly 
affect  the  sensor  characteristics.  The  only  significant  un¬ 
certainty  created  by  the  approximate  method  used  in  this  an¬ 
alysis  appears  in  the  magnitude  of  the  sensor  gain  at  Re£ 

This  gain  must  be  measured  anyway  in  the  semi-empirical  de¬ 
sign  equations  presented  because  of  the  incomplete  analysis 
of-  the  pickoff-drain  region.  There  is  some  uncertainty  in 
the  prediction  of  at  values  of  w'  near  i5  and  Rep  near  1, 
since  the  solutions  do  not  provide  bounds  in  those  ranges  of 
the  governing  parameters.  However,  in  the  present  design 
w:-th  the  pickoff  in  the  drain,  the  attenuation  at  the  corner 
is  ro  great  that  this  .uncertainty  is  not  of  great  consequence. 

In  summary,  from  a  control  systems  point  of  view 
the  dynamic  characteristics  of  a  rate  sensor  operating  in 
the  fully  developed  range  are  significantly  more  favorable 
than  those  of  an  inviscid  sensor.  Not  only  can  the  phase 
lag  in  the  output  be  much  less,  but  this  lag  is  accompanied 
by  an  attenuation  in  amplitude. 

If  only  the  height  of  the  present  design  is  changed  to 
induce  the  viscous  effects,  the  gain  is  reduced  and. the  temper¬ 
ature  sei.sitivity  may  increase  or  decrease,  depending  on  the 


initial  design.  However,  other  new  design  alternatives  can 
then  be  used  which  appear  beforehand  to  allow  improvement  of 
sensor  performance  in  all  areas,  as  will  be  discussed  next. 


7. 


FUTURE  MODELS 


7.1 


Experimental  Work 


7.1.1  Chamber  Height 

Sensors  can  currently  be  fabricated  with  chamber 
heights  as  small  as  several  thousandths  of  an  inch.  Anal¬ 
ysis  indicates  these  sensors  should  have  corner  frequen¬ 
cies  in  excess  of  1  kHz.  Models  are  being  fabricated  to 
test  this  prediction. 


7.1.2  Pickoff  Location 


1 


In  the  future  models,  the  signal -pickoff  will  be  lo¬ 
cated  in  the  chamber  just  upstream  of  the  drain  with  its  longi¬ 
tudinal  axis  parallel  to  the  drain  axis.  This  change  has  a 
wide  impact  on  sensor  performance.  Some  of  these  effects  will 
be  discussed. 

7 . 1 . 2 . 1  Noise  Reduction 


The  effect  of  drain  turbulence  on  the  output  can 
be  reduced  if  the  pickoff  is  located  in  the  chamber  upstream 
of  the  disturbance .  In  addition,  the  problem  of  designing 
a  drain  in  which  the  flow  is  stable  can  then  be  isolated 
from  problems  relating  to  the  pickoff  fabrication  and  mounting. 

If  the  flow  separates  on  the  surface  of  a  probe  in 
the  chamber  even  though  the  undisturbed  field  is  accelerating, 
the  problem  of  pickoff  turbulence  can  be  attacked.  This  tur¬ 
bulence  could  be  suppressed  with  the  top  and  bottom  walls  if 
the  chamber  is  sufficiently  thin  at  the  pickoff  location, 

A  pickoff  in  the  chamber  could  be  bonded  to  the  top 
and  bottom  walls,  so  that  a  pickoff-drain  insert  could  be 
fabricated  that  does  not  respond  to  thermal  or  mechanical 
disturbances.  The  signal  null  must  then  be  adjusted  by  some 
mechanism  external  to  the  pickoff.  If  this  can  be  achieved 
the  problem  of  adjusting  the  null  will  be  isolated  from  prob¬ 
lems  concerning  the  pickoff  design. 
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It  then  appears  that  the  problem  of  pickof f-drain 
noise  and  null  shift  could  be  attacked  systematically  if  the 
pickof f  is  located  in  the  chamber. 


7. 1.2. 2  Frequency  Response 


Relocating  the  pickoff  upstream  of  the  corner  would 
eliminate  corner  losses.  Figures  16  through  19  indicate 
that  frequency  response  would  be  about  doubled  for  (J  of  the 
order  of  15  or  less.  This  may  make  the  use  of  the  fully 
developed  sensor  in  a  control  system  even  more  attractive. 


It  can  be  shown  (App.  I)  that  if  there  are  no 
corner  losses,  for  small  Re^  and  large  oS  the  phase  lag  ap¬ 
proaches  a  limit,  <j>-max-,  ^where  135‘'<({>max  <225° .  This  maximum 
could  probably  be  reduced  by  changing  the:  pickoff  design  to' 
respond  :r.ore  strongly  to  the  information  near  the  wall. 


7. 1.2. 3  Signal  Gain 


The  relocation  of  the  pickoff  in  the  chamber  just 
upstream  of  the  drain  will  result  in  a  change  of  gain.  Equa¬ 
tion  (61)  indicates  that,  for  a  pickoff  in  the  chamber,  gain 
is  a  linear  function  of  pVr  and  Vg  in  the  undisturbed  flow 

field  at  the  pickof f's  longitudinal  axis,  as  it  is  of  pV2  ana 
Vq  for  a  pickoff  in  the  drain.  The  ratio  of  pvr  (at  the  cham¬ 
ber  location)  to  pVz  (at  the  drain  location)  might  realistically 
be  slightly  more  than  one-half  if  the  average  transport  veloc¬ 
ity  component  is  nearly  constant  in  the  vicinity  of  the  cor¬ 
ner,  as  it  is  in  the  typical  design.  The  ratio  of  the  tan¬ 
gential  velocities  might  also  be  about  one-half.  The  significant 
viscous  effect  at  the  center  of  the  drain  account  for  the 
difference  between  this  value  and  a  smaller  value  that  would 
be  predicted  by  inviscid  analysis.  The  combined  effect  is  a 
gain  reduction  by  about  4  for  a  typical  design. 


This  initial  reduction  of  gain  must  he  considered  a 
loss  in  sensor  effectiveness,  even  when  the  increased  fre¬ 
quency  response  is  considered.  However,  this  design  change 
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sets  the  stage  for  other  changes  which  can  result  in  an  in¬ 
creased  effectiveness,  the  first  step  of  which  is  to  reduce 
chamber  height.  The  effect  that  this  latter  change  has  on 
signal  gain  will  be  discussed  first.  For  sake  of  a  clear 
comparison,  the  gain  will  be  expressed  in  terms  of  that  of 
a  typical  present  sensor  operating  at  the  lower  limit  of  the 
inviscid  range  with  a  pickoff  in  the  drain. 

The  gain  before  and  after  the  relocation  is  indi¬ 
cated  by  points  A  and  B  respectively  in  figure  21.  Now  let 
h  be  reduced  while  the  rest  of  the  geometry  is  held  fixed. 

The  restrictions  on  power  supply  must  now  be  considered  in 
order  to  relate  A  to  h‘. 

Suppose  first  that  A  is  held  fixed  while  _h  is  re-^_ 
duced  with  the  pickoff  in  the  chamber ;  If  inviscid  gain  for 
these  conditions  is  multiplied  by  G0  from  the  analysis,  it 
can  be  shown  that  actual  gain  increases,  approaches  a  limit, 
and  is  nearly  constant  for  Re£<l,  as  indicated  in  figure  2(1. 
In  the  latter  "range  the  gain  is  nearly  equal  to  that  obtained 
from  a  pickoff  in  the  drain  with  Re£>  =  7. 


m= CONST 
^supply  =  CONST 


Reo  =  ? 


HEIGHT 


Figure  2i.  Effect  on  Gain  of  Relocating  Pickoff  and  Reducing  h 
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Consider  next  that  the  supply  pressure  is  held  con¬ 
stant  while  h  is  reduced  with  the  pickoff  in  the  next  chamber. 
If  one  assumes  that  the  curtain  area  (that  is,  2-rrRih)  is  the 
dominant  restriction,  then  it  can  be  shown  that  gain  is  nearly 
independent  of  h  for  Re£  nearly  7  and  proportional  to  h2  for 
Re£  less  than  1.  The  curve  representing  any  other  type  of 
power  supply  or  system  specification  on  consumption  must  lie 
between  these  two  curves. 

The  curves  for  the  same  two  cases  with  the  pickoff 
in  the  drain  are  represented  in  this  same  figure  by  the  curves 
that  meet  point  A.  It  appears  that  if  h  is  made  very  small 
on  a  model  with  a  pickoff  in  the  chamber,  the  gain  will  bo 
reduced  below  point  A  for  any  power  supply  that  responds  to 
downstream  pressure.  However,  this  reduction  is  less  than 
would  have  occurred  if  the  pickoff  had  remained  in  the  sink. 

In  addition,  this  sacrifice  in  gain  has  resulted  in  an  appre¬ 
ciable  increase  in  frequency  response,  as  will  be  discussed 
next. 

7 . 1 . 2 . 4  '  Gain-Bandwidth  Product 


A  quantity  which  is  often  of  interest  to  the  de¬ 
signer  of  a  control  system  is  the  gain-bandwidth  product  (or 
sensitivity-frequency  response  product,  since  this  is  a  sen¬ 
sor)  .  The  gain-bandwidth  product  before  and  after  the  trans¬ 
fer  of  the  pickoff  are  presented  in  figure  22  by  points  A  and 
B,  respectively.  For  the  present  typical  design,  the  reduc¬ 
tion  might  well  be  a  factor  of  about  two. 


Figure  22.'  Effect  on  Gain-bandwidth  Product  of  Relocating 

Pickoff  and  Reducing  h. 


69 


Consider  first  the  case  in  which  m  is  held  fixed 
while  h  is  reduced.  It  can  be  shown  that  frequency  response 
increases  by  about  a  factor  while  Re^  is  reduced  from  7  to 

1  and  is  proportional  to  1/ph  for  Re^  <  1.  Figure  21  can 
then  be  used  together  with  this  information  to  show  that  the 
gain-bandwidth  product  at  Re^  =  1  is  about  30  times  the  value 
at  point  A.  For  Re'  <1,  it  increases  as  1/ph2. 
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For  the  case  when  supply  pressure  is  held  constant, 
the  gain-bandwidth  product  increases  and  approaches  a  limit 
as  h  is  reduced.  It  can  be  shown  that  the  value  for  Re*  <  i 

Q 

is  about  three  times  the  value  at  point  A. 

Clearly,  then,  for  any  type  of  power  supply  condi¬ 
tions  the  gain-bandwidth  product  can  be  increased  signifi¬ 
cantly  over  the  value  in  the  present  inviscid  sensor  by 
placing  the  pickoff  in  the  chamber  and  reducing  the  height. 

For  small  h,  the  gain-bandwidth  product  is  greater  for  a  pick¬ 
off  in  the  chamber  than  for  one  in  the  drain  regardless  of  the 
power  supply  used.  The  improvement  is  especially  large  if  the 
mass  flowrate  is  not  a  strong  function  of  the  load  that  the 
curtain  area  at  the  drain  entrance  represents. 

7. 1.2. 5  Airfoil  Gain  and  Output  Impedance 


Although  the  transfer  of  the  pickoff  initially  re¬ 
sults  in  a  loss  of  gain  if  nothing  else  is  changed,  changes 
in  the  pickoff  design  that  may  improve  performance  can  then 
be  made. 
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Since  a  strong  faborable  pressure  gradient  exists 
in  the  flow  field  near  the  sink  due  to  acceleration,  pickoff 
designs  are  not  as  severely  restricted  by  flow  separation. 

It  may  be  possible  to  use  a  pickoff  then  whose  gain  is  higher 
than  those  presently  used;  perhaps  a  pickoff  design  that 
takes  an  appreciable  amount  of  flow  through  the  sensing  ports 
could  be  used  without  destroying  the  nearby  pressure  field. 

This  would  increase  the  gain  and  decrease  the  output 
impedance . 

7.1.3  Supply  Distribution  System 

A  model  is  being  designed  at  HDL  that  has  a  flow 
distribution  system  upstream  of  the  working  chamber  that  should 
reduce  the  magnitude  of  upstream  nonaxisymmetric  disturbances 
significantly  more  than  the  present  design.  This  should  allow 
the  influence  on  the  sensor  output  noise  level  of  disturbances 
near  the  drain  and  changes  in  the  mean  power  level  to  be 
studied. 
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7.1.4 


Multiple  Pickoffs 


It  can  be  shown  from  the  solution  for  V.obtained  in 

U 

this  analysis  that  for  Re^  <  1  the  lines  connecting  points 
near  the  pickoff  that  have  the  same  phase  with  respect  to  the 
wall  motion  are  nearly  parallel  to  the  wall  or,  in  other 
words,  that  the  signal  propagates  into  the  fluid  in  a  direc¬ 
tion  normal  to  the  walls.  Since  the  noise  propagate-  in  a 
direction  normal  to  that  of  the  signal,  then  the  time  his¬ 
tory  of  the  wall  motion  could  possibly  be  traced  independently 
of  that  of  the  noise. 

It  can  be  shown  that  the  output  of  n  pickoffs,  each 
of  which  measures  the  signal  at  a  different  elevation  in  the 
chamber,  can  measure  by  simple  arithmetical  operations  the 
first  through  the  nth  time  derivative  of  angular  motion.  In 
addition,  the  shape  of  the  frequency  response  curve  for  q 
can  be  shaped.  The  output  of  two  pickoffs  could  be  used, 
for  instance,  to  measure  both  6  and  #,  with  the  amplitude  of 
the  angular  rate  signal  relative  to  d-c  made  arbitrarily 
small  at  a  phase  lag  of  about  120  degrees. 

Some  careful  design  may  be  required  to  extract  the 
information  at  different  levels  in  the  fluid.  Placing  in  the 
fluid  a  multi-channeled  probe  that  has  a  number  of  pickoff 
ports  at  various  points  along  its  length  would  probably  cause 
such  a  severe  local  disturbance  that  the  variation  of  infor¬ 
mation  across  the  channel  is  washed  out,  an 2  the  output  from 
any  port  along  the  probe  may  be  nearly  the  same.  However, 
there  may  be  some  other  method  a  clever  designer  might  pro¬ 
pose  that  can  extract  a  flueric  signal.  In  addition,  other 
methods  such  as  anemometry  provide  ready  access  to  this  in¬ 
formation  if  the  systems  specifications  allow  an  immediate 
interface  to  electronics. 


7.1.5 


Chamber  Design 


Present  plans  also  call  for  the  testing  of  an  axi- 
symmetric  model  in  which  the  chamber  is  something  other  than 
cylindrical.  There  are  several  objectives  to  this  effort, 
one  of  which  is  to  design  a  model  whose  gain  does  not  change 
for  small  excursions  of  the  environmental  temperature  or 
power  supply  level,  or  some  combination  of  these  two  param¬ 
eters.  The  finer  details  of  this  design  will  hopefully  be 
specified  by  a  theoretical  program  presently  in  progress. 

7.2  Theoretical  Studies 

It  can  presently  be  shown  that  the  approximate 
governing  equations  used  in  this  analysis  can  also  be  solved 


for  a  thin  chamber  in  which  the  height  varies  in  any  mono¬ 
tonic  fashion  with  r.  This  solution  will  be  pursued. 

It  can  also  be  shown  that  for  the  cylindrical  cham¬ 
ber,  the  growth  or  decay  with  radius  of  perturbations  with 
respect  to  z,  B  and  t  can  be  treated  analytically.  In  ad¬ 
dition,  stability  in  the  supply  distribution  system  and  in 
the  drain  is  being  studied. 
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10.  LIST  OF  SYMBOLS 


=  3Q0/4irh 

=  Q  /2irh 

o 

=  \/\<Ji/V 


V  V  F: 


=  ratio  of  terms  in  two  series 

=  function  that  describes  the  propagation  of 
information  from  the  coupler 

=  function  that  describes  the  propagation  of 
information  upstream  from  the  drain 

—  eigenfunction  in  the  drain 

=  body  forces 
=  d-c  gain 

=  d-c  gain  relative  to  inviscid  case 
=  a-c  gain  relative  to  d-c  level 
=  chamber  height 


k  (r,  z) 


=  normalized  radial  velocity  profile 


=  entrance  length 

=  mass  flow  rate  through  the  sensor 
=  pressure 

=  volumetric  flow  rate  in  the  outer  portion 
of  the  chamber 

=  radial  distance 

=  chamber  radius 

=  drain  radius 

=  chamber  Reynolds  number  =  U  R  /v 

o  o 

=  modified  Reynolds  number  =  (U  R  /p)  (hJ/RQ2) 
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ZdrainlaVg> 


|Ri/j> 


=  nth  eigenvalue 


-  transport  time 

=  function  describing  the  propagation  of 
information  from  the  top  and  bottom 
plates 

=  function  describing  the  propagation  of 
information  from  the  coupler  and  drain 
into  the  chamber 

=  average  radial  velocity  at  Rq 

=  free  stream  velocity 

=  radial  velocity 

*  axial  velocity 


0rel 


=  tangential  velocity 

=  tangential  velocity  relative  to  the  wall 
=  distance  from  leading  edge 


=  distance  from  the  centerplane  of  the 
chamber 

=  angle  between  a  streamline  and  radial 
line 

=  nth  eigenvalue 

=  sink  tube  radius 
=  boundary  layer  thickness 
=  displacement  thickness 
=  nth  eigenvalue 
=  momentum  thickness 
=  nth  eigenvalue 
=  shape  factor 


dynamic  viscosity 
kinematic  viscosity 
density 

output  of  rate  sensor 

eigenfunction  m  the  chamber 

function  that  describes  propagation 
ot  information  down  the  drain 

oscillatory  frequency  of  rate  sensor 
2 

ail  /v 

angular  rate  amplitude 


APPENDIX  A.  VISCOUS  CORE  OF  A  2-D  VORTEX-SINK  FLOW  FIELD 


The  existence  of  a  rigid  body  core  in  two-dimension¬ 
al  vortex-sink  flow  is  usually  not  questioned.  Two  arguments 
imply  its  existence.  The  first  involves  discussion  of  the 
isentropic  energy  equation,  in  which  the  existence  of  a  limit 
circle  is  described.  The  second  is  a  description  of  an  ob¬ 
served  Rankine  (or  combined  forced  and  free)  vortex . 1 6 


The  first  discussion  indicates  that  energy  limita¬ 
tions  exist  in  a  potential  field.  If  a  set  of  stagnation 
conditions  is  established  and  the  strength  of  the  sink  and 
vortex  is  specified,  this  approach  shows  that  a  radius  exists 
at  which  the  absolute  pressure  is  zero.  The  region  inside 
this  limiting  radius  is  therefore  inaccessible  to  the  fluid 
'and  is  often  assigned  a  rigid  body  nature. 


What  this  approach  does  point  out  is  that  if  the 
outlet  is  placed  inside  the  limit  circle,  no  such  flow  field 
can  exist.  If  the  specified  features  of  the  external  flow 
are  somehow  established  with  the  outlet  placed  inside  the 
limit  circle,  the  flow  rate  cannot  be  maintained  and  there¬ 
fore  is  reduced.  The  field  eventually  adjusts  to  a  new 
sink  strength  as  defined  by  the  pressure  specified  in  the 
outlet  and  the  stagnation  conditions  in  the  external  field. 
The  fact  that  the  radius  of  the  limit  circle  changes  if  a 
swirl  is  introduced  in  the  external  field  indicates  that 
the  imposed  vortex  throttles  the  flow,  a  well-known 
phenomenon . 


The  nature  of  the  flow  field  in  the  interior  of 
the  limit  circle  is  then  at  best  academic  because  such  a 
field  cannot  exist  in  reality.  The  fact  that  it  is  often 
assigned  a  rigid-body  nature  may  lead  some  readers  to  be¬ 
lieve  that  a  rigid-body  core  does  exist  near  the  center  of 
a  two-dimensional  sink  flow  field. 


The  discussion  of  a  Rankine  vortex,  often  accom¬ 
panied  by  photographs  of  the  region  above  the  outlet  of  a 
hydraulic  sink -vortex  flow  field, may  lead  to  confusion  if 
it  is  not  made  quite  clear  that  this  rigid-body  core  is 
necessarily  a  feature  of  the  three-dimensional  nature  of 
the  flow  field  over  the  outlet. 


For  references,  see  pp.  72-73. 


Quite  often,  the  two-dimensional  vortex-sink  field 
is  discussed  without  any  apparent  concern  about  the  manner 
in  which  the  fluid  exists  r+-  the  center  of  the  field.  The 
popular  description1 7  is  tir:  the  field  is  a  combination  of 
free  vortex  in  the  outer  region,  a  solid  body  core  in  the 
center,  and  a  transition  region  between.  If  this  descrip¬ 
tion  is  applied  to  a  zruly  two-dimensional  vortex-sink  field, 
then  it  poses  several  physical  and  mathematical  absurdities . 
Instead  of  discussing  these,  however,  it  might  be  test  to 
attack  the  problem  directly. 

In  order  to  understand  the  two-dimensional  vortex- 
sink  field  more  clearly,  consider  the  field  shown  in  figure 
A-l.  The  physical  boundaries  consist  of  two  very  long  con¬ 
centric  cylinders  that  have  porous  walls.  The  outer  one 
turns  at  a  rate  U0;  the  inner  one  at  a  rate  ft2.  Fluid  enters 
the  region  between  the  two  in  an  axi symmetric  fashion  from 
the  outer  wall  and  leaves  through  the  inner  wall.  The  tan¬ 
gential  velocity  relative  to  either  bounding  wall  is  zero  at 
the  wall.  A  purely  two-dimensional  vortexrsink  field  can  be 
described  by  first  finding  the  solution  to* the  incompressible 
momentum  equations  for  the  region  between  these  two 
cylinders  and  by  then  examining  the  features  of  this  solution 
as  the  radius  of  the  inner  cylinder  is  allowed  to  approach 
zero.  The  implications  of  this  inner  radius  being  infini¬ 
tesimal  (in  relation  .to  the  previous  discussion  of  the  limit 
circle) ,  as  well  as  the  consequences  of  the  assumption  of 
incompressibility  will  be  discussed  after  the  solution  is 
found  so  that  the  ramifications  will  be  more  apparent. 

Equations  (17),  (18),  (19),  and  (20)  become 


Figure  A-l.  Two-dimensional  Vortex-sink  Field. 
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From  (A-l) ,  (A-3) ,  and  (A-4) ,  it  cam  be  shown  that 
the  pressure  gradient  is  given  by 


SM-W 


Since  the  flow  field  is  axisyimnetric  and  also  in¬ 
dependent  of  z,  the  partial  derivative  in  (A-4)  becomes  a 
total  derivative .  Therefore 


Vr-z!l 


f=2 tr  r  ' 


where  q  is  the  volumetric  flow  rate  per  unit  length  of  z. 

Equations  (A-2)  and  (?-6)  can  then-  be  combined  to 


obtain 
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where 


This  equation  has  the  solution 


v„  =  A  +  d 

6  r  rRe-l 

(Re^H) 
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(Re=2). 

If  the  boundary  conditions 
Vg(^o)=  R0a0 


ve  (R2)  =  r2£12 
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then  from  (A-15) 


R2 


But  since 


,im- (lif)1*  =  * 
€—0 


(A-15)  also  shows  that  R* 

-5i  =  6a»Re=2 
r2 


Thus,  (A-15)  can  be  used  for  all  Re. 


This  relationship  can  be  used  to  show  that  R3/R2  is 
bounded  fcr  all  Re>l.  In  particular. 


lim  /R3  N  _  , 
Re~®  Vr7 /  ' 


because 
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The  core  size  thus  approaches  that  of  the  sink  tube  as  Re 
becomes  large.  In  the  experiment  that  is  the  subject  of 
this  report.  Re  is  of  the  order  of  104. 


If  R*  is  reduced,  then  Re  will  change  unless  specific 
supply  and  sink  conditions  are  maintained.  If  the  pressure 
drop  between  the  two  cylinders  is  maintained,  for  instance,  the 
flow  rate  will  gradually  decrease.  (A-6)  and  (A-12)  with 
ft2=0  cam  be  used  in  (A-5)  to  show  that  at  some  R2  the  flow 
stops  if  the  pressure  drop  ii=  less  than  r02Q02  .  If  R2  is  de¬ 
creased,  the  flow  is  then  outward.  .  2 


If  R:  is  made  to  approach  zero  while  Re  is  held  con¬ 
stant,  (A-15)  shows  that  R3  also  approaches  zero  if  Re>l.  If 
Re£l  ,  SVg/Sr  is  not  zero  at  any  point  between  the  two  cylin¬ 
ders.  Thus,  as  R2  approaches  zero,  the  viscous  core  either 
gradually  shrinks  to  zero  or  continues  to  fill  the  entire 
flow  region.  The  behavior  of  the  fluid  in  the  core  is  de¬ 
finitely  not  that  of  a  rigid-body.  Notice  that  the  only  case  in 
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vhion  a  term  linear  in  r  does  appear  in  (A-12)  is  for  the 
case  Re  =  0,  or  when  the  flow  rate  is  zero.  Even  in  that 
case,  though,  the  1/r  dependence  makes  its  strongest  con¬ 
tribution  near  the  center  of  the  vortex,  and  the  rigid-body 
behavior  is  strongest  in  the  outer  portion  of  the  region, 
just  the  opposite  of  a  free  vortex  with  a  rigid-body  core. 

For  given  values  of  R0 ,  R2 ,  Re  and  ,  a  value  for 
ft2  exists  such  that  the  fluid  very  near  the  sink  wall 
behaves  as  a  rigid-body.  However,  this  value  increases 
without  limit  as  R2  is  decreased.  Thus,  it  is  not  possible 
to  construct*  a  situation  in  which  the  fluid  near  the  center 
approaches  a  specific  rigid-body  field  as  R2  is  decreased. 

It  was  specified  in  the  previous  discussion  that 
Re  be  constant.  This  is  equivalent  to  assuming  that  the 
strength  of  the  sink  and  that  the  stagnation  conditions  re¬ 
main  constant,  the  requirements  which  should  be  made  if  one 
were  to  consider  only  the  effect  of  shrinking  the  sink  tube. 
Re  can  be  maintained  constant  simply  by  adjusting  the 
pressure  difference  between  R0  and  R2  (in  keeping  with  the 
solution  of  (A-5))  as  R2  is  changed  so  that  the  mass  flow 
rate  through  the  region  remains  constant. 


If  the  pressure  ratio  between  R0  and  R2  varies  sig¬ 
nificantly  from  1  while  increasing  this  pressure  difference 
and  decreasing  R2 ,  corpressible  effects  must  be  considered. 
The  continuity  equation  will  require 

£(*)♦£■»  <A'17> 

The  term  (pVr)  on  the  left  side  of  (A-2)  remains  unchanged, 
because  (A-17 )  requires 


pVf-- 


constont 

r 


the  same  result  which  is  obtained  if  p  is  constant.  However, 
two  additional  effects  must  be  considered  in  the  viscous  term 
on  the  right  of  (A-2) .  In  deriving  the  incompressible  equa¬ 
tions,  it  is  assumed  that  y  is  constant  and  that  div  V  =  0. 
The  continuity  equation  requires  that  div  (pV)  =  0,  so  that 
div  V  is  not  zero  if  p  is  net  constant.  Also,  y  changes  with 
radius  due  to  a  change  in  temperature  induced  by  compressible 
effects.  However,  since  the  flow  can  at  best  become  sonic 
at  R2  if  nc  is  small,  these  two  effects  will  not  be  severe 
and  thus  will  not  drastically  alter  the  results.  One  should 
expect  the  previous  results  to  predict  the  qualitative  be¬ 
havior  of  the  viscous  core  even  in  this  extreme  case. 
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APPENDIX  B.  ORTHOGONALITY  OF  EIGENFUNCTIONS 


i 

i 

! 


The  complex  space  composed  of  functions  whose  square 
is  integrable  in  the  sense  of  Lebesgue  over  the  range  a<z<b 
becomes  a  Hilbert  space  if  the  scalar  product  of  two  elements 
u{z)  and  v(z)  is  defined18  as 

(u,v)  =«£“v*d2  (B-l) 


where  the  star  indicates  a  complex  conjugate.  This  implies 
that  any  function  that  is  Lebesgue  integrable  can  be  repre¬ 
sented  by  a  convergent  sequence  of  eigenfunctions  generated 
by  a  complex  differential  operator  L,  of  the  type  appearing 
in  (56).  In  order  to  find  the  coefficients  of  such  a'se-' 
guence,  one  needs  to  find  the  eigenfunctions  of  the  adjoint 
operator,  L+,  which  must  satisfy  the  condition1 9 

^um  j  Lvn^  -  ^L^um,vn^  (B— 2) 

where  t^eE,  the  set  of  solutions  of  L+,  and  vneE,  the  set  of 
solutions  of  L. 


It  will  be  shown  that  if  (B-l)  is  used,  (B-2)  is 
satisfied  if  we  choose 

Lt=  L*  (B-3) 

Consider  first  the  operator  of  (56) 

LVnS(^-"rLjVn=^nk.{z)vn  (b-4) 

The  conjugate  equation  is  then 


L*u 


=  ^  k,(z)u 


m 


(B-5) 


Clearly,  the  two  sets  are. related  pairwise  such  that 


um  =vm 


(B-6) 


where  both  sets  have  Oc.?n  indexed  by  the  same  rule.  The  eigen¬ 
values  ar-  also  related  »n  the  same  way;  that  is. 


ft'm  =  Pm 


(B-7) 


If  (B-4)  is  multiplied  by  vm  and  this  result  is  sub¬ 
tracted  from  an  equation  that  is  identical  except  that  the 
indices  are  exchanged,  the  resUi  ; 

vm  -  vn — -j— ES-  =  (£n-A>m)k|(z)vn  vm  (B~8) 

dzz 
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This  equation  can  be  integrated  once  to  obtain 
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becomes 


If  (B-l)  and  (E-4)  through  (B-7)  are  used,  this 


v  ii^iL.v  li u* 

m  dz2  n  dz2 


1  =  (“m.l-Vn)-(L*um,vn) 

Jo 


(B-10) 


It  should  be  noted  that  taking  the  complex  conjugate  of  an 
integral  like  that  in  B-9  is  equivalent  to  reversing  the  order 
of  the  terms  in  the  bracket  in  B-l. 

In  that  class  of  boundary  conditions  for  which  the 
left-hand  side  of  B-10  vanishes  (in  particular  for  homoge¬ 
neous  boundary  conditions,  as  in  the  case  being  considered  in 
this  report) ,  (B-2)  is  satisfied  for  the  relationship  chosen 
in  (B-3) . 

Specifically,  then,  the  complex  conjugate  of  the 
differential  operator  of  (55)  is  the  adjoint  operator,  each 
element  of  E  (the  set  of  solutions  of  L)  is  the  complex  con¬ 
jugate  of  one  element  in  E,  and  each  element  of  E  is  orthogo¬ 
nal  (that  is,  the  scalar  product  is  zero)  to  every  element 
of  E  except  its  conjugate  and  is  not  orthogonal  to  any  ele¬ 
ment  of  E  unless  w  =  0.  When  w  =  0,  the  sets  E  and'E  coincide 
and  L+  =  L. 

A  function  f(z)  that  is  square  integrable  on  the 
interval  (a,b)  can  then  be  represented  by  a  sequence 


Hz)  =  2Cn  vn 


(B-ll) 


where 


r  v0 


Ia  f(z)k,(2  )vndz 

*'0vnk|l2)vndz 


The  function  f(z)  could  also  have  been  represented 
by  a  sequence  in  which  all  the  v's  are  replaced  by  u's.  Al¬ 
though  the  coefficients,  of  these  two  biorthogonal  sequences 
are  not  identical  term  by  term,  the  complete  series  are. 

This  duplicity  is  a  direct  result  of  the  fact  that  L  is  not 
self-adjoint,  that  is,  L+  ^  L. 

Note  ai~o  that  no  solution  of  the  adjoint  equation 
appears  in  (B-ll).  Although  one  must  pair  the  elements  of  E 
and  E  to  generate  an  orthogonal  series ,  the.  conditions  that 
1)  these  sets  are  conjugates  of  each  other  and  that  2)  the 
complex  conjugate  of  one  solution  is  used  in  calculating  the 
scalar  product  combine  to  cause  solutions  from  only  one 
set  to  be  used  in  the  mechanics  of  determining  the  coeffici¬ 
ents  of  the  orthogonal  series.  This  is  of  more  than  passing 
interest,  since  only  one  equation  need  be  solved.  No  diffi¬ 
culty  arises  as  it  does  in  matrix  problems,  where  one  occas¬ 
ionally  must  operate  on  the  left  and  the  right  of  a  quantity 
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using  matrices  that  have  different  eigenvalues,  because  the 
definition  of  scalar  product  used  her®  is  not  analogous  to 
the  definition  of  dot  product  used  in  those  matrix  problems. 


APPENDIX  C. 


RELATIONSHIP  BETWEEN  THE  TWO  HYPOTHETICAL  SOLU¬ 
TIONS  AND  THE  STEPWISE  SOLUTION 


What  will  be  discussed  here  are  the  ranges  of  the 
governing  parameters  for  which  the  two  hypothetical  solutions 
can  be  considered  as  bounding  expressions  for  the  stepwise 
solution.  Before  this  is  undertaken,  it  is  necessary  first 
to  discuss  some  features  of  the  profile  ki . 

Consider  the  flow  in  the  annular  region  about  a 
typical  ri  used  in  the  stepwise  solution.  The  ki  valid 
there  is  a  solution  of  (30).  The  nature  of  this  equation 
is  such  that  the  value  of  this  particular  ki  in  the  range 
o<|z|<h/2  is  greater  than  that  of  the  solution  for  any  r.j 
where  ri'>ri  and  less  than  the  solution  for  any  ri*' where 
ri"<  n .  Furthermore,  the  shape  of  this  solution  changes 
continuously  if  n  is  varied  continuously.  Therefore,  the 
value  of  k\  at  a  particular  z  varies  continuously  ard  in  a 
monotonic  fashion  as  r  is  varied  over  the  full  extent  of 
its  range. 

Now  consider  the  annulus  that  is  bounded  on  the 
upstream  side  by  the  coupler.  Suppose  that  the  pickoff  were 
located  at  the  downstream  limit  of  this  annulus  and  further 
that  it  can  be  shown  that  the  output  from  the  pickoff  varies 
in  a  monotonic  fashion  as  the  profile  is  varied  in  any  arbi¬ 
trary  but  monotonic  fashion  from  parabolic  to  uniform.  One 
solution  would  then  provide  an  upper  bound  for  tne  output 
at  the  downstream  limit  of  the  annulus,  and  the  other  a  lower 
bound. 

If  one  were  to  proceed  with  the  step-wise  solution, 
these  solutions  could  now  be  used  to  provide  an  upper  and  a 
lower  bound  for  the  upstream  boundary  condition  on  the  next 
inner  annulus.  For  instance,  the  upper  bound  for  sensor  out¬ 
put  if  the  pickoff  were  located  »*t  the  downstream  limit  of 
this  second  annulus  could  be  obtained  if  the  solution  which 
provided  the  upper  bound  in  the  previous  annulus  were  con¬ 
tinued  across  this  second  annulus.  To  understand  this,  re¬ 
member  that  a  "maximum"  function  provides  an  upper  bound  for 
the  solution  at  the  downstream  limit  of  a  given  annulus  for 
a  given  upstream  boundary  condition,  and  that  a  "minimum" 
function  gives  a  lower  bound.  If  this  lower  bound  is  used 
in  the  next  inner  annulus  as  the  upstream  boundary  condition 
for  the  minimum  function,  the  result  obtained  is  thus  the 
lower  bound  for  the  downstream  limit  of  this  second  annulus. 
Likewise,  the  upper  bound  at  the  downstream  limit  of  the 
first  annulus  can  be  used  with  the  maximum  function  to  ob¬ 
tain  the  up]-er  bound  at  the  lower  limit  of  this  second  an¬ 
nulus.  This  procedure  is  equivalent  to  extending  the  maxi¬ 
mum  and  minimum  functions  across  both  annuli  without 
bothering  to  evaluate  the  results  at  the  interface  of  the 
annuli . 
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The  same  argument  can  be  extended  for  both  bounding 
expressions  through  the  entire  chamber  to  the  pickoff.  As  a 
result,  it  can  be  seen  that  these  two  hypothetical  solutions 
provide  bounding  expressions  for  the  sensor  gain  if  it  can  be 
established  that  the  expression  for  gain  behaves  in  a  mono¬ 
tonic  fashion  as  ki  varies  in  any  arbitrary  but  monotonic 
fashion  from  parabolic  to  uniform  where  ki  is  mad'-'  to  apply 
over  the  entire  range  of  r. 

This  is  an  unnecessarily  strong  condition.  It  does 
not,  for  instance,  take  into  account  the  fact  that  ki  varies 
in  a  monotonic  fashion  with  rx.  As  such,  this  method  could 
be  used  to  specify  output  bounds  if  ki  were  allowed  to  vary 
in  an  arbitrary  manner  with  n .  The  governing  equation  speci¬ 
fies  a  specific  manner,  not  an  arbitrary  one,  in  which  the 
profile -changes  from  parabolic  to  uniform  as  n  decreases.  But 
although  this  condition  is  not  a  necessary  one  for  the  hypo¬ 
thetical  solutions  to  be  bounds,  it  is  certainly  a  sufficient 
one.  It  is  the  one  pursued  here  because  it  was  felt  that  no 
other  method  could  be  used  which  would  not  require  the  solving 
of  (30) .  Such  an  effort  would  partially  defeat  the  purpose 
of  this  abbreviated  approach  to  finding  expressions  for  the 
sensor  output. 

It  will  not  be  possible  in  all  cases  to  establish  the 
monotonic  behavior  for  an  arbitrary  variation  of  kx ,  particu¬ 
larly  because  the  form  of  Yn  depends  on  ki .  In  those  cases, 
some  insight  can  be  gained  by  allowing  ki  to  vary  in  some  con¬ 
venient  manner  from  parabolic  to  uniform.  This  approach  is 
not  rigorous,  but  should  be  of  some  value  for  the  purpose  of 
obtaining  engineering  design  equations. 

Consider  first  the  expression  for  Gc  obtained  from 
(63),  (29),  (38),  (43),  (46^  ,  (52),  (53)  and  (59)  with  o>=0. 

Only  the  first  term  in  the  series  generated  by  U2  will  be 
considered,  since  it  is  shown  in  Section  2.3.2  of  the  text 
that  the  higher  order  terms  are  very  small.  The  coefficient 
C0  is  found  here  by  a  technique  analogous  to  that  used  when 
kx  =  kro.  That  technique  will  be  explained  when  equation 
(D-21)  is  found  in  the  following  appendix. 

J  i  h  rj 
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where 


%.o 


After  some  examination,  it  is  evident  that  the  com¬ 
bined  contribution  of  I2  and  the  exponential  term  diminishes 
as  ki  changes  from  parabolic  to  uniform,  as  in  the  case  of 
the  first  term  in  (C-l) .  However,  since  the  signs  of  the 
first  and  second  terms  are  opposite,  it  is  necessary  to 
establish  their  relative  magnitude  in  order  to  determine  the 
sign  of  the  net  change  in  G0  as  kx  is  varied. 

The  non-rigorous  method  of  choosing  kx  equal  to  some 
convenient  function  of  z  will  be  used  in  this  e'ffort.  This 
lack  of  rigor  may  not  be  so  disturbing  if  it  is  recalled 
that  we  have  established  from  intuitive  means  that  the  most 
attractive  operational  range  is  that  in  which  the  entrance 
effects  can  be  ignored,  or  that  range  in  which  the  second 
term  in  (C-l)  can  be  ignored. 

The  solution  of  (C-2)  will  now  be  given  for  the  range 
0  <  z  <  h/2  when  kx  is  chosen  to  be  some  convenient  func¬ 
tion  of  z.  The  solution  for  negative  z  can  be  found  from 
this  solution  by  symmetry.  The  boundary  conditions  are 


=0 

dz  z  =o 


(C-7) 
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Now  leu 


(C-9) 


where  1/2 >  m>  0.  The  upper  limit  of  m  has  been  chosen 
so  that  the  area  under  kx  for  m  =  1/2  is  equal  to  the  area 
under  the  parabolic  profile.  The  mean  deviation  of  this 
profile  and  the  parabolic  is  about  5  percent.  Comparable 
deviations  for  other  kx’s  might  be  expected,  yielding  an 
approximation  accurate  to  the  first  order. 


.First  let 


(C-10) 
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Then  (C— 2 )  becomes 

Now  let 

and 
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/— 7—  m  +  2 
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Equation  (C-ll)  becomes 


w2ilU  +w-^-  + 

w  dW2  +  w  dW 


••-ter)’ 


u  =  0. 


This  equation  has  the  solution 


2  0 


U=A'J  (W)+B'Y  (W), 

(■JS+F )  1  m+2  ) 


(C-ll) 
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(C-14) 


(C-15) 


where  A*  and  B'  are  constants  and  J  ?nd  Y  are  Bessel  functions. 
If  (C-12)  and  (C-13)  are  applied,  (C-15)  becomes 


+0 = A/y* 


m+  2 
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i  /2V 
+B'y^Y,_L_.(-;rr 
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(C-16) 


and 


The  boundary  conditions  are  given  by 
d^0 


dy  | y  —  1 


=  0 


*0  (0)  =  0. 


(C-17) 


(C-18) 


Equation  (C-18)  requires  B’  =  0.  Equation  (C-17)  can 
now  be  used  to  determine  Bq. 

If  the  series  representation20  of  the  Bessel  function 
is  used,  (C-16)  becomes 
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where 


f/(y)  = 


m  4-  2 
m  +  I 


,  m  +  2 


m+2 


■] 


(C-24) 


The  function  f ' (y)  has  a  value  of  one  at  the  center  plane  for 
all  m  and  zero  at  the  vail,  but  its  average  decreases  over  the 
range  of  y  is  less  than  two  percent  as  m  increases  from  zero 
to  one-half. 


Finally 


h  _ 
2(m  + 1) 


(C-25) 


If  (C-6)  and  (C-22)  through  (C-25)  are  used,  the 
second  term  in  (C-l)  becomes 
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*) 


where  X0  is  found  from  (C-21) . 

Both  and  f 1  are  normalized  functions  whose  maximum 
value  is  independent  of  m  and  whose  shape  depends  only  weakly 
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on  m.  Thus,  the  ratio  of  the  integrals  may  be  considered  a 
constant.  Since  the  above  approach  is  identical  when  m  =  0 
to  the  method  used  to  obtain  the  uniform  solution,  the  value 
of  this  ratio  can  be  found  from  (78)  and  (C-26)  with  m  =  0. 
It  can  be  shown  that  for  m  =  0, 

x  .  it  (027) 

Ao  "  2 

and 


(The  behavior  of  the  second  term  in  (C-l)  can  then  be 
found  from  (C-26)  and  (C-28) .  It  was  calculated  for  Re^  =7, 
the  worst  case.  It  was  found  to  vary  from  .081  h  for  m  =  0  to 
.  0°9  h  for  m  =  1/2.  The  behavior  was  very  nearly  linear  with 
m.  Although  this  behavior  is  somewhat  weak,  it  should  be 
compared  to  the  behavior  of  the  first  term  in  (C-l)  to  estab¬ 
lish  whether  the  entire  expression  increases  as  m  varies 
from  0  to  1/2. 

Now  consider  the. first  term. in  (C-l).  If  this  is  re¬ 
written  in  terms  of  ki,  then 


(C-29) 
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It  can  be  shown  that  (C-30)  behaves  in  a  nearly  lin¬ 
ear  fashion  with  m,  and  that  this  expression  varies  from 
.167  h  to  .225  h  as  m  varies  from  0  to  1/2.  It  can  be  shown 
from  (78)  that  (C-30)  should  be  .167  h  when  the  profile  is 
uniform  and  from  (64)  that  it  should  be  .243  h  when  the  pro¬ 
file  is  parabolic.  This  approach  would  then  seem  sufficient¬ 
ly  accurate.  It  shows  that  the  increase  with  m  in  the  first 
term  in  (C-l)  is  about  four  times  greater  than  the  increase 
in  magnitude  of  the  second  term.  Thus,  Go  very  probably 
increases  monotonically  as  ki  changes  from  uniform  to  para¬ 
bolic  even  in  the  upper  limit  of  the  fully  developed  range. 

Let  us  next  consider  when  Re'  is  very  small.  If 
(52)  is  used,  one  then  has 

cosh(o  Z  )  ,4-  I  [V 

I  2  k,  cosh/ ah  \  I  cosh2  an  f  kicosh  (o£)d£d^d« 

„  I  Jo  1  ixil  z _ /J0 _ (C-31) 
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For  small  w ' ,  this  may  be  expressed  as 
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The  phase  lag,  <j>,  of  Giw  at  small  w’  is  thus  given 
in  radians  approximately  by 


♦  -T1* 


(C-34) 


and  the  amplitude,  jGiJ,  is  given  approximately  by 


|G 


(C-35) 


Consider  first  (C-34) .  It  can  be  established  in¬ 
tuitively  that  the  value  of  I3  decreases  as  ki  changes 
from  parabolic  to  uniform.  To  see  this,  it  may  help  to  think 
of  I 3  in  the  form 


n_  n_  y 

Jz  k l{z)p[z)  JZ  f  k,(e)df d,dz 

=  ~7h  75  75  (C-36) 
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4z2 

where  p(z)  3  1 - r-5-  is  looked  on  as  a  weighting  function.  If 

p(z)  were  a  constant,  .then  the  value  of  I3  would  be  independent 
of  ki (z) .  But  since  p(z)  decreases  monctonically  from  the 
centerplane  to  the  wall  in  this  case,  the  value  of  I3 
increases  if  ki  is  weighted  more  strongly  toward  the  center- 
plane.  Thus,  I 3  increases  as  the  profile  changes  from 
uniform  to  parabolic.  This  implies  that  the  parabolic  solu¬ 
tion  provides  an  upper  bound  for  the  phase  lag  when  o>'  and 
Re 5  are  small,  and  the  uniform  provides  a  lower  bound. 

It  would  appear  that  (C-35)  cannot  be  handled  in  any- 
convenient  manner  for  a  general  k3 .  Let  us  then  attack  the 
problem  in  a  less  rigorous  manner  by  assuming  that 
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where  2*  m<  » .  The  upper  limit  for  m  represents  a  uniform 
profile. 

1 1,  is  then  given  by_ 

/'  (l-2'ra)(l-Z'2)(l-  Jz'2)  fV(\~  ('"idt'dlj'dz'  (C 

^4~  si  r\  rv'  _ 
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where  all  tho  independent  variables  have  been  normalized  with 
respect  to  h/2.  It  can  then  be  shown  that 
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It  can  also  be  shown  that  if  C-37  is  usea  r„  becomes 
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and  that  I 3  becomes 
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If  (C-39) ,  (C-40)  and  (C-41)  are  used,  (C— 3 5 )  can  be 
evaluated.  It  can  be  shown  that  the  bracket  to  the  right  in 
(C-35)isa  weakly  increasing  function  of  m.  Thus,  the  attenu¬ 
ation  predicted  by  either  expression  can  be  considered  as  a 
good  estimate  of  the  behavior  at  small  u>*  and  Reo,  although 
it  would  appear  that  the  parabolic  solution  provides  a  max¬ 
imum  estimate  for  the  amplitude  of  Gi^  and  the  uniform 
provider,  a  minimum 

When  Re*  is  very  small  and  to1  is  large,  (C-32)  can 
then  be  approximated  by 
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The  influence  of  the  functions  of  w  in  this  expression  is  now 
more  apparent.  Clearly,  the  exponential  functions  with  a 
real  exponent,  which  may  be  regarded  as  weighting  functions, 
cause  the  integrands  to  be  weighted  quite  heavily  toward 
the  values  near  the  wall.  Then  if  the  value  of  k*  is  made  to 
increase  near  the  wall,  the  magnitude  of  the  numerator  in¬ 
creases  more  than  that  of  the  denominator,  which  has  no  such 
weighting  function.  Thus,  the  amplitude  of  G ^  increases  as 
the  profile  changes  from  parabolic  to  uniform  for  large  to' 
and  very  small  Re£. 

The  phase  can  be  discussed  by  considering  only  the 
numerator,  since  the  denominator  is  real.  For  this  dis¬ 
cussion,  it  may  help  to  consider  the  function  ki  as  the 
weighting  function.  If  the  value  of  ki  near  the  wall  in¬ 
creases,  then  the  contributions  due  to  the  exponential  that 
has  an  imaginary  exponent  are  weighted  more  heavily  near  the 
wall.  Since  the  phase  lag  relative  to  the  wall  decreases  as 
z  approaches  the  wall,  this  implies  that  the  phase  lag  de¬ 
creases  as  the  profile  changes  from  parabolic  to  uniform. 

It  has  been  shown  that  when  Reo  is  very  small,  the 
amplitude  of  G^  for  small  u'  does  not  vary  significantly  as 
the  radial  profile  is  changed  and  that  for  larqe  u'  the 
amplitude  is  bounded  above  by  the  uniform  solution  and  below 
by  the  parabolic.  Since  the  amplitude  of  the  functions  of  a 
(that  is,  w ')  in(C-31)  is  monotonic  with  to',  it  would  then 
seem  reasonable  to  conclude  that  for  all  id',  either  the  ampli¬ 
tude  is  bounded  above  by  the  uniform  solution  and  below  by 


the  parabolic  or  the  difference  in  the  two  expressions  for 
amplitude  is  negligibly  small. 


If  Re£  is  not  small,  then  the  contributions  due  to 
U2  must  be  considered.  Before  this  problem  is  attacked  head- 
on,  a  little  discretion  can  be  used.  The  form  of  can  be 
represented  as 


where  Bq  is  the  first  eigenvalue  of  the  equation 

d2*« 

*(- 
V  n 

The  function  F0(Reo,oO  represents  the  output  due 
to  the  \p0  contribution  to  the  boundary  condition  at  the 
coupler,  and  the  exponential  term  specifies  the  attenuation 
and  phase  shift  that  results  when  this  signal  is  trans¬ 
ported  to  the  pickoff .  The  function  F0  (Rep, in')  will  be  dis¬ 
cussed  later . 

First  consider  equation  (C-44) .  The  only  signifi¬ 
cant  variation  in  ki (2)  occurs  near  the  wall,  where  is 
small.  The  effect  in  this  equation  of  changing  the  radial 
profile  from  parabolic  to  uniform  is  then  small.  Thus,  the 
shape  of"  iJ-q  and  the  value  of  Bo  are  relatively  independent 
of  the  shape  of  ki .  The  solution  for  a  uniform  radial  pro¬ 
file  (to  be  shown  later  in  appendix  F)  can  be  used  then  to 
show  that  the  imaginary  part  of  Bq  is  positive  and  nearly  a 
linear  function  of  u'  and  that  the  real  part  is  positive 
and  nearly  independent  of  in'. 

Some  features  of  F0(Reo,cO  can  now  be  discussed. 
First,  since  the  shape  of  v0  does  not  change  greatly  with 
kj  (z) ,  F0(Reo,w')  does  not  depend  significantly  on  ki (z) 
either.  Its  behavior  can  then  be  treated  simultaneously 
for  all  ki  (z) .  To  treat  it,  we  musr.  understand  what  the 
contributions  due  to  u2  represent. 

Notice  that  the  boundary'  conditions  on  i{in  are 
homogeneous.  It  is  as  if  some  signal  were  introduced  at 
Ro  in  a  stationary  sensor  and  made  to  propagate  toward  the 
pickoff.  The  time  dependent  motion  of  the  sensor  would 
effect  the  output  due  to  this  signal  only  in  that  the  input 
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signal  would  vary  due  to  a  change  in  the  function  which  U2  must 
satisfy  at  the  coupler.  This  input  signal  is  just  the  value  of 
ui  relative  to  the  wall. 

Now  consider  u2,  as  given  by  (46),  (52)  and  (53).  If 
a)'  is  very  large,  (53)  becomes 

'  g,(2)=^  — 

1  Jz  Jo 

Since  the  dummy  variable  n  is  the  upper  limit  for  the  dummy 
variable  £,  the  value  of  is  always  negative,  and 

_  V  2  '  h  ~  h ' 

in  particular,  increases  in  magnitude  as  w'  increases.  The 
exponential  function  with  the  real  exponent  then  causes  the 
integrand  to  diminish  rapidly  as  o' 1  increases.  Thus,  gi  van¬ 
ishes  for  large  ui’,  and  Ui  takes  on  the  behavior 


f 

5 


(C-46) 


For  z  ^  h/2,  g2  vanishes  for  large  os'.  For  points  away  from 
the  wall,  then,  ui  goes  to  zero  as  o>'  increases.  In  that 
case,  u2  must  then  take  on  the  value  Rq  at  the  coupler  for 
points  away  from  the  wall. 

It  will  be  shown  in  appendix  D  that  if  the  radial 
profile  is  assumed  parabolic  throughout  the  chamber  when 
oi '=0,  Uj  is  given  by 


u,= 


5  Re  * 


16  V  5h2  /  r  +Rn 


(CHI  7) 


In  appendix  F  it  will  be  shown  that  if  the  radial  profile  is 
assumed  uniform,  when  oi  ’=0 


u.= 


4  V  h2  i  r  t  R0 


(CH8) 


Thus,  the  magnitude  of  the  boundary  condition  at  the  center- 
plane  which  u2  must  satisfy,  changes  from  about  (-Rep)Rp  for 
small  oi '  to  Rq  for  very  large  oi  * .  4 

These  facts  can  be  used  to  approximate  the  magnitude 
of  Fc'(Rep,  oi* ),  as  indicated  in  Figure  C3.  The  magnitude  at 
small  oi'  is  a  linear  function  of  Rep,  but  approaches  a  value  of 
1  at  large  oi’  for  all  Reg. 

Now  consider  the  exponential  term  It  was  noted 
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previously  that  the  real  part  of  Bois  nearly  independent  of 
to*.  Then  the  magnitude  of  the  product  of  the  exponential  and 
Fc(Re'o»w‘)  approaches  a  constant  value  as  to*  increases.  Thus, 
the  signal  contribution  due  to  U2  vanishes  for  large  co*  only 
if  the  exponential  term  vanishes.  This  implies  that  Reo  must 
be  sufficiently  small  to  ensure  sufficient  transport  attenua¬ 
tion  to  reduce  the  signal  contribution  below  the  noise  level 
for  all  to 1 . 

Since  both  Bo,  and  i Pq  are  nearly  independent  of  kx  (z) , 
the  behavior  of  the  signal  contribution  due  to  u2  can  be  dis¬ 
cussed  simply  by  noting  the  effect  on  the  exponential  term  of 
changing  the  integral  in  the  exponent.  The  integral  increases 
monotonically  as  kx (z)  changes  from  parabolic  to  uniform. 

If  the  behavior  of  Bo  with  u> '  is  recalled,  it  can 
then  be  shown  that  the  parabolic  solution  provides  a  mini¬ 
mum  expression  for  phase  lag  and  a  maximum  for  amplitude  in 
.  the  signal  contribution  due  to  u2  .  The  uniform  solution 
has  the  complementary  features. 

It  should  be  noted  that  the  sign  of  the  contribution 
due  to  u2  is  negative,  while  that  due  to  ui  is  positive.  This 
fact  must  be  recalled  carefully  when  the  behavior  of  G^u  is 
discussed  in  a  range  where  both  contributions  must  be 
considered. 

Now  the  behavior  of  G can  be  discussed  if  several 
other  points  are  mentioned.  '  As  n)r  increases,  the  amplitude- 
due  to  Uj  diminishes  relative  to  dc  .  and  approaches  zero, 
while  the  phase  lag  increases  and  approaches  some  finite 
limit.  The  phase  of  the  signal  due  to  u2  increases  without 
limit  as  a) '  increases,  and  the  amplitude  approaches  a  finite 
limit.  This  vaiue  may  be  either  greater  or  less  than  the 
value  at  dc  ,  depending  on  the  value  of  Re0 . 

These  results  can  be  summarized  by  the  statement  in 
Section  2.3  of  the  main  text. 


APPENDIX  D.  SOLUTION  FOR  A  PARABOLIC  RADIAL  PROFILE 


integral 


if  (33)  is  used,  (53)  can  be  found  by  evaluating  the 
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The  following  evaluations*  are  needed 

♦The  number  to  the  right  of  each  evaluation  identifies  the  equation  in 
reference  21. 


J cosh2 v  dw  =y+  [sinh(2w)j  , 
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where  B2k  is  Bernoulli's  coefficient, 
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/tanhwdw=^n(cosh  w)  , 
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,  a>  22k  (22k  -I)  B-k  w2k+2 

Jw2  tanh  w  dw=  £ - p— — =: - 

k=l  (2k+2)[(2k)!j 


as  well  as  the  relationships 


sinh  (2w) = 2  sinh(w)  cosh(w) , 
cosh  (2w)  =  2  cosh2  (w)  - 1 . 


2.479  #7 


It  can  be  shown  that 


From  (D-l)  and  (D-2) ,  as  well  as  the  relationships 
arid  evaluations  above,  it  can  be  seen  that 


If  this  result  and  (34)  are  used,  (53)  becomes 
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The  value  of  0gi,  that  is,  the  value  of  91  as  u  goes  co  zero, 
can  be  found  by  reevaluating  Ii  with 


The  result  is 


cosh 2 (<>17)=' cosh2  (n(  )  =  I  , 


(D-^5) 


Equation  (D-4)  can  be  checked  by.  showing  that  (B-5)  'can  ba- 
obtained  also  by  finding  the  limit  of  gx  as  a  approaches  zero, 
in  (D-4).  In  this  limiting  process,  only  the  first  two  terms 
of  the  series  representation  of  the  hyperbolic  tangent 
need  be  considered  as  significant: 

♦anh  x  —  x  — ^  , 

In  order  to  solve  (56)  when  kx  is  given  by  (33) , 
consider  the  substitution 


(D-6) 


where  Yn  is  some  constant  as  yet  unspecified.  If  (D-6)  is 
differentiated  twice  with  respect  to  z,  the  result  is 

ft2  llr  r  #  *  .t  v  —2 


Equation  (56)  will  be  satisfied  if 
<l2(9s,n  ...  , 
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(D-8) 


0n  +o2  =  2yn  +  £n 


(D~9) 


and 
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wherd  £n  is  another  constant  as  yet  unspecified.  It  will 
be  seen  that  it  is  the  eigenvalue  of  (D-8)  defined  by  the 
boundary  condition . 


(93)n(±!)=0- 
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Equations  (D-9)  and  (D-10)  can  be  used  to  find  t*n 
andYnin  terms  of  a  and£n-  The  result  is 
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The  positive  sign  has  been  chosen  in  front  of  the 
radical  in  (D-12.)  and  (D-13)  for  convenience.  The  solution 
obtained  in  this  manner  must  be  unique,  since  the  original 
problem  is  well  posed.  The  solution  obtained  for  (g3)n  does 
depend  on  the  sign  chosen;  however,  when  the  accompanying 
change  of  Yn  is  considered  in  (D-6) ',  the  solution  obtained 
for  ijjn  by  either  method  should  be  formally  identical.  No 
effort  is  made  to  pursue  this  point,  though. 

Now  consider  (D-8) .  If  the  substitution 


*=-2V‘ 


is  used,  this  equation  becomes 


(D-14) 


(D-I5) 


This  equation  is  the  confluent  hypergeometric  equation.  The 
general  solution  is  given22  by 


g3=  K'M 


Cn  I 

■  t  —  t  y 

8/n  2 


K**2  m  (-4fn+i),l,x 


(IV 16) 
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where 


(D-17) 


u  r  r(V-w)r(ag)*m 

[at.  a2,x]  r  (a2*  m)r(o,)  mi 


Since  the  boundary  conditions  are  symmetric  with  respect  to 
the  centerplane ,  K"  =  0. 

Equation  (D-ll)  can  now  be  used  to  find  the  eigen¬ 
values  of  (D-16) .  It  is  then  found  that 


*n=  e 


4  (*)'  i 


M  n  1  IJL  Z  \2 

ML'pn’2’-T-liry  J 


(D-18) 


where 


yn  =  -{l+/«  +  '«,+  Sn') 


p  _ _ 

8  |l  +  A/l+iu'+Sn' ) 


and  the  are  the  ordered  set  of  values  of  S  which  satisfy 
the  equation 


M  __ — i+VTw^s1)  -°> 

8(l  +  vm«+S  ) 


(D-X9) 


If  (40),  (43),  (46),  (52),  (59),  and  (D-12)  are 
used,  it  can  be  seen  that 


r  n  >‘-[?J>L  +  r  ^  cosh (02) 
R°n°  ‘  V  R0  °/  Cosh(-^) 


I  2A©h2 

+  -k  l  c„e  *„ 


(D-20) 


where 


i«'  +  Sn  +  2(l  +  ./i+iw'+Sn  ) 


xii®, 


If  both  sides  of  (D-20)  are  multiplied  by  (l-4z2/h2)i|»n 
and  integrated  with  respect  to  z  from  0  to  h/2,  che  orthogo¬ 
nality  condition  provides  an  evaluation  for  Cn,  namely 


Since (R2/3ReoRo ) is  of  the  order  of  10  ,  the  expo¬ 

nential  in  the  series'above  can  be  set  equal  to  1  if  An  is  of 
the  order  of  102  or  less.  Since  it  will  be  seen  later  that 
only  the  first  term  in  this  series  need  be  considered  and 
that  A  is  about  11  when  w'=0 ,  the  exponential  in  (D-23)  is 
set  equal  to  1.  The  term  linear  in  Ri  can  be  neglected,  and 
thus  (D-23)  becomes 


(D-24) 


one  obtains 


i 


! 

I 


I 


-X, 


c»r('-^)^=-45'3Re° 


cosh 


(f Jus) 


/  (l-w2)  <MW 


3Re'0  /  . 

~ 2 —  $2  (*,«')  + 1 


-I 


(l-w2)fndw 


(D-30) 


*ndw 


where 


.  X‘(— •) 


-2iT('-*::)+(3+t) 


I  j  j  \  tanh  (?y^) 

vA? 


27 

and  ij/n  is  given  by  (D-18)  with  the  substitution  W  =  -r—. 


If  (D-18),  (D-22) ,  (D-25) ,  (D-23) ,  (D-28)  and  (D-30) 
are  used,  (61)  can  be  evaluated.  If  the  limit  for  this  ex¬ 
pression  as  oj  approaches  zero  is  found  (with  the  help  of 
(D-29))  (63)  can  be  evaluated.  This  evaluation  is  given  by 
(64).  If  (63)  used  to  specify  Xi  such  that  $  approaches  1 
as  a)  goes  to  zero,  then  (65)  is  obtained. 
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APPENDIX  E.  COMPUTER  PROGRAMS  FOR  EQUATIONS  (64)  AND  (65) 

The  following  are  the  computer  programs  used  to 
evaluate  *.(64)  and  (65)  numerically. 

The  program  to  evaluate  Gq  is  listed  first.  One 
data  card  is  required  to  specify  ReZS  (starting  value  of 
Re^) ,  REZF  (final  value  of  Re©) ,  and  DR  (the  increment  be¬ 
tween  the  values  desired) .  Total  running  time  of  this 
program  for  thirteen  values  of  Re^,  for  instance,  was  1.4 
minutes  on  an  IBM  7094.  The  results  will  be  plotted  by  a 
CalComp  plotter  if  the  cards  calling  for  the  subroutines 
PLOT,  PLOTS,  SCALE,  AXIS,  and  FLINE  axe  retained. 


The  program  to  evaluate  Gj.^  is  listed  next.  Double 
precision  is  required  in  evaluating  $2  (w,w  *)  because  the 
series  representation  for  the  hyperbolic  tangent  must  be 
carried  out  to  an  equivalent  of  at  least  the  ninth  degree 
of  the  p.olynomial  expression  to  obtain  any  contributions  to  . 
Giu.  The  subroutine  HOME  searches  for  the  first  eigenvalue 
of  (D-19) . 


Two  data  cards  are  required,  one  for  the  title 
(called  "Frequency  Response"  here)  and  one  to  specify 
REZERO  (value  of  Re'),  WPS  (starting  value  of  w'),  and 
DWP  (the  increment  between  the  values  desired) .  The  re¬ 
sults  will  be  plotted  by  a  CalComp  if  the  call  cards  for 
the  same  subroutines  previously  mentioned  and  also  the 
subroutine  SYMBOL  are  retained. 

Typical  running  time  of  this  program  is  about 
3  minutes  on  an  IBM  7094. 
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PROGRAM  TO  PLOT  GO  VS  REZERO 

DIMENSION  I8UFI 1024)«R( 100), GO! 100) 

CALL  PLOTS! IBUF* 1024) 

CALL  PL0TI0.,-30.,-3l 
CALL  PLOTtl.r .75,-3) 

REAL  L AMO AO 
LAMDA0=11 .311052 
V 6*. 08936417 
V7=. 5054014 
'  V8*. 4208766 
C6*. 026984127 

9  READ  (5,10)  REZS,REZF,DR 
10  FORMAT  (3F10.0) 

REZS  IS  STARTING  VALUE  OF  REZERO 
REZF  IS  FINAL  VALUE  OF  REZERO 
DR  IS  THE  INCREMENT  FOR  REZERO 

Insert  a  blank  card  after  last  data  card 

IFIDR.EQ.O.)  GO  TO  50 
N*  !  REZF-REZS 1 /DR* 1. 5 
DO  20  1*1 tN 

R  ! 1 l*REZS*FLO AT  f 1-1 1 *0R 
C1*3.*R! I) 

C3*EXP (— LAMDAO/C 1 ) 

C4*C3*V6*!V7/V8) 

C5*C4*.25 

C7*C6-C5 

20  G0(I)=R1I)*9.*C7 

"WRITE  16,21)  !R!D,60U),I*1,N) 

21  FORMAT! 1H1,5X«7HRE-ZER0»9X,6HG-ZERQ//!  IX,  IP2E15.5) ) 

"  PLOT  SECTION  " 

CALL  SCALE !R, 8. ,N, 1 I 
CALL  SCALE  !G0,  9.,N,1) 

CALL  AXIS!0.,0.,7HRE-ZER0,-7,8.,0.,R!N*l),R!N+2)) 
CALL  AXIS!0.,0.,6HG-ZER0,  6,  9.,90.,G0!N+l),G0!N+2) ) 
CALL  FLINE!R,GO,-N, 1,0,0) 

CALL  PLOT ( 10. ,0. »-3) 

GO  TO  9 

50  CALL  PLOT(0.,0.,999) 

STOP 

ENO 
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£ 

I 
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C  PROGRAM  TO  PLOT  OB  ANO  THETA  VS  OMECA 
C 

DIMENSION  I8UF( 1024) , HP! 100) ,R! 100) .DEGREE! 100) 

CALL  PLOTS! 1BUF, 1024) 

CALL  XMAX! 125.) 

CALL  PLOTIO.,-30.,-3) 

CALL  PL0TI1... 75,-3) 

EXTERNAL  EVAL 

COMPLEX  Z 1, Z2, Z3.Z4.Z5.Z6.Z7, LAMBDA, P, GAMMA, Z8, 29, Z 10, Z11.Z12.M 

COMPLEX  Z 13, Z 14, PS! ,Z15,Z37  , Z16.Z17, Z18,Zl9,CC0SH,Z20,Z2l, Z2? 

COMPLEX  Z23,PHI2,Z24,Z25,Z26,Z27,Z28,Z29,Z30,Z31,Z32,Z33,Z34,Z35 

COMPLEX  PH11, THETA 

COMPLEX  Z36 

COMPLEX  S 

COMMON  SPART 

COMMON  OMEGA 

COMMON  A1,A2 

REAL  LAMDAO 

DIMENSION  SPART !2),DELTA!2) 

DIMENSION  LABEL  16) 

DIMENSION  Y1!101),V2I101),Y3!101),Y4!101),Y51101),Z! 101 ) 

DIMENSION  Y9! 101 ) 

51  FORMAT  !6A6) 

52  FORMAT  I IH1.48X.6A6///////) 

53  FORMAT  11H  ) 

54  FORMAT  !1X,57X,8HREZER0  =,F10.6) 

55  FORMAT  11X,36X.5H0MEGA,BX,1HX,11X,1HY,10X,3HRH0,8X,2HDB7X,5HPHASE/ 
2/) 

56  FORMAT  I IX,35X,F6.3,3F12.5,F10.3,F10.1) 

10  FORMAT  I8F10.0) 

H=0.0l 

LAMDAO* 11. 31 1052 
V 6*. 08936417 
V7*« 5054014 
V8*. 4208766 
C6=. 026984127 
READ  15,51)  LABEL 
9  READ  15,10)  REZERO.MPS.HPF.DHP 
C 

C  UPS  IS  STARTING  VALUE  OF  OMEGA 
C  UPF  IS  THE  FINAL  VALUE  OF  OMEGA 

C  DUP  IS  THE  INCREMENT  FOR  OMEGA 

C 

C  INSERT  A  BLANK  CARO  AFTER  LAST  OATA  CARO 

C 

IF  iOMP  .EQ.O.O)  GO  TO  100 

11  C1*3.0*REZER0 
C2*2.0/Cl 

C3*EXPI -L AMOAO/C 1 ) 

C4=C3*V6*IV7/V8) 

C5*C4/4.0 
C7*C6— C5 

MRITE  16,52)  LABEL 
URITE  ! 6, 54 )  REZEKQ 
HRITE  16,53) 

URITE  !6,55) 

L* ! UPF-UPS ) /OUP» l . 5 
00  8  J=I,L 

UP 1 J ) =UPS»FLOAT i J-l ) *DUP 
OMEGA*MP! J) 

SPARTI1 >*4.584669 
SPART! 2) *0.0 
DELTA! l)*0.l 
DELTA12 1*0. 1 
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^  T-yfr-jt  j-j  - 


aT.»  JifAAG. 


0*1.0E*36 

N*2 

I *3000 

CAU.  H0HE(EVAt*SPARTtDELfA*M»I*0) 
S*CHP1X( A1.A2 ) 

Z1*CMPLX( 0.0*  OMEGA) 

Z2=Z1*S 

23*1.0*22 

Z4=CSQftT(Z3> 

25*1.0*24 

Z6*2.Q*Z5 

Z 7*8. 0*25 

LAMBDA=22*Z6 

P=S/Z7 

GAMMA*-25 

Z8=-P 

Z9*-GAMNA/2.0 
210*-29/2.0 
1*1 
W*0.0 
5  W2=M**2 
M3=1.0~W2 
Zll=Z9*W2 
Z12=M(28,0.S,21l) 

2l3*2lO*H2 

Z14=CEXP(213) 

PSI*212*214 
Z15*W3*PS1 
Yl  ( 1 1**5EAL(2151 
Y2t I )=AIMAG( 215) 

216=215*PSI 
Y3f I )=REAL(Z16) 

Y9( I )*AIMAG( Z16) 

Z17=CS0»TI21) 

218*217/2.0 

219=CC0SH(218) 

220*218*11 
Z2l=CCUSH(Z20 ) 

222=221/219 
223=222 *PH 12 fH, OMEGA) 

224=222-1.0 
225=C2*224 
226=223*225 
Z27=226*H3*PSl 
Y4C I  )*REAC(227) 

Y5 ( I )=AIMAG( 227) 

IF  (I.EQ.101)  GO  TO  6 
1*1*1 
M=M*H 
GO  70  5 

6  CAU.  QSr ?H,Y1,2, 101) 

VWU>  ! 

CAU  U'<  «H,Y2,2,101) 

92*2  HOD 

C&U  0SF(H*Y3»2.  101) 

Yi*2(101) 

CAU  0SP(H.Y4,2,101) 

.  .=2(101) 

C-U  OSFIH. Y5 1 2.  101 ) 

95=2(101) 

CAU  QSFIH.Y9, 2,101) 

V9*2 (101) 

228*C«PLX(W1,V2) 

237=CMPLX(V3,V9) 

Z29=LAHB0A/C1 
230*CEXP(— Z29) 

231*Z2«/Z37 
Z32=CMPLX( V4, V51 


’f 

g. 


s 


-c 

i 

ii 
8C 


I 

§ 


§1 


c 

c 

c 


100 


Z33=Z30AZ32*Z3l 
Z3A=Z33/4.0 
Z35-PHI 1 ( ONEGA) 

Z36=Z35-Z3* 

THETA*Z36/C7 
X=»REAL(  THETA) 

Y«AIMAG(THETA) 

RHO*CABS( THETA) 

RU)*20.0*AL0G10(RH0) 

ANGLE* ATAN2(Y,X) 

DEGREE ( J ) *57 • 2957795*ANGLE 

IFt ANGLE. 6E.0.)  DEGREE) J )=OEGR£E( J )-360. 

WRITE  (6.56)  OMEGA. X,Y,RHO,R(J),DEGREE( J) 

PLOT-  SECTION 

CALL  SCALE  (WP.8..L.1) 

CALL  SCALE  <R.  9..L.1) 

CALL  SCALE  (OEGREE,  9..L.1) 

CALI  AXIS  (0. .0. »5H0MEGA,-5»  8.  ,0.»  WP  (L*-l ) »  WPtL+2 ) ) 

CALL  AXIS  I0.,0.,2HDB,2,  9. ,90. ,R(L+1) ,R(L*2 >) 

CALL  AXIS  (B.,0.,5HPHASE,-5,9.  ,90. .DEGREE IL*1) ,DEGREE(L+2 > ) 
CALL  FLINE(WP,R,-L, 1.1,0) 

CALL  FLINE(WP.OEGREE,-L* 1.1,11) 

CALL  SYMBOL  ( .5,9.8,.1,0,0.,-1) 

CALL  SYN80L  ( 999. ,999., . 1,3H  OB.O.,3) 

CALL  SYMBOL  ( .5,9.5, .1, 11,0. ,-l ) 

CALL  SYMBOL  ( 999. ,999., . 1.6H  PHASE, 0., 6) 

CALL  SYMBOL (2., 9. 8, . 1, LABEL, 0«, 36) 

CALL  SYMBOL  (2.,9.5,.1,8HREZER0=  ,0.,8) 

CALL  NUMBER  ( 999. ,999., . 1, REZERO, 0., 6) 

CALL  PLOTUO.,0.,-3) 

GO  TO  9 

CALL  PlOT (0.,0.« 999) 

STOP 

END 
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SUBROUTINE  EVAUQA.QB,*) 
COMMON  SPART 
COMMON  OMEGA 
COMMON  A1,A2 
DIMENSION  SPART (2) 

COMPLEX  G 
COMPLEX  S,T 

S=CMPLX  t  SPART ( 1 ) , SPART ( 2 ) ) 
T=G(S, OMEGA) 

OB=CABS(T) 

IF  (QB.GE.dA)  RETURN 
A1=SPART( 1) 

A2=SPART(2) 

RETURN  l 
ENO 
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SUBROUTINE  HOMEIQSET, EX, OS, N.1MAX.Q) 


DIMENSION  EX( 1) , 

IF  (N.LE.O)  RETURN 

1MAX?*IMAX 

FACT=.5/FL0AT(N) 

DO  10  1*1, N 


DS(1),P(9),0D(9), XTEM(9) ,PP(9) ,MUG(9) 
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HUGH  1*0 
P(i)*OS(I) 

10  0D(l)*4.*DStI) 

12  *  3 

13  *  N 
1*1 

00  150  17*1, IHAX7 

IMAX*17 

DEL0*0. 

00  90  16*1, N 

I  *  1-1 

IF(I.EQ.O)  1*N 

IF(MUG(I)  .EQ.  0  .OR.  N  .EQ.  1)  GO  TO  15 

HUGd>*0 

GO  TO  55 

15  XTEN(l)  *  EX(  I ) 

IF( I2.NE.3)  GO  TO  20 
OS(I)  *  .25*00(1) 

00(1)  *  0. 

20  Pd)*SIGN(DS(I), Pd)) 

?icrc*i 

60  XP*XTEN( I )*P( I) 

IF(XP.EQ.XTEM(I ) )  GO  TO  50 
EX( 1 )aXP 

CALL  OStT(0,QT,S60) 

GO  TO  ( 35,70,25) «NCYC 
25  EX( 1 )*XTEH( I ) 

SO  TO  90 
35  Q1  *  QT 

P(I>  *  -P(I) 

NCYC*2 
SO  TO  60 

70  P(I)=P(l)*.5*(Ql-0T)/;01*QT-2.*0) 

NCYC*3 
GO  TO  60 
50  NUG(I)*1 

IF(NCYC.EQ.3>  GO  TO  90 
55  13*13-1 

IF( I3.EQ.0)  RETURN 
SO  TO  90 

60  IF1QT.EQ.0. )  RETUP.'( 

13=N 

IF(Q-QT.LE.OELQ.OR.N.EQ.l)  GO  TO  80 
OELQ  =  Q-OT 

II  =  1 
80  Q  =  QT 

90  00(1)  =  DD(I)+ABS(FLOAT  (6-I2)*P(I)) 

I3*N 

12  =  12-1 
IF( I2.EQ.0)  12*3 
IF(N.EQ.l)  GO  TO  150 
PR00*1. 

00  100  15*1, N 
PP( I )*P( I )*PROO 
PR00*PR00-FACT 
IF(I.EQ.N)  1*0 
100  1*1*1 

00  160  15*1,10 
00  120  1*1, N 
110  XTEM( 1 )*EX( I ) 

120  EX( 1 )  *  EX( I )*PP( I ) 

CALL  QSET(Q,QT, $130) 

00  165  1*1, N 
165  EX( 1 )  *  XTEM(I) 

SO  TO  150 

130  IF( QT.EQ.O. )  RETURN 

IFIO-QT.GT.OELQ)  GO  TO  160 


116 


SffiSKT&& 


Q*0T 

GO  TO  150 
140  0«QT 
150  1*11 
RETURN 
ENO 
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COMPLEX  FUNCTION  CCOSH(Z) 
COMPLEX  Z 

CCOSH* (CEXPI Z )*CEXP(— Z) 1/2.0 

RETURN 

ENO 
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COMPLEX  FUNCTION  Ml AZERO.BZERO. ZZEROJ 

COMPLEX  AZEROtZZERO«A«Z  tCiLOGZRO.LOG At LOGZjLOGCtLOGX* Xt CLOG 
COMPLEX  CT 

REAL  LOGKt LOGS, LOGO *LOGN 

REAL  K,N 

REAL  1MAGA.IMAGC 

M*i.O 

X*l.O 

4*0 

A*AZERO 

B*BZERO 

Z-ZZERO 

IF  (CABS(ZZERO).EQ.O.O)  RETURN 
X*(A*£)/8 
M*M*X 
C*A*l.O 
0*8+1. 0 
N*2.0 
GO  TO  2 
4  M*M+X 

IF  (CABS(X).LT.1.0E-091  RETURN 

N*N»1.0 

C*C+1.0 

0*0*1. 0 

GO  TO  3 

2  L0G2RQ*CL0G( ZZERQ) 

IF  ICABSI A) )  12*13.12 
12  REAlA*REAL(A) 

'1MAGA*AIMAG(A1 
IF  (REALA.LT. 0.0)  GO  TO  16 

18  LOGA*Cl.p'  :a) 

LOGK*AL JG(K) 

L0G8*AL0G(8) 

L0GZ*CL0G(Z) 

3  LOGN*ALQG(N) 

IF  (CA8SIC) )  14.13,14 

14  REALC*REAL(C) 

IMAGC*AIHAG(C) 

IF  (REALC.LT.O.O)  GO  TO  17 
LOGC*CLOGtC) 

19  L0GD*AL0G(0) 

LOGA-LOGA+LOGC 

15  L0GB*L0GB+L06D 
LOGZ*LQGZ+LOGZRO 
LCGK*LOGK+LOGN 


LOGX=LOGA+LOGZ-LOGB-LOGK 
X=C£XPIL0GX>*U-1.0)**J) 
GO  TO  4 

16  REALA=-REALA 
IHAGA=-IMAGA 
A=CMPLX ( RE ALA  * I MAGA ) 
J=J*1 

SO  TO  18 

17  CREAL*-REALC 
IMAGC=-IMAGC 
CT=CMPLX(CREAL.INAGCI 
LOGC=CLOGCCT) 

J=J»1 
GO  TO  19 
13  X=0.0 
GO  TO  4 
END 
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COMPLEX  FUNCTION  GIS, OMEGA) 
COMPLEX  I0MEGA,Z1.P«Z2. 23, GAMMA, M 
COMPLEX  S 
8=0-5 

I OMEGA=CMPLX( O-Ot OMEGA) 

21  =  1.0'f  CCSQRT  (l.O+IOMEGA+S) ) 
P=S/(8.0*2l) 

Z2=-P 
GAMMA=-21 
Z  3= ( -GAMMA 1/2-0 
G=M( 22«B>23) 

RETURN 

END 
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COMPLEX  FUNCTION  PHI1< OMEGA) 

DOUBLE  PRECISION  D0MEGA,XI,X2.X3,X4,X5,X6»X7,U1,V1,U2.V2»U3, V3,X 
OOUBLE  PRECISION  TW0X.X8.X9, X10.X1 1» X13.U7.X14, V7,U8, V8, U4, V4.U9 
DOUBLE  PRECISION  V9, U10, V10.U11 » V11.U5, U6, V6.U12, V12.U13, V13 
DOUBLE  PRECISION  X12.V5 
DOME  GA= OBLE ( OMEGA ) 

X=OSQRT I DOMEGA 1/2.8284271247461 901 

TW0X=2.0*X 

U5=TM0X 

V5=TM0X 

X1=00MEGA**2 

X2=3.0*X1 

X3=4«0/X2 

X4=165.0/00NEGA 

X5=32«0/Xl 

X6=-55.0/C4.0*DOMEGA* 

X7*3.0/00MEGA 


Ul=1.0 

V1=X4 

U2=X3*U1 

V2=X34V1 

U3=l-0 

V?=X6 

X8-DEXPCTM0X) 
X9*'jcXP(— TMOX) 


■"Jr  i  JaaSsKt  ■ 


m 


X1Q* (X8-X9)/2.0 
Xll*tX8*X9)/2.0 
X12=DC0S(TW0X ) 

X13=X11+X12 
U7-X10/X13 
X14=0SIN(TW0X ) 

V7*X14/X13 

U8= ( <  U7*U5 ) ♦ ( V7* V5 ) I / C ( U5**2  J ♦ ( V5**2  > ) 
V8*I (U5*V7)-<  U7*V5) )/ti U5**2 )♦( V5**2 ) ) 
U4*X5*U3 
V4*X5*V3 

U9*(U4*U3)-(V4*V8) 

V9*(U4*V8)-MU8*V4) 

U10*1.0 

V10=X7 

U1X=(U7*U7)-(V7*V7> 

V11*(U7*V7)*( U7*V7) 

06=X3*U10 

V6*X3*V10 

UX2=(U64U11)-(V6*V11> 

vi2=(u6*vn)+(un*w6) 

U13*-U2-U9+U12 

V13=-V2-V9*V12 

U14=SNGL(U13) 

Vl4= SNGLI V13) 

PH11*CMPLX(U14,V14I 

RETURN 

ENO 


PHI2K  -  EFN  SOURCE  STATEMENT  -  IFN(S)  - 


COMPLEX  FUNCTION  PHI2(M, OMEGA) 

DOUBLE  PRECISION  DOMECr,Ul,VltU2.V2  ’3,V3tU4, V4.U5, V5,U7,V7.U8, VB 

OOUBLE  PRECISION  U9,V9,X,TM0)',X8,X9,xl0,Xll>X12,X13.X14,U10.V10 

DOUBLE  PRECISION  Ull.VlltU12'V12,U13,V13.U14.V14,U15,V15,U16.VX6 

DOUBLE  PRECISION  U17.V17,U18»V18tU19,V19 

DOHEGA=  DBLE ( ONEGA I 

U1*0.0 

Vl^l.O 

U2=UX/0NEGA 

V2=VX/ONEGA 

U3*0.5*U2 

W3=V2 

U4=U3*W 

V4=V3*W 

U5=-U4» ( I W**3 >/6.0) 

V5=-V4 

U7=OSQRT(OOMEGA) /I. 4142135623730950 

V7=U7 

U8=U7/2.0 

V8=V7/2.0 

U9-U84U 

V9=V8»H 

4*1 

X=U9 

3  TW0X=2.0*X 
X8=DEXP{TM0X) 

X9=DEXP (-TNQX ) 

X10=(X8-X9 1/2.0 
XllMX8*S9)/2.0 
X12=DC0SITH0X) 

X13*X11*X12 

XX4=DSIN(TM0X) 

GO  TO  (X,2),J 
X  UX0*XX0/X13 
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onooooooooononoooonnonoonoonnonnno 


V10*X14/X13 

J*2 

Ull= ( t U10*U7 1 *1 V10*V7 > ) / 1 1 U7**2  > ♦<  V7**2 )> 
V 1 1= ( C  U7*V10 ) -C  U10*V7 ) ) / ( ( U7*»2  )  ♦  <  V7 **2  J ) 
012=(U5*U11 )- ( V5*V11 ) 
V12=<U5*V11)*CUU*V5) 

U13=U1/(2.0*0MEGA> 

V13=V1/C2.0*0M£GA) 

U14=U13*< 1.0- IM**2) ) 

V14=V13*(1.0-IW**2> J 
U 15=. 33333333333 333333*U2 
V15*V2 
X*U8 
SO  ro  3 
2  U16=X10/X13 
V16=X14/XI3 

U17=l  (U16*U7)*(V16*V7)  J / U  U7  **2 ) ♦ I V7**2 ) ) 
VI 7= ( (U7*Vl6)-I U16*V7) )/ ( (U7**2 )♦( V7**2 ) > 
U18=IUl5*U17)-( V15*V17) 
V18=CU15*V17)*IU17*V15) 

019=012-014*018 
V 19= V12-V14*V 18 
U20=SNSL( U19) 

V2G*SN3L(V19) 

PHI2=CHPLX(U20,V20J 

RETURN 

END 
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-  EFN  SOURCE  STATEMENT  -  IFN(S)  - 


SUBROUTINE  QSF 


PURPOSE 

TO  COMPUTE  THE  VECTOR  OF  INTEGRAL  VALUES  FOR  A  GIVEN 
EQUIDISTANT  TABLE  OF  FUNC 'ION  VALUES. 


USAGE 

CALL  QSF  (H.Y.Z.NDIH) 


DESCRIPTION  OF  PARAMETERS 

H  -  THE  INCREMENT  OF  ARGUMENT  VALUES, 

r  -  THE  INPUT  VECTOR  OF  FUNCTION  VALUES. 

I  -  THE  RESULTING  VECTOR  OF  INTEGRAL  VALUES.  Z  MAT  Be 

IOENTICAL  WITH  Y. 

NOIM  -  THE  DIMENSION  OF  VECTORS  Y  AND  Z. 


REMARKS 

NO  ACTION  IN  CASE  NOIM  LESS  THAN  3. 


SUBROUTINES  ANO  FUNCTION  SUBPROGRAMS  REQUIREO 
NONE 


METHOD 

BEGINNING  WITH  Z(1I=0.  EVALUATION  OF  VECTOR  Z  IS  DONE  BY 
MEANS  OF  SIMPSONS  RULE  TOGETHER  WITH  NEWTONS  3/8  RULE  OR  A 
COMBINATION  OF  THESE  TWO  RULES.  TRUNCATION  ERROR  IS  OF 
ORDER  H**5  (I.E.  FOURTH  ORDER  METHOD).  ONLY  IN  CASE  NDIM=3 
TRU" CATION  ERROR  OF  Z(2)  IS  OF  ORDER  K<**4. 

FOR  REFERENCE.  SEE 

(1)  F.B.HILDE8RAND.  INTRODUCTION  10  NUMERICAL  ANALYSIS. 
MCGRAW-HILL.  NEW  YORK/TORONTO/LCSDON,  1956,  PP.71-76. 


,,.  ft  ZIXIMUEHL*  PRAKTISCHE  HATHENATIK  FUER  INGENltURE  UNO 
PHYS1 KER*  SPRINGER*  BERLtN/GQETTINGEN/HElOELBERG,  1963* 
PP. 216-221. 


SUBROUTINE  QSFlH*Y*Z»NPIN) 


DIMENSION  Y(l>»Zll> 

HT*.33i3333*H 
IFINDIH-5 >7,8,1 

NOIN  IS  GREATER  THAN  5.  PP.cPARATIONS  OF  INTEGRATION  LOOP 
1  SUN1*Y12)+Y(2) 

SU«1=SUN1-»SUM1 
SUH1=HT*( Y ( 1 > ♦SUHl*Y<  3) > 

AUXl*Y(4>*Y(4> 

JlU*l=AUXl+AUXl 

AUX1*SUH1+HT*(Y(3)*AUX1*-Yt5) )  . 

AL'X2»HT*tYU>+3.875*IY(2>«Y(5m2.625*(Y<3>*Yt4>>4Yl6>> 

SUH2=Y<S>-*Y(5> 

SUM2=$UN2*SUN2 

SUM2='AUX2-HT*(Y(6)+SUM2*Yt6)  ) 

zm*o. 

AUX=Yl3>*Yt3> 

AUX=AUX*AUX 

Z ( 2 > *SUN2-HT* ( Y ( 2 1 ♦AUX>Y 1 4> > 

Z ( 3) *SUNl 
Z  C4  >  — SUN2 
IFINDIN— 6)5»5*2 

INTEGRATION  LOOP 

2  00  6  I*7,N0IM,2 
SUN1-AUX1 
SUH2=AUX2 
AUXl=*YU-l»+YtI-l> 

AUX1*AUX1*AUX1 

AUX1=SUN1*HT*(YI I-2J*AUX1*YI I)> 

ZU-2I=SUN1 
IF( I-NO IN >3*6*6 

3  AUX2=YI I )*Y ( I ) 

AUX2=AUX2»AUX2 

A*IX?=SUN2»HT*!Y( I-1>4AUX2«7I 1*1 >  > 

6  Z ( I-1>=SUN2 

5  ZCNOIM-l)*AUXl 
Z(NDIN) =AUX2 
RETURN 

6  ZINDIN-l>*SUH2 
ZlNOIN>=AUXl 
RETURN 

END  OF  INTEGRATION  LOOP 

7  l FI  NOIN-3 >12* 11*8 

8  SUN2*1»125*HT*Iy!1)*YC2>«-VI2>*VI2>*YI3>«YC 3>*VI3>«-YI4>) 

SUNl=YI2>4YI2) 

SUNI=SUNUSUK1 
SUNl*HTMYtl>*SUNl+YI3>  > 

zin*o. 

AUX1*YI 3)*Y( 3 ) 

AUX1«AUX1*AUX1 

Z(2>*SUN2-HT*(YI2>»-AUXl«-YI4)  ) 

IF t NOIN-5) 10*9*9 
9  AUX1=YI4>*YI4> 

AUXl*AUXl»AUXl 
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036 
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Z(5)*SUNl+HT*(Y(3)+AUXl*Yt5) ) 

QSF 

103 

10 

zm=suni 

QSF 

104 

Z (4)*SUH2 

QSF 

105 

RETURN 

QSF 

106 

QSF 

107 

NOIN  IS  EQUAL  TO  3 

QSF 

108 

11 

SUHl=HT«‘U.254Y<l)*YI2)+YC2)-.25*Y(3n 

QSF 

109 

SUN2*Y ( 2 ) +Y ( 2 ) 

QSF 

110 

SUH2»SUM2+SUM2 

QSF 

111 

Z ( 3 ) *HT  *  t  Y { 1 ) +SUM2+ Y( 3)1 

QSF 

112 

Z(1)«0. 

QSF 

113 

Z (2)=SUM1 

QSF 

114 

12 

RETURN 

QSF 

115 

ENO 

QSF 

116 

£ 

I 
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APPENDIX  F.  SOLUTION  FOR  A  UNIFORM  RADIAL  PROFILE 

If  (36)  and  (37)  are  used  with  the  same  procedure 
as  that  described  in  Appendix  D  one  finds  that 


?■ 


-^-1 
2vhv0 


•j  ^  tanh(^p)  -z  tanh(az) 


and 


o9,= 


ReoRo /,  4z2\ 
~4  '  ' 


(F-i) 


(F-2) 


Equation  (F-l)  can  be  checked  by  comparing  its  limit  as  a 
approaches  zero  (let  tanh  x=x)  with  (F-2) . 

Now  consider  (56) : 


d2»n_ 

dz2 


(a2+^nV„ 


This  has  a  solution 


where 


(F-3> 


(F-4) 


The  boundary  conditions  given  by  (57)  require  that  K'" 
and 

*  2n-H 
«»■ — — »• 


-  o, 

; (f-5) 
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< 

*r 

m 


fV-TT 


The  coefficients  Cn  can  then  be  found  if  four  integrals 
(a,  b,  c,  d)  are  evaluated*  . 

Integral  (a)  is  given  by 
jj '*  [cosh(az)|[cos(fn2)Jd2=.  — — _  jjsinh  (az  )j  Jcos(Cnr  >j 


+  -gin__j^cosh  (a2)J  jsin(£nz)j 


• (p-10) 


o  2.671  #4 


Integral  (b)  is  given  by 


where 


Therefore 


f*  z  jsitvh (a i )j|cos(£n z)J  dz  =uz  j*-  JT$  udz  . 
>/[sinh  (az)][co5(en2)]  dz 

' aTfjl  [c«hto*)][co,(enz)]+-is_ [,inh(oz)][*in(fB  z)] 


[cosh(az)][cos(£nz}]dz 


(F-llj 


(F-12) 
2.671  #2 


(P-13) 


X2  [“>'■(«)][*»({„«)] 


(F-14) 
2.671  #1 


*  “  eValMti°n  “fies  «>e 
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From  (F-10)  through  (F-14)  one  finds  that 


j 2  i  jsinh  (oz)|  |cos  (£n*)]d* 

■?Sf*HT)]-^?['-(f)] 


(F-1S) 


Furthermore,  one  can  show  that  the  integral  (c)  is  given  by 


(F-1G) 


and  that  the  integral  (d)  is  given  by 


(F-17) 


(F-5)  has  been  used  in  all  of  these  evaluations. 

It  can  be  shown'  from  (F-9) ,  (F-10) ,  (F-l5) ,  (F-16) 
and  (F-17)  that 


4  Ro(-l)n*T  Qp  Cn  o2 

n  h  €n  (o*+tf)*  “2  +  C„* 


^sL  (o2+e2) 

2ReA  '  ' 


(F-18) 


To  find  G0  and  Gj^,  one  must  first  evaluate  the 
integral  in  (61) :  h 


1  *Jo  [w*] 


9 1  cosh(oz)  .  „  ^osh_(az__ 

r>  .  ,  .  '  Ri  .Ah. 


One  has  ®rellr=R|  |R|  cosh^-)  1  cosh^j 


h2  R* 


+  l*;SCne2R0°  Ro  ""  [cos(Cnz)J  floe' 


(F-X9) 


(F-20) 


>  ~.f- 


tiitiiu&jti&fi*  Vft\ 


E  ! 


I  ! 


% 


m 


^  i 


B- 


Just  as  in  appendix  D  the  linear  term  in  Ri  can  be  ignored 
and  the  exponential  term  in  the  series  set  equal  to  1.  Thus, 


®rel  r=R(  R ,1 


iut 


(F-21) 


Therefore,  from  (28),  (36),  (37),  and  (F-21),  (F-19)  becomes 


I  =- 


Qo 


2»hR* 


_ 1  f  2  a,  fcoshfaz)"]  dz 

cosh  0 


(F-22) 


^C./1  [C0Stfn2)]d2* 

Now  frOm  (F-l) 

ji  ,,  [cosh  (02)]*  «  i  ["■""(¥)]  f*  [“sh  <“»]* 

2  2  £sinh(az)j 


dz 


(F-23) 


Qpcosht  2  )  J 


4#wo2 


tanh 


+'s[,0"h(f)]-1 


Equation  (F-23)  can  be  checked  by  comparing  its 
limit  as  a  approaches  zero  with  the  result  obtained  by 
using  (F-2)  for  gi  in  the  integral  on  the  left  of  (F-23). 
One  finds  that 


r~  r~ 

JQ  oV^o-ojf  9i  [cosh(oz)jdz  = 


Qnh‘ 


(F-24 ) 


24™ 


where  the  relationships  tanh  x=x-  and  tanh4  x=x2  are 
used.  J 


If  (F-16) r  (F-18)  and  (F-23)  are  substituted  into 
(F-22) ,  one  finds  that 

•k\ 

dz* 


Qoh*  «  [lanhz( 

)+- 

ton 

4ww  olh2  L  '  2  ) 

'  a  h 

\  2  / 

- 

-2 


16 


»ha*uRp"j 

- 1 

o 

o 

1 

2R« 


t(*«) 


(F-25) 
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Equation  (63)  can  be  evaluated  from  (F-25) .  If  the  limit  of 
this  expression  as  a  approaches  zero  is  found  (with  the  help 
of  (P-24)),  (78)  can  be  determined.  Equation  (79)  can  be 
obtained  by  using  (F-15)  in  (60)  and  requiring  that  ki  be 
chosen  such  that  <J>=1  at  u=0. 


APPENDIX  G.  ERROR  DUE  TO  NEGLECTING 

cT  \  cr 

For  the  case  when  ki=k0=l,  (51)  becomes 


(G-l 

\ dr  r  / 

dr 

\  dr  T  /  d 2* 

If  (54)  and  (56)  are  used,  one  finds  that 


IF  kj 


(G-2 


The  coefficient  of  the  first  derivative  can  be  simplified 
if  one  defines  a  new  dependent  variable  (f3)n  such  t>iat 


f 


n 


r 


Ao 

2u  Cf3)n  . 


(G-: 


Then 


(^(^3)"  .  1  4  (f 3)r 


dr‘ 


dr 


(4=—)' 


0 


Finally,  let 


(G-4 


=  rv^-A 


(G-5 


where  (F-6)  has  been  used  to  evaluate  8n. 
Then  (G-4)  becomes 


d7»* 


+1? 


4(f»)n 

dij 


(M„  =  0 


(GH 


The  solution  of  this  equation  is  given23  by 


Equations  (H-l) ,  (H-3)  and  (H-4)  cam  be  combined  to  obtain 


1  .1  (t^z)=  ± 

1  dr  \  dr  /  ft  d* 


(H-5) 


If  this  equation  is  integrated  twice  and  the  boundary 
conditions  are  imposed,  one  obtains 


(H-6) 


where  A  is  determined  by  the  total  flow  rate  through  the 
drain.  Since  there  are  two  drains  in  this  model,  one  finds 
that 


(H-7) 


where  Qd  is  the  volumetric  flow  rate  specified  by  the  total 
mass  flow  rate  through  the  sensor  and  the  static  pressure 
in  either  drain. 


Now  let 


vas(v.+v*)a  o«lw,> 


(H-8) 


where  Vx  is  the  solution  that  satisfies  the  inhomogeneous 
boundary  conditions  at  the  drain  wall  and  V2  is  the  solu¬ 
tion  that  is  used  to  satisfy  the  inhomogeneous  boundary  con¬ 
ditions  on  the  upstream  and  downstream  boundaries  of  the 
drain  region.  Clearly,  Vi  is  a  function  of  r  only.  Thus, 
from  (H-2) ,  (H-7)  and  (H-8)  one  finds  that 


(H-9) 


with  the  boundary  conditions 


VI(R,)  =  R, 


(H~10) 


V,  (01*0, 


(H-ll) 


The  solutions  of  (H-9)  are  the  modified  Bessel  functions.  3  If 
the  boundary  conditions  are  applied,  it  is  found  that 


it 


15 


1 


W 


I 


£  (wJn(r)e 


<nnz 


(H-14) 


where  Real  (Yi)n  and  Real  (Y2)n  >0* 

The  first  series  on  the  right  side  of  this  equation 
describes  how  information  introduced  at  the  upstream  bounding 
surface  decays  as  the  fluid  flows  down  the  drain.  The  second 
series  describes  the  propagation  in  the  upstream  direction 
from  the  downstream  boundary,  and  will  be  ignored  since  it 
can  be  shown  that  the  propagation  in  the  upstream  direction 
is  of  extremely  short  range. 

In  the  actual  model  tested,  the  pickoff  was  located 
just  outside  the  drain  exit.  For  convenience  the  pickoff 
is  treated  here  as  if  it  were  located  inside  the  drain  tube 
so  that  upstream  propagation  can  be  ignored.  This  latter 
model  does  not  properly  account  for  changes  in  output  brought 
on  by  changes  in  Mach  number  near  the  actual  pickoff,  for  it 
does  not  allow  the  streamline  pattern  to  change  with  m.  How¬ 
ever,  the  viscous  losses  predicted  by  such  a  model  should  be 
fairly  accurate  since  the  losses  due  to  upstream  propagation 
will  be  small  if  the  pickoff  is  submerged  in  the  exhaust  jet 
at  a  point  where  the  spread  would  not  be  great  if  the  pickoff 
were  not  there. 

In  the  actual  model,  the  pickoff  holes  were  less 
than  a  drain  diameter  downstream  from  the  exit.  Thus,  the 
pickoff  can  be  considered  as  being  in  the  drain. 

If  (H-13)  and  (H-14)  are  used  and  the  second  series 
in  (H-14)  is  ignored,  then  one  has 
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k 


Finally,  let 


\  =  ~  2  c, 


Then  (H-20)  becomes 


d2t*4>n  = 

1,2  d*,2  d72  4^4)n“  * 


This  is  the  confluent  hypergeometric  equation, 
is  given22  by 


(w4)n  =  a'm  .  o  .  v2  ] 

+  B'i,2  M  f(l-  -|i)  ,  2  ,  ,f2J  • 


(H-23) 


(H-24) 


The  solution 


(H-25) 


If  (H-16)  ,  (H-18)  ,  (H-19) ,  (H-21)  and  (H-23)  are 
used,  (H-25)  becomes 


Cf* 


(H-26) 


Since  Vg  i  be  bounded  at  r  =  0,  A  =  0.  It  can  then  be 
shown  that 

_in  I2.  r  .  , 


»«-  a'w‘  2 


'  4S^Sn  1 6 Ref  4^  nlRi/ 


(H-27) 


where  6n's  are  eigenvalues  that  satisfy  the  boundary  condition 


M[(l  +  — ^r-  + — -  ~r).  2  ,  8nl» 
_V  4S/  Sn  16 Re,2  4  /’  " 


(H-28) 


Equation  (85)  can  be  found  by  using  equation  (H— 8) , 
(H-12) ,  (H-14),  and  (H-27),  and  by  noting  that 


(Yl)n  3Re7  Ri 


(H-29) 
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